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1 . 0  Research  Description 


The  main  objective  of  this  research  is  to  carry  out  theoretical  and  experimental 
investigation  on  optical  wave  mixing  and  phase  conjugation  in  various  nonlinear  media, 
including  photorefractive  crystals  and  optical  fibers,  as  well  as  novel  applications  using  these 
nonlinear  optical  effects.  During  the  research  period,  we  have  carried  out  intensive  theoretical  and 
experimental  investigation  in  several  important  areas  in  wave  mixing  and  phase  conjugation, 
including  the  problem  of  multigrating  effects  in  MPPC,  fanning  noise  reduction  by  using 
partially  coherent  beams,  grating  enhancement  and  restoration  in  photorefractive  media,  fidelity 
of  image  restoration  in  multimode  fibers,  reflection-based  phase  conjugators,  cross-talk  noise  in 
volume  holographic  storage,  optical  phase  conjugation  and  pulse  restoration  in  single  mode 
fibers,  photorefractive  spatial  mode  converter  for  multimode-to-single-mode  fiber-optic  coupling, 
energy  coupling  by  partially  degenerate  four-wave  mixing  in  multichannel  lightwave  Systems, 
form  birefringence  of  layered  media  and  volume  gratings,  the  effect  of  beam  coherence  on 
mutually  pumped  phase  conjugators  (MPPC),  polarization-dependent  mechanism  transformation 
during  self-pumped  phase  conjugation  in  Ce:BaTi03,  intensity-dependent  absorption  and 
photorefractive  properties  of  some  new  crystals,  including  Rb:KNb03,  and  CerBaTiOs. 

1.1  Scientific  Problem 

The  unique  distortion  correction  property  of  optical  phase  conjugation  continues  to  be 
extremely  attractive  in  many  advanced  applications.  These  include  frequency  conversion  and 
pulse  restoration  in  fiber  optical  communication,  wavefront  correction  in  free  space  optical 
communication,  high  power  lasers,  optical  image  processing,  etc.  There  are  many  approaches 
which  can  be  employed  to  achieve  optical  phase  conjugation.  These  include  optical  four-wave 
mixing,  stimulated  Brillouin  scattering  (SBS),  self-pumped  phase  conjugation,  etc.  In  nonlinear 
media  such  as  photorefractive  crystals  and  Kerr  materials,  it  is  possible  to  achieve  mutually 
pumped  phase  conjugation  (MPPC).  Such  phase  conjugators  are  ideal  for  optical 
communications  through  scattering  media  such  as  atmosphere  or  sea  water.  MPPC  can  also  be 
employed  for  the  phase  locking  of  an  array  of  lasers.  It  is  known  that  both  SPPC  and  MPPC  czin 
involve  several  index  gratings  (reflection  gratings,  transmission  gratings  and  2k  gratings)  which 
may  compete  for  charge  carriers  resulting  in  possible  instability  in  the  phase  conjugation  process. 
In  addition,  the  coherence  of  the  beams  plays  an  important  role  in  the  strength  of  the  gratings.  In 
the  area  of  fiber  optical  communication,  it  is  desirable  to  have  frequency  conversion  devices 
which  can  provide  a  frequency  translation  for  the  signal  carrier  in  WDM  systems.  Using  optical 
four-wave  mixing  it  is  possible  to  achieve  the  frequency  translation  in  optical  fibers  via  the  Kerr 
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effect.  An  important  issue  is  the  energy  efficiency  and  the  response  time  of  the  system  (latency). 
Using  optical  four-wave  mixing,  it  is  also  possible  to  achieve  forward  phase  conjugation  which 
can  be  utilized  for  pulse  restoration  in  high  speed  optical  communication  system. 

1 . 2  Scientific  and  T echnical  Approach: 

Mutually-pumped  phase  conjugation  (MPPC)  is  a  nonlinear  optical  process  in  which  two 
incoherent  laser  beams  mutually  pump  each  other  and  generate  two  phase  conjugate  beams 
simultaneously.  Most  theoretical  models  for  MPPC  consider  only  transmission  gratings.  In 
practice,  when  the  laser  beams  become  coherent  (e.g.  laser  phase-locking)  or  the  nonlinear 
material  has  fast  response  time  (e.g.  Kerr  media),  the  formation  of  reflection  gratings  and/or  2k 
gratings  is  inevitable.  Although  there  have  been  several  models  for  MPPC,  the  effect  of  reflection 
and  2k  gratings  in  MPPC  has  never  been  investigated.  In  an  effort  to  understand  the  dynamic 
competition  between  transmission  gratings,  reflection  gratings  and  2k  gratings,  we  theoretically 
investigated  the  coupling  between  all  the  waves  by  using  a  set  of  multi-wave  coupled-wave 
equations.  To  vary  the  strength  of  reflection  gratings  and  2k  gratings,  we  used  a  laser  with  a 
finite  coherence  length.  By  varying  the  path  difference,  it  is  possible  to  control  the  mutual 
coherence  between  the  coupled  beams.  Experimental  investigations  using  BaTi03  crystals  have 
also  been  carried  out  to  validate  our  theory.  In  our  investigation  into  the  possibility  of  optical 
pulse  restoration  using  phase  conjugation,  we  theoretically  analyze  the  coupling  between  optical 
waves  in  an  optical  fiber  via  the  Kerr  effect.  By  solving  the  coupled  wave  equations,  we  obtained 
an  exact  solution  for  the  intensity  of  all  the  four  waves  involved.  The  solution  is  useful  to 
provide  the  energy  conversion  efficiency  as  well  as  the  phase  conjugation  efficiency. 

Optical  four-wave  mixing  in  nonlinear  media  such  as  optical  fibers  can  be  employed 
canbe  employed  for  the  compression  of  optical  pulses  as  well  as  the  restoration  of  pulse  shape 
due  to  group  velocity  dispersion.  Holograms  and  gratings  are  responsible  for  wave  mixing  and 
phase  conjugation  in  nonlinear  media,  including  photorefractive  crystals.  The  strength  of  these 
holograms  and  gratings  depends  on  the  coherence  of  the  writing  beams.  In  SPPC  and  MPPC, 
the  writing  beams  may  depend  on  the  crystal  orientation  and  the  interaction  configuration.  During 
this  period,  we  have  investigated  theoretically  and  experimentally  various  models  of  SPPC  in 
Ce:BaTi03  crystals.  To  investigate  the  effect  of  2k-gratings,  we  used  an  erasing  beam  to  weaken 
the  2k-gratings  in  SPPC.  The  specific  model  for  SPPC  also  depends  on  the  fanning  pattern.  By 
using  various  polarization  states,  we  are  able  to  vary  the  fanning  patterns  leading  to  different 
models  of  operation.  We  also  investigate  the  dispersion  of  form  birefringence  due  to  layered 
structures. 
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1.3  Publications: 


1.3.1  Papers  Published: 


283.  R.  F.  Nabiev,  P.  Yeh,  and  D.  Botez,  "Spatial  gap  solitons  in  periodic  nonlinear 
structures,"  Opt.  Lett.,  IS,  1612-1614  (1993). 

284.  S.  M.  Zhou,  Q.  B.  He,  P.  Yeh,  and  H.-K.  Liu,  "Spatial  fidelity  in  photorefractive  image 
amplification,"  Opt.  Comm,  99,  18-24  (1993). 

285.  Q.  B.  He,  P.  Yeh,  L.  J.  Hu,  S.  P.  Lin,  T.  S.  Yeh,  S.  L.  Tu,  S.  J.  Yang,  "Shift-invariant 
photorefractive  joint-transform  correlator  using  Fe:LiNb03  crystal  plates,"  Appl.  Opt., 
32,  3113-3115  (1993). 

293.  Y.  Zhu,  D.  P.  Zhang,  P.  Yeh,  and  C.  Gu,  "High  Efficiency  Photorefractive  BaTiOs 
Plate  for  Optical  Image  Processing,"  Chinese  Science  Bulletin,  3Si  1424-1427  (1993). 

298.  Q.  B.  He  and  P.  Yeh,  "Fanning  noise  reduction  in  photorefractive  amplifiers  using 
incoherent  erasures,"  Appl.  Opt.,  33,  283-287  (1994). 

304.  C.  Gu  and  P.  Yeh,  "Partial  Phase  Conjugation,  Fidelity,  and  Reciprocity,"  Opt. 
Commun.,  107,  353-357  (1994). 

310.  Y.  Zhang,  S.  Campbell,  P.  Yeh,  D.  Shen,  X.  Ma  and  J.  Chen,  "Submillisecond 
photorefractive  response  time  of  KNb03:Rb''',"  Opt.  Lett.,  12,  1397-1399  (1994). 

311.  X.  Yi,  P  Yeh  and  C.  Gu,  "Statistical  analysis  of  cross-talk  noise  and  storage  capacity  in 
volume  holographic  memory,"  Opt.  Lett.,  19,  1580-1582  (1994). 

314.  P.  Yeh  and  A.  Chiou,  "Optical  Phase  Conjugation  for  Interconnection  and  Image 
Processing,"  Chapter  12  of  Real-Time  Optical  Information  Processing,  Ed.  B.  Javidi,  p. 
475-531  (Academic  Press,  Inc.,  1994). 

315.  Q.  B.  He,  S.  Campbell,  P.  Yeh,  X.  Y.  Ma  and  D.  Shen,  "Compact  Reflection-grating- 
based  phase  conjugator  in  Fe:  IG4b03  crystals,"  Appl.  Opt.  33, 4320-4322  (1994). 

318.  Q.  B.  He,  H.  K.  Liu  and  P.  Yeh,  "Asymmetric  Photorefractive  Fabry-Perot  Etalons," 
Appl.  Phys.  B-Lasers  and  Optics,  52, 467-470  (1994). 

319.  W.  Wu,  P.  Yeh  and  S.  Chi,  "Phase  Conjugation  by  Four-Wave  Mixing  in  Single-Mode 
Fibers,"  Photonics  Tech.  Lett.,  1448-1450  (1994). 

324.  P.  Yeh  and  C.  Gu,  "Photorefractive  Nonlinear  Optics  and  Applications,"  Series  in 
Nonlinear  Optics  and  Optical  Physics,  Vol.2,  Chapter  12, 341-388  (1994). 


328.  C.  Yang,  Y.  Zhang,  P.  Yeh,  Y.  Zhu  and  X.  Wu,  "Photorefractive  properties  of 
Ce:BaTi03  crystals, "  Opt.  Comm.,  113.  416-420  (1995). 
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329. 


S.  Campbell,  P.  Yeh,  C.  Gu  and  Q.  B.  He,  "Fidelity  of  image  restoration  by  partial 
phase  conjugation  through  multimode  fibers, "  Opt.  Comm.,  1 14. 50-56  (1995). 

330.  Q.  B.  He  and  P.  Yeh,  "Photorefractive  mutually  pumped  phase  conjugation  with  partially 
coherent  beams,"  Appl.  Phys.  B  60, 47-50  (1995). 

334.  S.  Campbell,  P.  Yeh,  C.  Gu,  S.  H.  Lin,  C-J  Cheng,  and  K.  Y  Hsu,  "Optical  Restoration 
Of  Photorefractive  Holograms  Through  Self-Enhanced  Diffraction,"  Opt.  Lett.,  20,  330- 
332  (1995). 

336.  X.  Yi,  S.  Campbell,  P.  Yeh,  and  C.  Gu,  "Statistical  Analysis  Of  Cross-Talk  Noise  And 
Storage  Capacity  In  Volume  Holographic  Memory  -  Image  Plane  Holograms,"  Opt. 
Lett.,  20  779-781  (1995). 

338.  A.  E.  T.  Chiou,  P.  Yeh,  C.  Yang  and  C.  Gu,  "Photorefractive  spatial  mode  converter  for 
multimode-to-single-mode  fiber-optic  coupling,"  Opt.  Lett.,  2Q,  1 125-1127  (1995). 

341.  W.  Wu  and  P.  Yeh, "  Energy  Coupling  By  Partially  Degenerate  Four-Wave  Mixing  In 
Multichannel  Lightwave  Systems,"  IEEE  Photonics  Technology  Letters,  Vol.  7,  585- 
587  (1995). 

346.  A.  Chiou,  P.  Yeh,  C.  Yang,  C.  Gu,  "Photorefractive  coupler  for  fault-tolerant  coupling," 
IEEE  Photonics  Tech.  Lett.,  7, 789-791  (1995). 

348.  C.  Gu  and  P.  Yeh,  "Form  birefriogence  of  layered  media  and  volume  gratings,"  J.  Opt. 
Soc.  Am.,  B,  12,  1094-1099  (1995). 

35,1.  Y.  Lian,  S.  Lin,  S.  Campbell,  K.  Y  Hsu,  et  al,  "Polarization-Dependent  Mechanism 
Transformation  During  Self-Pumped  Phase  Conjugation  In  Ce:BaTi03,"  Optics  Lett.  20, 
1683-1685  (1995). 

352.  P.  Yeh,  C.  Gu,  C-J.  Cheng,  K.  Y  Hsu,  "Hologram  Enhancement  In  Photorefractive 
Media,"  Optical  Engineering,  34. 2204-2212  (1995). 

357.  S.C.  Delacruz,  S.  Marcorma,  J.  Feinberg,  Q.  B.  He,  H.  K.  Liu,  and  P.  Yeh,  "Effect  of 
Beam  Coherence  on  Mutually  Pumped  Phase  Conjugators,"  J.  Opt.  Soc.  Am.,  B12, 
1363-1369  (1995). 

358.  .X.  Yi,  P.  Yeh,  and  C.  Gu,"Cross-Talk  Noise  In  Volume  Holographic  Memory  With 

Spherical  Reference  Beams,"  Optics  Lett.,  20, 1812-1814  (1995). 

359.  C.  Yang,  Y.  H.  Zhang,  X.  Yi  X,  P,  Yeh,  et  al,  "Intensity-Dependent  Absorption  and 
Photorefractive  Properties  in  Cerium-Doped  Batio3  Crystals,"  Journal  Of  Applied 
Physics,  2S,  4323-4330  (1995). 

361.  P.  Yeh,  C.  Gu,  C.  J.  Cheng  and  K.  Y.  Hsu,  "Optical  Restoration  of  Photorefractive 
Holograms,"  Applied  Physics  B-Lasers  and  Optics,  M  ,511-514  (1995). 

366.  C.  Gu  and  P.  Yeh, "Form  birefringence  dispersion  in  periodic  layered  media,"  Opt.  Lett., 
11,  504-506  (1996). 

380.  X.  Yi,  S.  H.  Lin,  P.  Yeh  and  K.  Y.  Hsu,  "Contradirectional  Two-Wave  Mixing  With 
Partially  Coherent  Waves  In  Photorefractive  Crystals,"  Opt.  Lett.,  Vol.  21  No.  15, 1 123- 
1125  (Aug  1,  1996). 
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381.  S.  H.  Lin,  Y.  Lian,  P.  Yeh,  K.  Y.  Hsu  et  al,  "2k-Grating-Assisted  Self-Pumped  Phase 
Conjugation  -  Theoretical  And  Experimental  Studies,"  J.  Opt.  Soc.  Am.  B-Optical 
Physics,  Vol.  13,  No.  8,  1772-1779,  (Aug  1996). 

395.  C.  Yang  and  P.  Yeh,  "Form  Birefringence  Of  Volume  Gratings  In  Photopolymers," 
Applied  Physics  Letters,  V69  N23:3468-3470  (Dec.  2, 1996) 

398.  C.  Yang  and  P.  Yeh,  "Artificial  Uniaxial  And  Biaxial  Dielectrics  With  The  Use  Of 
Gratings,"  J.  Appl.  Phys.,  V81  Nl:23-29  (Jan.  1997). 


1.3.2  Conference  Papers; 

264.  R.  F.  Nabiev,  P.  Yeh  and  D.  Botez,  "Spatial  gap  soliton  in  periodic  nonlinear 
structures,"  in  Integrated  Photonics  Research  Conference  Technical  Digest,  1993  (Optical 
Society  of  America,  Washington,  DC,  1993),  Vol.  10,  pp.  312-315. 

271.  Q.  B.  He,  P.  Yeh,  X.  Ma  and  D.  Shen,  "Compact  self-pumped  phase  conjugator  in 
photorefractive  FerKNbOs,"  Conference  on  Lasers  and  Electro-Optics,  1993  Vol.  11, 
OSA  Technical  Digest  Series  (Optical  Society  of  America,  Washington,  DC  1993),  pp. 
526-527. 

272.  A.  Chiou,  P.  Yeh,  and  C.  Gu,  "Spatial  mode  conversion  in  photorefractive  fibers,"  in 
OSA  Annual  Meeting  Technical  Digest,  1993  (Optical  Society  Of  America,  Washington, 
DC  1993),  Vol.  16,  p.  7. 

278.  Q.  He,  X.  Yi,  P.  Yeh,  S.  Chi,  "Optical  pulse  compression  by  means  of  a  single  mode 
fiber  and  phase  conjugation,"  in  OSA  Annual  Meeting  Technical  Digest,  1993  (Optical 
Society  of  America,  Washington,  DC  1993),  Vol.  16,  p.  182. 

290.  A.  Chiou,  P.  Yeh,  C.  Gu,  and  R.  Neurgaonkar,  "Beam  Coupling  and  Spatial  Mode 
Conversion  in  Photorefractive  Planar  Waveguide,"  Conference  Proceedings  of 
lEEE/LEOS  Annual  Meeting,  319-320  (San  Jose,  California,  1993). 

299.  S.  Campbell,  P,  Yeh,  Q.  B.  He,  and  C.  Gu,  "Studies  of  image  fidelity  during  partial 
phase  conjugation,"  Conference  on  Lasers  and  Electro-Optics,  Vol.  8,  1994  OSA 
Technical  Digest  Series  (Optical  Society  of  America,  Washington,  DC  1994),  pp.  12-13. 

300.  S.  Campbell,  P.  Yeh,  and  Daofan  Zhang,  "Photorefractive  and  thermal  fixing  properties 
of  4  mol%-doped  lithium  niobate,"  Conference  on  Lasers  and  Electro-Optics,  Vol.  8, 
1994  OSA  Technical  Digest  Series  (Optical  Society  of  America,  Washington,  DC  1994), 

p.  80. 

301.  Q.  B.  He,  H.-K.  Liu  and  P.  Yeh,  "Reflection  gratings  in  the  photorefractive  mutually 
pumped  phase  conjugators,"  Conference  on  Lasers  and  Electro-Optics,  Vol.  8,  1994 
OSA  Technical  Digest  Series  (Optical  Society  of  America,  Washington,  DC  1994),  p. 
97. 

305.  Y.  Zhang,  S.  Campbell,  P.  Yeh;  D.  Shen,  X.  Ma,  J.  Chen,  "Grating  response  time  of 
photorefractive  KNb03:Rb+,"  Proc.  1994  IEEE  Conference  on  Nonlinear  Optics, 
Materials,  Fundamentals,  and  Applications  (Waikoloa,  Hawaii,  July  25-29,  1994),  pp. 
284-286. 
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306.  W.  Wu,  P.  Yeh,  S.  Chi,  "Frequency  conversion  by  four-wave  mixing  in  single-mode 
fibers,"  Proc.  1994  IEEE  Conference  on  Nonlinear  Optics,  Materials,  Fundamentals,  and 
Applications  (Waikoloa,  Hawaii,  July  25-29,  1994),  pp.  332-334. 

307.  P.  Yeh,  C.  Gu,  C.-J.  Cheng,  K.  Y.  Hsu,  "Hologram  restoration  and  enhancement  in 
photorefractive  media,"  (Invited  Paper)  Proc.  1994  IEEE  Conference  on  Nonlinear 
Optics,  Materials,  Fundamentals,  and  Applications  (Waikoloa,  Hawaii,  July  25-29, 
1994),  pp.  421-423. 

308.  S.  Campbell,  P.  Yeh,  C.  Gu,  S.  H.  Lin,  C.  -J.  Cheng,  K.  Y.  Hsu,  "Optical  self¬ 
enhancement  of  photorefractive  holograms,"  Proc.  1994  IEEE  Conference  on  Nonlinear 
Optics,  Materials,  Fundamentals,  and  Applications  (Waikoloa,  Hawaii,  July  25-29, 
1994),  pp.  430-432. 

309.  X.  Yi,  P.  Yeh,  C.  Gu,  "Cross-talk  noise  and  storage  density  in  holographic  memory," 
Proc.  1994  IEEE  Conference  on  Nonlinear  Optics,  Materials,  Fundamentals,  and 
Applications  (Waikoloa,  Hawaii,  July  25-29,  1994),  pp.  436-438. 
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The  following  is  a  list  of  significant  progresses  we  have  made  during  this  period  under 
the  support  of  this  program: 

•  Photorefractive  MPPC  with  partially  coherent  beams 

•  Fanning  noise  reduction  in  photorefractive  amplifiers  using  incoherent  erasures 

•  Compact  Reflection-grating-based  phase  conjugator  in  Fe:  KNb03  crystals 

•  Asymmetric  Photorefractive  Fabry-Perot  Etalons 

•  Submillisecond  photorefractive  response  time  of  KNb03:Rb+ 

•  Statistical  analysis  of  cross-talk  noise  and  storage  capacity  in  volume  holographic 
memory  -both  image  plane  and  Fourier  plane 

•  Photorefractive  Properties  of  Ce:BaTi03 

•  Phase  Conjugation  by  Four-Wave  Mixing  in  Single-Mode  Fibers 

•  Image  restoration  by  partial  phase  conjugation  through  multimode  fibers 

•  Partially  Degenerate  Four-Wave  Mixing  in  Single-Mode  Fibers 

•  Form  birefringence  in  layered  media 

•  Polarization-Dependent  Mechanism  Transformation  During  Self-Pumped  Phase 
Conjugation  In  Ce:BaTi03 

•  Hologram  Enhancement  and  restoration  In  Photorefractive  Media 

•  Effect  of  Beam  Coherence  on  Mutually  Pumped  Phase  Conjugators 

•  Cross-Talk  Noise  In  Volume  Holographic  Memory  With  Spherical  Reference  Beams 

•  Intensity-Dependent  Absorption  and  Photorefractive  Properties  in  Cerium-Doped 
Batio3  Crystals 

•  Form  birefringence  dispersion  in  periodic  layered  media 


2.2  Progress  details: 

Since  most  of  the  details  of  research  results  have  been  reported  earlier  in  the  annual 
reports,  we  only  describe  the  details  of  the  research  results  during  the  period  between  the  last 
annual  report  and  the  final  report.  Detail  discussion  on  some  of  the  important  research  results  for 
the  whole  period  can  be  found  in  Section  3  which  contains  reprints  of  these  papers. 

Self-focusing  effects  in  periodic  leaky-wave  coupled  (antiguided)  and  evanescent-wave 
coupled  array  structures  are  analyzed  for  what  is  to  our  knowledge  the  first  time.  The  results 
obtained  can  be  explained  in  terms  of  spatial  gap  solitons.  For  infinite-extent  arrays  of  antiguides 
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it  is  found  that  (for  interelement-spacing  values),  below  the  inphase-mode  resonance,  self- 
focusing  occurs  with  increasing  light  intensity  even  when  the  medium  is  intrinsically  self- 
defocusing  (i.e.,  negative  nonlinearity).  If  the  medium  has  positive  nonlinearity  the  analysis 
shows  that  antiguided  arrays  possess  either  self-focusing  or  self-de focusing  properties 
depending  on  the  parameters  of  the  structure,  whereas  evanescent-wave  in-phase  coupled  arrays 
always  have  self-focusing  properties.  (Opt.  Lett.,  18,  1612-1614  (1993)) 

The  spatial  fidelity  of  amplified  images  in  photorefractive  materials  due  to  spatial  gain 
variations  with  computation  and  experiment  based  on  BaTiOi  is  investigated.  The  results  indicate 
how  the  spatial  fidelity  of  the  amplified  images  decreases  as  the  relative  aperture  of  the  Fourier 
transform  lens  used  in  the  system  increases,  and  varies  with  the  angle  between  the  pump  beam 
and  the  axis  of  the  signal  converging  wave.  (Opt.  Comm,  99, 18-24  (1993)) 

The  authors  report  the  results  of  their  experimental  investigation  on  a  shift-invariant 
photorefractive  image  correlator  that  uses  a  thin  crystal  plate  of  Fe:LiNbO/sub  3/,  which  operates 
in  the  Raman-Nath  regime  of  diffraction.  (Appl.  Opt.,  22,  3113-3115  (1993)) 

We  propose  and  demonstrate  a  new  scheme  of  beam  fanning  reduction  in  photorefractive 
amplifiers  by  using  incoherent  erasure.  This  simple  and  effective  approach  employs  a  pump 
beam  with  a  mixed  polarization  state  or  introduces  a  third  incoherent  beam  for  selective  erasure. 
The  approach  is  analyzed,  and  the  experimental  results  are  presented  and  discussed.  (Appl.  Opt., 
33,  283-287  (1994)) 

The  fidelity  of  partially  phase  conjugated  images  is  considered  taking  into  account  the 
finite  aperture  of  phase  conjugate  mirrors.  The  analysis  is  based  on  the  conservation  of  energy 
and  the  reciprocity  theorems.  In  the  absence  of  absorption,  the  conservation  of  energy  leads  to  a 
simple  relationship  between  the  image-to-background  ratio  (I/B)  and  the  aperture  of  the  phase 
conjugator.  In  the  presence  of  absorption,  the  reciprocity  theorem  leads  to  an  expression  of  the 
image-to-background  ratio  (I/B)  that  depends  on  the  fraction  of  power  intercepted  by  the  phase 
conjugator.  The  results  are  presented  and  discussed.  (Opt.  Commun.,  107.  353-357  (1994)) 

We  investigate  the  photorefractive  response  time  of  diffraction  gratings  in  Rb/sup  +/- 
doped  KNbO/sub  3/  crystals.  Two  components  of  charge  carriers,  with  different  response  times, 
are  found  to  contribute  to  the  photoinduced  index  gratings.  One  component,  which  contributes  to 
approximately  70%  of  the  index  grating,  shows  a  submillisecond  rise  time  when  the  total 
intensity  of  the  writing  beams  is  larger  than  5  W/cm/sup2/.  The  dependence  of  the  response  time 
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on  the  light  intensity  and  the  grating  wave  vector  is  also  investigated.  (Opt.  Lett.,  1^,  1397-1399 
(1994)) 


We  present  a  statistical  analysis  of  cross-talk  noise  and  consider  the  cross-talk-limited 
storage  capacity  of  volume  holographic  memory.  Random-walk  approximation  is  employed  to 
estimate  the  intrapage  interpixel  cross-talk  noise  and  the  signal-to-noise  ratio.  We  obtain  simple 
expressions  of  the  storage  capacity  in  terms  of  number  of  bits  per  volume  for  a  given  signal-to- 
noise  ratio  or  bit-error  rate.  The  results  indicate  that  the  storage  density  is  reduced  from  the 
ultimate  density  of  lambda/sup  -3/  by  a  factor  related  to  the  signal-to-noise  ratio.  (Opt.  Lett.,  19, 
1580-1582(1994)) 

A  brief  review  of  phase  conjugators  that  use  photorefractive  materials  as  the  nonlinear 
media  is  presented.  An  introduction  to  phase  conjugation  is  followed  by  a  description  of  the 
photorefractive  effect  and  wave  mixing  in  these  media.  Various  applications  using 
photorefractive  phase  conjugators,  especially  optical  interconnection  and  real-time  imaging 
processing  are  discussed.  (Chapter  12  of  Real-Time  Optical  Information  Processing,  Ed.  B. 
Javidi,  p.  475-531  (Academic  Press,  Inc.,  1994)) 

We  present  the  design  and  the  results  of  experimental  investigations  of  a  self-pumped 
phase  conjugator  through  a  reflection  grating  in  a  special-cut  photorefractive  KNb03:Fe  crystal. 
High  reflectivity,  good  phase-conjugate  fidelity,  and  fast  response  time  are  obtained.  (Appl.  Opt. 
33,  4320-4322  (1994)) 

We  theoretically  analyze  the  properties  of  Fabry-Perot  etalons  which  consist  of 
photorefractive-crystal  plates  with  asymmetric  mirror  reflectance.  Peak  transmissions  and  non¬ 
reciprocities  are  enhanced  significantly  compared  with  the  symmetric  photorefractive  etalons.  The 
effect  of  material  absorption  is  also  considered.  (Appl.  Phys.  B-Lasers  and  Optics,  52, 467-470 
(1994)) 


We  present  an  exact  solution  to  the  coupled-mode  equations,  governing  four-wave 
mixing  in  single-mode  fibers,  under  the  perfect  phase-matching  condition.  The  solution  accounts 
for  pump  depletion  and  fiber  absorption.  Optimum  length  of  fiber  is  derived  and  phase 
conjugation  efficiency  is  calculated  and  discussed.  The  calculated  results  are  in  good  agreement 
with  the  reported  experimental  results.  (Photonics  Tech.  Lett.,  5, 1448-1450  (1994)) 

We  investigate  the  absorption  and  photorefractive  properties  of  Ce:BaTi03  (cerium  doped 
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barium  titanate)  crystals.  The  absorption  coefficient,  electro-optic  gain  coefficient,  and  total 
effective  trap  density  at  wavelength  5 14.5  nm  increase  systematically  with  cerium  doping 
concentration.  The  absorption  and  the  electro-optic  gain  coefficients  exhibit  a  dependence  on  the 
optical  intensity.  This  indicates  that  there  are  at  least  two  active  photorefractive  species.  One  of 
these  crystals  has  the  highest  beam-coupling  gain  coefficient  reported  to  date,  gamma  =8.06  cm" 
1  for  ordinarily  polarized  light.  The  response  time  of  as-grown  Ce:BaTi03  with  15  ppm  cerium 
is  220  ms  at  an  intensity  of  800  mW/cm^.  (Opt.  Comm.,  113. 416-420  (1995)). 

We  investigate  the  fidelity  of  image  restoration  under  conditions  of  partial  phase 
conjugation  through  a  random  medium.  For  the  first  time,  the  conjugated  image-to-background 
ratio  is  experimentally  measured  as  an  evaluation  of  the  phase  conjugate  mirror's  fidelity.  Output 
signal-to-noise  ratio  and  information  throughput  are  also  considered  as  functions  of  aperture 
limitations  at  the  various  planes  of  the  optical  system.  (Opt.  Comm.,  114, 50-56  (1995)) 

We  investigate  the  mutually  pumped  phase-conjugation  process  in  photorefractive  media 
with  partially  coherent  pump  beams.  We  consider  the  effects  of  transmission,  reflection,  and  2k- 
reflection  gratings  as  well  as  the  depletion  of  both  pump  beams.  The  theoretical  results  can  be 
employed  to  qualitatively  explain  the  experimental  observations.  (Appl.  Phys.  B  60,  47-50 
(1995)) 


For  the  first  time  to  our  knowledge,  we  demonstrate  the  use  of  self-enhanced  diffraction 
to  optically  restore  weak  dynamic  volume  holograms.  This  restoration  process  is  achieved 
through  the  use  of  a  novel  beam-toggling  scheme  that  significantly  enhances  the  steady-state 
diffraction  efficiency  of  volume  index  gratings  stored  in  photorefractive  media.  Experiments 
demonstrating  and  confirming  theory  are  presented  and  discussed.  (Opt.  Lett.,  2^,  330-332 
(1995)) 

We  investigate  the  cross-talk  noise  in  optical  storage  based  on  angle-multiplexed  image 
plane  volume  holograms.  Simple  expressions  for  the  signal-to-noise  ratio  and  the  storage  density 
are  obtained.  The  cross-talk  noise  is  found  to  limit  the  size  of  the  pixels  and  the  number  of 
recorded  holograms.  The  cross-talk-limited  storage  density  of  image  plane  holographic  storage  is 
found  to  be  close  to  that  of  Fourier  plane  holographic  storage.  (Opt.  Lett.,  20  779-781  (1995)) 

We  report  what  is  to  our  knowledge  the  first  experimental  demonstration  of  a 
photorefractive  spatial  mode  converter  (based  on  mutually  pumped  phase  conjugation)  that 
couples  light  efficiently  from  a  multimode  fiber  into  a  single-mode  fiber  with  an  extremely  large 
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degree  of  tolerance  to  misalignment.  Using  an  argon  laser  (514.5  nm)  and  a  barium  titanate 
crystal,  we  have  demonstrated  that  the  laser  light  can  be  coupled  from  a  multimode  fiber  (core 
diameter  100  mu  m,  numerical  aperture  0.37)  into  a  single-mode  fiber  (core  diameter  2.9  mu  m, 
numerical  aperture  0.1 1),  with  an  efficiency  of  approximately  15%  and  an  alignment  tolerance  of 
approximately  100  mu  m.  The  coupling  efficiency  is  more  than  two  orders  of  magnitude,  and  the 
tolerance  to  misalignments  is  more  than  30  times  better  than  the  corresponding  values  achievable 
by  conventional  techniques.  (Opt.  Lett.,  20, 1125-1127  (1995)) 

We  investigate  the  energy  coupling  among  different  frequency  carriers  by  partially 
degenerate  four-wave  mixing  in  single-mode  fibers.  We  show  that  for  a  multichannel  lightwave 
system  operated  near  the  zero-dispersion  frequency  of  the  single-mode  fiber,  four-wave  mixing 
may  cause  significant  signal  depletion  in  a  channel  whose  power  acts  as  a  pump  in  the  four-wave 
mixing  process.  Analytical  solutions  are  derived  and  the  signal  depletion  due  to  four-wave 
mixing  is  obtained.  (IEEE  Photonics  Technology  Letters,  Vol.  7,  585-587  (1995)) 

We  report  the  first  experimental  demonstration  of  a  fault-tolerant  laser-to-fiber  coupler 
that  is  insensitive  to  both  angular  and  lateral  misalignments.  FauU-tolerance  is  achieved  by 
sending  a  laser  beam  into  a  photorefractive  crystal  that  automatically  senses  and  tracks  the 
location  of  the  fiber  facet,  and  diffracts  the  input  beam  (via  self-generated  adaptive  dynamic 
holograms)  into  the  fiber.  Using  an  argon  laser  (514.5  nm),  a  barium  titanate  crystal,  and  a 
multimode  fiber  (core  diameter=400  um  m,  NA=0.16),  we  have  achieved  a  coupling  efficiency 
of  58%  and  a  tolerance  of  more  than  0.5  mm  of  lateral  misalignment.  Techniques  to  improve  the 
coupling  efficiency  and  the  degree  of  fault  tolerance,  and  to  extend  this  method  for  the  coupling 
of  a  laser  diode  output  into  a  single-mode  fiber  are  discussed.  (IEEE  Photonics  Tech.  Lett.,  7, 
789-791  (1995)) 

We  analyze  the  form  birefringence  dispersion  of  periodic  dielectric  layered  media.  We 
investigate  the  shape  of  the  normal  surface  and  its  dependence  on  wavelength.  Approximation 
gives  an  explicit  expression  for  the  additional  dispersion  of  the  form  birefringence  that  is  due  to 
the  finite  layer  thickness.  (J.  Opt.  Soc.  Am.,  B,  12,  1094-1099  (1995)) 

We  report  the  results  of  experimental  investigations  of  the  mechanism  transformation 
from  four-wave  mixing  and  stimulated  photorefractive  backscattering  to  four-wave  mixing  and 
total  internal  reflection  during  self-pumped  phase  conjugation  in  a  photorefractive  crystal  by 
varying  the  plane  of  polarization  of  the  input  beam.  The  use  of  a  mixed-polarization  input  beam 
can  change  the  fanning  pattern,  which  leads  to  the  mechanism  transformation  and  an 
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enhancement  of  the  phase-conjugation  reflectivity.  The  temporal  dynamics  of  the  buildup  process 
indicate  that  such  an  enhancement  is  due  to  the  presence  of  both  types  of  mechanism  during 
phase  conjugation.  In  addition,  we  compare  the  contributions  of  2k  gratings  with  the  phase 
conjugation  during  both  mechanisms.  (Optics  Lett.  20,  1683-1685  (1995)) 

We  consider  the  possibility  of  restoration  and/or  enhancement  of  decaying  holograms  in 
photorefractive  media  by  using  a  simple  optical  readout  in  conjunction  with  a  phase  conjugator. 
The  results  indicate  that  extremely  weak  holograms  can  be  enhanced  provided  that  the  two-beam 
coupling  is  sufficiently  strong.  Steady-state  photorefractive  holograms  can  be  maintained 
continuously  without  decay  by  using  a  properly  designed  readout  scheme.  The  result  also 
provides  an  explanation  for  the  formiation  of  mutually  pumped  phase  conjugation  in  terms  of  the 
amplification  of  initial  noise  gratings.  (Optical  Engineering,  M.  2204-2212  (1995)) 

The  performance  of  four  different  mutually  pumped  phase  conjugators  in  barium  titanate 
depends  on  the  mutual  coherence  of  the  two  input  beams.  In  three  of  the  conjugators  the  use  of 
fully  mutually  coherent  input  beams  enhances  the  overall  phase-conjugate  reflectivity.  The 
authors  find  that  back-scattering  and  transmission  gratings  in  the  photorefractive  crystal  both 
contribute  to  the  phase-conjugate  signal  and  that  their  relative  strength  depends  on  the  relative 
coherence  of  the  two  input  beams.  We  numerically  solve  coupled-wave  equations  that  include  all 
gratings  and  find  reasonable  agreement  between  our  theory  and  our  experimental  data.  (J.  Opt. 
Soc.  Am.,  B12,  1363-1369  (1995)) 

We  investigate  angle-multiplexed  volume  holographic  memory  with  spherical  reference 
beams,  for  which  the  spherical  approximation  is  made  to  model  the  wave-front  distortion  in 
general.  We  find  that  the  angular  selectivity  and  the  cross-talk  noise  with  spherical  reference 
beams  are  close  to  those  with  planar  reference  beams.  The  results  indicate  that  angle-multiplexed 
volume  holographic  memory  can  be  realized  in  compact  systems  for  which  large  wave-front 
distortion  is  expected.  (Optics  Lett.,  20,  1812-1814  (1995)) 

We  report  the  optical  absorption  and  photorefractive  properties  of  a  series  of  cerium- 
doped  barium  titanate  crystals.  The  light-induced  absorption  and  the  intensity-dependent  two- 
beam  coupling  gain  coefficient  indicate  that  there  are  at  least  two  active  trap  levels  in  cerium- 
doped  BaTiOs  crystals.  The  total  effective  trap  density,  the  Debye  screening  wave  vector,  and 
the  intensity-dependent  factor  increase  with  cerium-doping  concentration  as  well  as  the  pump 
beam  intensity.  The  positive  light-induced  absorption  and  the  negative  absorption  gain  coefficient 
as  the  grating  wave  vector  approaches  zero  indicate  that  the  excitation  cross  section  of  shallow 
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traps  is  larger  than  that  of  deep  traps.  One  of  these  crystals  exhibits  a  two-beam  coupling 
coefficient  of  8.06  cm’l  for  ordinarily  polarized  beams,  which  is  the  highest  ever  reported  to 
date.  Both  electro-optic  index  grating  and  absorption  grating  are  contributing  to  the  coupling 
coefficient.  The  response  time  of  as-grown  cerium-doped  BaTiOs,  with  15  ppm  cerium  is  220 
ms  at  an  intensityof  800  mW/cm^.  (Journal  Of  Applied  Physics,  78, 4323-4330  ( 1995)) 

We  consider  the  possibility  of  restoration  and/or  enhancement  of  decaying  holograms  in 
photorefractive  media  by  using  a  simple  optical  readout  in  conjunction  with  a  phase  conjugator. 
The  results  indicate  that  extremely  weak  holograms  can  be  enhanced  provided  that  the  two-beam 
coupling  is  sufficiently  strong.  Steady-state  photorefractive  holograms  can  be  maintained 
continuously  without  decay  by  using  a  properly  designed  readout  scheme.  The  result  also 
provides  an  explanation  for  the  formation  of  mutually  pumped  phase  conjugation  in  terms  of  the 
amplification  of  initial  noise  gratings.  (Applied  Physics  B-Lasers  and  Optics,  ^  ,511-514 
(1995)) 


We  analyze  the  form  birefringence  dispersion  of  periodic  dielectric  layered  media.  We 
investigate  the  shape  of  the  normal  surface  and  its  dependence  on  wavelength.  Approximation 
gives  an  explicit  expression  for  the  additional  dispersion  of  the  form  birefringence  that  is  due  to 
the 

finite  layer  thickness.  (Opt.  Lett.,  21,  504-506  (1996)) 

We  investigate  contradirectional  two-wave  mixing  with  partially  coherent  waves  in 
photorefractive  crystals  in  the  nondepleted  pump  regime.  Equations  governing  the  propagation  of 
the  self-coherence  function  and  the  mutual-coherence  function  of  the  signal  wave  and  the  pump 
wave  are  derived  and  simulated  numerically.  Numerical  solutions  of  these  equations  are  in 
excellent  agreement  with  the  experimental  measurements.  (Opt.  Lett.,  Vol.  21  No.  15,  1123- 
1125  (Aug  1, 1996)) 

We  investigated  2k-grating-assisted  self-pumped  phase  conjugation  in  photorefractive 
Ce:BaTiO/sub  3/  crystals.  The  phase-conjugation  process  involves  a  combination  of  four-wave 
mixing  and  stimulated  photorefractive  backscattering.  An  approximation  involving  separate 
interaction  regions  is  used  to  theoretically  calculate  the  reflectivity  of  phase  conjugation  as  a 
function  of  the  coupling  strength  of  four-wave  mixing  and  stimulated  photorefractive 
backscattering.  In  our  experiments,  grating-erasure  techniques  are  employed  at  the  interaction 
regions  to  investigate  the  dependence  of  phase-conjugate  reflectivity  on  the  coupling  strength  of 
four-wave  mixing  and  stimulated  photorefractive  backscattering.  The  experimental  results  are  in 
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good  agreement  with  the  theoretical  prediction.  (J.  Opt.  Soc.  Am.  B-Optical  Physics,  Vol.  13, 
No.  8,  1772-1779,  (Aug  1996)) 

We  investigate  experimentally  the  form  birefringence  due  to  photoinduced  one-  and  two- 
dimensional  volume  index  gratings  in  DuPont  photopolymers.  Experimental  results  on  the 
dispersion  of  the  form  birefringence  versus  the  illumination  wavelength  and  the  grating  period 
for  one-dimensional  gratings  are  presented.  We  also  investigate  the  form  birefringence  due  to 
two-dimensional  volume  gratings.  Positive  form  birefringence  and  biaxial  form  birefringence  can 
be  obtained  by  recording  two  volume  gratings  with  different  grating  wave  vectors.  (Applied 
Physics  Letters,  V69  N23:3468-3470  (Dec.  2, 1996)) 

We  investigate  the  form  birefringence  due  to  photoinduced  one-  and  two-dimensional 
volume  index  gratings  in  DuPont  photopolymers.  The  effective  medium  theory  has  been  applied 
to  the  continuous  volume  subwavelength  gratings.  With  the  presence  of  a  one-dimensional 
volume  subwavelength  grating,  the  originally  isotropic  photopolymer  exhibits  a  negative  form 
birefringence  and  behaves  like  a  uniaxial  dielectric.  A  form  birefringence  of  -1.2*10‘3  at 
>,=632.8  nm  is  measured  at  a  grating  period  of  A  =447  nm.  Experimental  results  on  the 
dispersion  of  the  form  birefringence  versus  the  illumination  wavelength  and  the  grating  period 
are  presented.  With  the  presence  of  a  two-dimensional  grating,  in  general,  the  originally  isotropic 
photopolymer  exhibits  a  biaxial  form  birefringence  and  behaves  like  a  biaxial  dielectric.  Positive 
form  birefringence  can  also  be  obtained  by  recording  two  volume  gratings  with  different  grating 
wave  vectors.  By  recording  two  photoinduced  volume  gratings  we  can  synthesize  biaxial 
dielectrics  with  the  principal  axes  along  any  desired  directions.  We  also  report  the  measured 
principal  refractive  indices  of  the  medium  due  to  the  presence  of  two-dimensional  volume 
gratings  in  DuPont  photopolymer.  (J.  Appl.  Phys.,  V81  Nl:23-29  (Jan.  1997)) 


Two- wave  mixing  in  photorefractive  crystals  is  a  very  useful  nonlinear  optical 
phenomenon  in  many  applications,  such  as  image  amplification,  laser  beam  cleanup,  spatial  light 
modulators,  thresholding  and  power  limiting  devices.  Most  of  the  theoretical  works  in  this  area 
are  on  the  two  wave  interaction  with  monochromatic  beams.  Two  wave  mixing  with  partially 
coherent  beams  has  been  studied  recently  for  the  case  of  transmission  grating  interaction.  In  our 
publications  (e.g.  Optics  Letters,  1  Aug.  1996,  vol.21,  (no.  15):  1123-5),  we  present  the 
theoretical  analysis  and  experimental  investigation  for  the  case  of  reflection  grating  interaction. 

In  the  case  of  contradirectional  two  wave  mixing  with  purely  diffusive  photorefractive 
medium,  the  coupled-wave  equations  for  the  slowly  varying  amplitudes  E,(t,z)  and  £2(1,2)  can 
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be  written  as 


aE,  ^  1  8E,  ^  Y  Q  E, 
dz  V  at  2 1,  +  Ij  2  ‘ 

aE,  laE,  _yQ-  E,  ^ 
dz  V  at  2  I|  + 12  2  ' 


where  y  is  the  intensity  coupling  constant,  a  is  the  intensity  absorption  coefficient,  v  is  the 
group  velocity,  Eiexp(-ia)t  +  ikz)  and  E2exp(-io)t  -  ikz)  are  the  coupled  quasi-monochromatic 
waves  which  interact  through  a  dynamic  photorefractive  grating  5n(z,t)  «=  Qexp[2ikz]  +  c.c.. 
The  dynamics  of  the  photorefractive  grating  can  be  written  as 

t^  +  Q  =  E,E,-  (3) 

where  T  is  the  relaxation  time  constant  determined  by  the  average  total  intensity, 

r~((E,y-<-(E!)T 

Consider  that  the  temporal  behavior  of  each  wave’s  complex  amplitude  is  a  stationary  * 
random  process,  with  coherentce  time  6co‘‘  being  substantially  less  than  relaxation  time  of  the 
material  8q)i:»  1.  When  the  optical  path  difference  of  the  two  beams  is  smaller  than  the 
coherence  length  of  the  laser,  a  dynamic  photorefractive  grating  is  recorded  in  the  medium.  Its 
position  and  profile  are  nearly  temporally  constant.  Then,  as  an  approximatioin,  we  can  replace 
the  dynamic  grating  amplitude  Q  in  Eqs.(l)  and  (2)  with  its  ensemble  average  as  Q  «  (EjE2*). 

Since  the  complex  amplitudes  E,(t,z)  and  £2(1,  z)  are  stationary  random  processes,  we 
can  define  some  of  their  ensemble  averages  as  r,2(z,At) s ^E,(z,t,)E2*{z,t2)), 
r,,(z,At)  =  (E,(z,t,)E,*{z,t2)),  and  r22(z,At)  =  (E2(z,t,)E2*(z,t2))  where  At  =  t,-t2.  We 
refer  to  T,,  and  Tjj  as  the  self  coherence  of  the  pump  beam  Ei(t,z)  and  the  signal  beam 
E2(t,z)  respectively,  and  r,2  as  the  mutual  coherence  between  the  signal  beam  and  the  pump 
beam.  With  these  notations,  we  can  immediately  write  down  Q  =  r,2(z,0),  I,(z)  =  r,,(z,0)  and 
l2(z)  =  £22(2,0).  Eqs.(l)  and  (2)  can  therefore  be  reduced  to  a  system  for  these  average  values 
as 

__2 3ri2(z,At)  Y ri2(z,o)r  /^AO+r  (zAt)l  (4) 

dz  V  dAt  +2  I,  +  l2  ^  ’  ^  ’ 
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These  solutions  are  consistent  with  the  classical  theory  on  optical  wave  interference  in  linear 
medium.  The  mutual  coherence  function  is  presented  as  a  propagation  wave  with  the  boundary 
conditions  at  z=o  plane.  The  self  coherence  decrease  along  the  z  axis  due  to  absorption.  When 
there  is.  coupling,  the  photorefractive  grating  implies  that  the  both  beam  are  rescattered  by  the 
grating  into  each  other.  As  a  result,  part  of  pump  wave  branches  off  in  the  direction  of  the  signal 
waves.  It  causes  that  the  mutual  and  self  conherence  functions  of  the  two  waves  are  not  the  same 
as  simple  propagation  of  the  boundary  condition.  However,  we  note  that  Eqs.(4)-(6)  are  self- 
consistent.  They  can  be  solved  numerically  if  we  have  the  proper  boundary  conditions.  In  the 
presence  of  pump  depletion,  we  can  not  obtain  sufficient  boundary  conditions  for  typical 
situations.  Therefore,  although  Eqs.(4)-(6)  are  general,  we  can  not  use  them  to  solve  the  case  of 
pump  depletion.  For  the  case  of  nondepleted  pump,  sufficient  boundary  conditions  can  be  easily 
obtained. 

Consider  a  two  wave  mixing  configuration  shown  in  Fig.l.  Both  the  signal  beam  and  the 
pump  beam  come  from  the  same  laser  source.  If  the  laser  source  is  assumed  to  have  a  Gaussian 
lineshape  (fit  from  the  experimental  data)  with  a  linewidth  of  Av,  then  the  normalized  self 
coherence  of  the  laser  source  can  be  written  as 
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At  the  signal  beam  input  plane  z=0,  the  intensity  ratio  between  the  signal  beam  and  the  pump 
beam  is  assumed  to  be  p,  and  the  time  delay  between  the  signal  beam  and  the  pump  beam  is 
assumed  to  be  5t.  In  our  simulation,  the  intensity  of  the  pump  beam  is  taken  to  be  unit  one. 
Then  the  boundary  conditions  at  z=0  can  be  written  as 

r,2(z  =  0,At)  =  Vpr,(At  +  6t)  (11) 

r„(z  =  o,At)sr,(At)  (12) 

r22(z  =  o,At)spr,(At)  (13) 

With  these  boundary  conditions,  Eqs.(4)-(6)  become  an  initial  value  problem,  tinder  the 
approximation  of  nondepleted  pump. 

The  results  of  our  calculations  are  show  in  Figs.  2,  in  which  we  have  plotted  the  mutual 
coherence  function  of  the  two  beams  and  the  self  coherence  function  of  the  signal  beam  as  a 
function  of  distance  for  various  coupling  constant  y  of  1,  3,  5,  and  10.  The  initial  optical  path 
difference  is  taken  to  be  8  cm  at  the  signal  incident  plane  z=0,  the  line  width  of  the  laser  Av  is 
assumed  to  be  1.83  GHz,  and  the  input  intensity  ratio  p  is  assumed  to  be  lO^.  For  a  small 
coupling  constant,  the  mutual  coherence  function  is  hardly  affedted  by  the  presence  of  the 
coupling.  The  curve  is  almost  the  same  as  the  propagation  of  wave  equation  using  Eqs.(7)  and 
(11).  For  a  lager  coupling  constant,  as  a  result  of  the  coupling,  part  of  the  pump  wave  branches 
off  in  the  direction  of  the  signal  beam.  When  the  weak  beam  is  amplified,  it  gains  substantially 
by  the  instantaneous  convergence  of  the  tempered  profile  of  the  strong  pump  beam,  increasing  the 
strong  beam’s  coherence  with  weak  beam.  In  other,  the  increase  of  the  mutual  coherence  of  the 
two  waves  also  enhances  the  amplitude  of  the  grating.  The  interaction  length  of  the  two  beams 
could  therefore  become  longer  than  the  coherent  legth  of  the  laser  source.  It  leads  also  to  a  sharp 
increase  of  the  gain  of  the  signal  beam  when  the  coupling  constant  is  large  (y  >  5,  as  shown  in 
Fig.  2(b)). 

Our  experimental  arrangement  for  this  study  is  shown  in  Fig.l.  Two  partially  coherent 
beams,  obtained  by  splitting  an  argon  laser  beam  with  linewidth  of  1.83  Ghz,  passed  through 
the  signal  beam  and  the  pump  brenches  of  the  system  respectively  and  then  were 
contradirectionally  incident  into  a  KNb03:Co  crystal  (y  «  3.3cm‘‘ ,  a  *=  0.5cm"‘  and  thickness 
L=0.72  cm).  The  optical  path  difference  of  the  two  beams  at  the  signal  incident  plane  (i.e.  z=0) 
was  set  to  be  AL  =  ^2  -  To  monitor  the  mutual  coherence  between  the  signal  beam  and  the 

pump  beam  at  the  output  plane  z=L,  another  reference  beam  ( E2ref )  was  splitted  from  the  pump 
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beam  £2-  The  optical  path  difference  of  Hi  and  E2ref  beams  was  adjusted  to  be  the  same  as  that  of 
El  and  £2  beams  at  the  output  plane  z=L.  The  visibility  of  interference  fringes  generated  by  Ei 
and  E2ref  beams  was  observed  by  a  CCD  camera  at  the  output  plane  PI.  Therefore  the  mutual 
conherence  r,2  can  be  etsimated  in  terms  of  the  visibility  as  (l^  -  Imi„)/(4-x/Ii  I2)’  where 
(Imax  “  Imin)  ^be  amplitude  of  the  fringes.  In  the  experiment  we  monitor  the  interference  pattern 

with  and  without  pump  beam  £2-  Fig.3  shows  pictures  with  an  initial  mutual  coherence 
rj2(0,0)«0.43  (i.e.  AL  =  4  cm)  and  intensity  ratio  P  =  664.  Notice  the  increase  of  fringes 
visibility  due  to  the  coupling.  The  measured  mutual  coherence  r,2(L,0)  increases  from  0.19  to 

0.7.  To  futher  verify  the  theoretical  prediction  of  the  increase  of  the  mutual  coherence,  we  then 
measure  the  mutual  coherence  r,2(L,0)  as  a  function  of  the  artifical  optical  path  difference  AL. 

Fig.  4  shows  a  complete  family  of  data  obtained  in  this  measurement  (dots)  and  the  theoretical 
curve  for  the  sjune  parameters.  Owing  to  the  reflected  wave  by  the  grating  being  completely 
coherent  with  the  pump  wave  at  the  reflected  face,  it  leads  to  an  significant  increase  in  the  mutual 
coherence.  The  curve  shows  this  effect  which  makes  the  mutual  conherence  increase  very 
sharply  as  the  optical  path  difference  decreases. 

Another  effect  on  the  increase  of  the  mutual  coherence  is  the  enhancement  of  the 
amplitude  of  the  grating  and  interaction  length  of  the  two  beams  inside  crystal.  This  enhancement 
influences  the  achievable  intensity  gain  of  the  signal  wave.  Fig.  5  shows  the  measurement  of  the 
two-wave-mixing  gain  (dots)  of  the  signal  wave  at  z=L  plane  as  a  function  of  the  path  difference 
AL.  Along  with  the  datas  is  the  theoretical  curve  for  the  same  parameters.  The  curve  shows 
asymmetric  with  respect  to  the  maximum  point  of  the  intensity  gain  and  the  shift  of  the  point 
inside  crystal  due  to  the  different  tails  of  the  decrease  of  the  coherence  function.  The  experimental 
results  are  excellent  agreement  with  the  theoretical  prediction. 

In  conclusion  we  have  demonstrated  that  the  contradirectional  two  wave  mixing  of  partial 
coherent  beams  inside  a  photrefractive  crystal  with  an  inertial  nonlinear  response  leads  to  a 
increase  of  the  mutual  coherence  subatantially  and  the  interaction  legth  of  the  two  beams.  A  set  of 
the  coupled  equation  has  been  derived  to  describe  the  propagation  of  the  mutual  coherence  and 
the  self  coherence  functions  of  the  two  beams.  The  experimental  results  have  demonstrated  the 
increase  of  the  mutual  coherence  and  the  enhancement  of  the  two-wave-mixing  gain  of  the  signal 
beam  due  to  the  coupling. 

Figures: 
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Fig.  1.  Schematic  of  the  two-wave-mixing  configuration  used  in  our  calculations  and 
experiments.  The  distance  and  £2  optical  path  length  of  the  signal  beam  and 

pump  beam  from  laser  source  to  the  signal  beam  incident  plane  z=0,  respectively.  The 
£2tef  is  the  optical  path  length  of  reference  beam  E2ref  from  laser  source  to  the  signal 

output  plane  z=L. 
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Fig.2.  (a)  The  mutual  coherence  function  r,2(z,0)  and  (b)  the  self  coherence  function  of  the 
signal  beam  ri,(z,0)  as  a  function  of  z  for  various  coupling  constant  y  of  1,  3,  5,  and 
10. 
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Fig.  3.  Interference  pattern  of  the  signal  beam  and  the  reference  beam  at  the  output  plane  PI  (a) 
without  (b)  with  coupling.  Notice  the  increase  of  fringe  visibility  due  to  the  coupling. 
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Fig.  5.  The  measurement  of  the  two-wave-mixing  gain  of  the  signal  wave  at  z=L  plane  as  a 
function  of  the  path  difference  AL  of  the  two  beams. 
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Energy  Coupling  by  Partially  Degenerate  Four- Wave 
Mixing  in  Multichannel  Lightwave  Systems 

Weishu  Wu  and  Pochi  Yeh,  Fellow.  IEEE 


Abstract — We  investigate  the  energy  coupling  among  different 
frequency  carriers  by  partially  degenerate  four-wave  mixing 
in  single-mode  fibers.  We  show  that  for  a  multichannel  light¬ 
wave  system  operated  near  the  zero-dispersion  frequency  of  the 
single-mode  fiber,  four-wave  mixing  may  cause  significant  signal 
depletion  in  a  channel  whose  power  acts  as  a  pump  in  the  four- 
wave  mixing  process.  Analytical  solutions  are  derived  and  the 
signal  depletion  due  to  four-wave  mixing  is  obtained. 

FOUR- WAVE  MIXING  (FWM)  is  an  important  nonlinear 
process  for  optical  fiber  communication  systems  [1H3]. 
Applications  of  four-wave  mixing  such  as  frequency  conver¬ 
sion  [4]  and  phase  conjugation  [5],  [6]  have  been  analyzed 
and  demonstrated.  On  the  other  hand,  FWM  can  degrade  the 
performance  of  a  multichannel  lightwave  system  [7].  One 
problem  related  to  FWM  is  the  cross  talk  resulting  from 
the  generation  of  new  frequency  components  for  channels 
supposed  to  be  zero.  Another  problem  associated  with  FWM 
is  the  signal  depletion  due  to  energy  coupling  between  dif¬ 
ferent  channels.  Although  the  effects  of  cross  talk  can  be,  to 
some  extent,  alleviated  by  designing  a  system  composed  of 
unequally  spaced  frequency  channels  [7],  signal  depletion  is 
inevitable  for  both  systems  with  equal  and  unequal  frequency 
spacings  if  the  input  powers  are  sufficiently  large. 

FWM  in  single-mode  fibers  has  been  extensively  studied 
under  the  assumption  of  pump  nondepletion  [l]-[4],  [8]. 
Analytical  solutions  taking  into  account  of  pump  depletion 
has  been  obtained  only  for  cases  where  fiber  absorption  [2], 
[9]  and/or  group  velocity  dispersion  [8]  can  be  neglected. 
However,  these  assumptions  are  not  valid  for  cases  in  which 
long  interaction  length  of  the  fiber  is  utilized  to  obtain  efficient 
FWM  [5]-[7].  In  this  paper,  we  present  analytical  solu¬ 
tions  to  partially  degenerate  four-wave  mixing  (PDFWM)  in 
single-mode  fibers  near  the  zero-dispersion  frequency.  These 
solutions  are  then  used  to  analyze  the  energy  coupling  in  both 
systems  with  unequal  and  equal  frequency  spacings. 

To  describe  the  energy  coupling,  we  start  from  the  coupled¬ 
mode  equations  for  PDFWM, 

^  =  -i27.4IA2.<43exp(tAA:z)  - 
az  l 

exp{-iAkz)  -  ^^2 
dz  2 
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—  =  ^i^AiAiAl  exp{-iAkz)  -  -A3  (1) 

where  a  is  the  absorption  coefficient  of  the  fiber,  7  is  the 
nonlinearity  coefficient  [1]  of  the  fiber,  =  2^7  -  /:>  - 
is  the  propagation  constant  difference  describing  the  phase 
mismatch  of  the  PDFWM  process,  and  A  1,^2,  and  A3  are 
the  amplitudes  for  the  center  frequency  component  and  two 
side  components,  respectively.  In  arriving  at  (1),  we  assumed 
that  the  intensity-dependent  phase  shift  is  small  enough  to  be 
neglected.  The  phase  matching  condition  for  PDFWM,  Ait  = 
0,  can  be  satisfied  if  the  central  frequency,  /^  =  (/2  -i-  /3)/2, 
is  set  to  the  zero-dispersion  frequency,  /o,  as  well  as  all  the 
frequencies  are  near  fo  so  that  the  higher  order  dispersion 
terms  can  be  neglected  [10].  Under  the  phase  matching  con¬ 
dition,  the  above  equations  can  be  rewritten,  by  substituting 
-  \/^exp(i(^j)  =  y/Q~exp{-az/2)exp{i(pj),  as 


dz 

dQ2 

dz 

=+27  \/QiQi<32Q3  exp(-Q;z)sin  Ao 

dA<i 


=-47v/Qi(3iQ2Q3exp(-ar)sin  A<? 
=  +27\/QiQiQ2Q3 exp(-Q:2)sin  A<p 


dz 


■  exp(-ar)  cos  A0, 


(2) 


f3) 


where  Ap  =  2pi  -  02  -  03  is  the  relative  phase  of  the  three 
waves. 

It  can  be  easily  seen  that  the  direction  of  energy  coupling 
is  determined  by  the  sign  of  the  relative  phase,  Ap.  In  cases 
when  ^2(0)  or  P^iO)  is  absent,  the  FWM  process  will  result 
in  a  relative  phase  such  that  energy  will  be  transferred  to  P2 
or  P3.  If  the  energy  is  coupled  from  Pi  to  P2  and  P3,  the 
solution  to  (2)  can  be  written  as 


Qi=Pmy 

Q2=0.5Pi(0)(l~y)-f  P2(0) 
Q3=0.5Pi(0)(l-y)  +  P3(0)  (4) 

where  y  is  a  function  of  z  to  be  determined. 

Analytical  expression  for  y  has  been  derived  for  initial 
condition  PaCO)  =  0  [6].  When  P^iO)  =  0,  it  has  been  shown 
that  the  relative  phase  is  always  7r/2  due  to  the  fact  that 
v/CjiOiOzOa  cos  A0  is  a  constant.  In  that  case,  y  is  given  by 

^  1  /CN 

y  = - - — --  -  (5) 

rcosh''(\/rTT/)  -h  1 

where  r  =  2P2(0)/Pi(0)  and  /  =  7Pi(0)[l  -  exp{—az)]/a. 
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Fig.  1.  Power  coupling  by  PDFWM  as  a  function  of  fiber  length  r.  The  dashed  lines  are  the  results  obtained  using  the  analytical  solution  without  the 
effects  of  SPM/XPM.  while  the  solid  lines  are  the  results  obtained  using  numerical  methods  taking  into  account  the  effect  of  SPM/XPM  for  input  powers 
of  (a)  16  dBm  and  (b)  18  dBm.  The  nonlinear  coefficient  -  is  1.853  x  lO”^  cm”^  W-^.  and  fiber  absorpuon  coefficient  is  0.23  dB/km  for  a  fiber  with 
effective  area  of  70  /im*  and  nonlinear  refractive  index  n2  of  3.2  x  10”^®  m^/W  at  A  =  1.55  mm. 


Fig.  1  shows  the  energy  coupling  by  PDFWM  in  a  50  km 
long  fiber.  The  dashed  lines  are  the  results  obtained  using  the 
analytical  solution  derived  above,  while  the  solid  lines  are  the 
results  obtained  using  numerical  methods  taking  into  account 
the  phase  shift  caused  by  the  intensity-dependent  refractive 
index.  The  input  powers  are  16  dBm  and  18  dBm  in  Fig. 
1(a)  and  (b),  respectively.  It  shows  that  for  moderate  input 
powers,  the  effect  of  the  intensity-dependent  phase  shift  can 
be  reasonably  neglected. 

It  is  worth  noticing  that  although  (5)  was  derived  under 
the  boundary  condition  P3(0)  =  0,  it  is  still  valid  for  cases 
where  a  small  PaCO)  is  present.  Although  the  presence  of  PaCO) 
may  result  in  an  arbitrary  value  for  the  initial  relative  phase 
the  relative  phase  will  quickly  approach  7r/2  before  any 
significant  energy  coupling  occurs  because  is  very 

large  (roughly  proportional  to  When  PaCO)  can  not 

be  ignored  compared  with  the  other  two  inputs,  which  is  the 
case  for  multichannel  lightwave  systems  with  equal  frequency 
spacings  and  with  all  input  powers  present,  the  energy  coupling 
depends  on  A<p,  which  determines  the  interference  between 
the  input  signals  and  the  newly  generated  FWM  components. 
By  inspecting  (2)  and  (3),  we  find  out  that  maximum  energy 
coupling  occurs  when  A(p(0)  =  ±7r/2.  For  these  two  limiting 
cases,  analytical  solutions  can  be  derived.  . By  examining  (3), 
we  can  easily  find  that  A<^  =  ±7r/2  are  solutions  to  (3).  The 
function  y  can  be  solved  by  substituting  (4)  back  to  the  first 
equation  of  (2)  and  is  given  by 

_ 4r^  exp(2f/) _ 

^  r^exp(4f/)-|-2(2+r2+r3)  exp(2f/)-hrQ ^(ra-ra)^ 

(6) 


where  /  is  the  same  as  in  (5),  r2  =  2P2(0)/Pi(0),  ra  = 
2Pa(0)/Pi(0)>  r  =  \/(l  +  r2)(l -f  ra),  and  tq  = 
\/Ta(r2  +  1)  +  v^T2(r3  -h  1).  When  r2  =  ra  =  r,  (6) 
reduces  to 

_  _ ^  ^ _  /-TN 

rexp[2(H-r)/]  +  l‘ 


In  this  case  (A<^  =  7r/2),  Pi  is  a  decreasing  function  of  2,  and 
P2,  Pa  are  increasing  functions  of  2,  according  to  (2).  In  other 
words,  the  energy  is  transferred  from  Pi  to  P2  and  Pa. 

When  A^(0)  =  '-7r/2,  analytical  solution  can  be  similarly 
obtained.  In  this  case,  the  energy  is  transferred  from  Po  and 
Pa  to  Pi.  The  energy  transfer  will  cease  if  either  P2  or  Pa 
is  depleted.  If  we  assume,  without  loss  of  generality,  that 
P2(0)  <  Pa(0),  then  the  solution  is  in  form  of 

Qi=2P2(0)u  +  Pi(0) 

Q2=-P2(0)u  +  P2(0), 

Q3=”/^2(0)u  +  P3(0)  (8) 


with  u  given  by  (9)  below,  where  ri2  =  Pi(0)/2P2(0).ra2  = 
Pa(0)/P2(0),/2  =  27P2(0)[1  -  exp(-Q2)]/a.  c  =  (1  + 
ri2)(Ta2  +  7*12)»  and  a  =  2yjcr^  +  +  ri2ra2  +  2r32.  It 

is  interesting  to  note  that  when  Pi(0)  =  0,  u  =  0,  which 
means  that  only  amplification,  instead  of  generation,  of  the 
center  component  is  possible. 

Equation  (9)  can  also  be  simplified  when  P2(0)  =  PaCO). 
In  this  case,  r32  =  1,  and  (9)  reduces  to 


1  -  exp(-2(l  +ri2)/2] 

. . . . -  '■■'■■■■  I  Ti  2 ♦ 

exp[-2(l  +  ri2)/2]  +  ri2 


(10) 


For  A0(O)  other  than  ±7r/2.  exact  analytical  solutions  in 
simple  form  are  not  available.  The  interference  between  the 
input  signals  and  the  newly  generated  FWM  components  may 
result  in  a  complicated  energy  coupling,  i.e.,  periodic  energy 
transfer.  However,  when  P3{0)  <  Pi(0),  we  expect  that 
the  energy  coupling  can  be  well  described  by  (4)  and  (5), 
as  pointed  out  earlier.  In  addition,  if  A0  does  not  change 
significantly  over  the  fiber  length,  an  approximate  analytical 
solution  can  be  easily  obtained.  The  presence  of  sin  A</)  in 
(2)  is  equivalent  to  a  simple  modification  of  the  nonlinear 
coefficient  7, 

The  analytical  solutions  can  be  used  to  investigate  energy 
coupling  and  depletion  in  fiber  communications  systems.  For  a 


_ 4cri2exp(~2v/c/2) _ 

aexp(-4v^/2)  -  2a(l  +  2ri2  +  r32)ri2  exp(-2yc/2)  +  r?2(r32  -  1)2/“ 


-  ri2. 


(9) 
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Fig.  2.  Output  powers  and  the  RVM-induced  losses  due  to  PDFW'M  in 
channel  fi  as  functions  of  input  power  for  multichannel  lightwave  systems. 
The  solid  lines  arc  for  systems  with  unequal  frequency  spacings.  and  dashed 
lines  are  for  systems  with  equal  frequency  spacings.  Equal  input  powers  for 
all  channels  are  assumed  for  both  cases. 


multichannel  lightwave  system  with  unequal  frequency  spac¬ 
ings,  energy  coupling  can  be  described  by  (4)  and  (5),  The 
coupling  will  lead  to  a  significant  energy  depletion  for  signal 
in  channel  /l,  thus  imposing  an  upper  limit  on  the  input 
power.  Further  increase  of  input  power  will  result  in  a  smaller 
output  power  in  channel  fi  due  to  the  FWM  process.  To 
see  this,  let  us  study  the  signal  evolution  after  transmission 
through  a  50  km  long  fiber.  For  simplicity,  we  assume  that 
all  channels  have  the  same  input  power.  We  plot  the  output 
power  Pi{L)  as  a  function  of  Pi(0)  in  Fig.  2  (solid  lines). 
We  note  that  for  input  power  greater  than  12  dBm,  the  output 
power  at  c  =  50  km  begins  to  drop.  Furthermore,  even  in 
regions  where  output  power  does  increase  with  the  increase 
of  input  power,  the  FWM-induced  loss  of  Pi  (defined  as 
l01og[Pf'^  is  considerable.  For  fibers 

with  larger  nonlinearity,  the  FW'M-induced  loss  will  be  even 
more  significant. 

For  multichannel  lightwave  systems  with  equal  frequency 
spacings,  we  can  use  ('6)-(l0)  to  describe  energy  coupling 
if  each  channel  has  a  finite  power  at  its  input  port.  For 
the  puipose  of  comparing  systems  with  unequal  frequency 
spacings,  we  limit  our  attention  to  the  case  of  A</>  =  7r/2, 
which  corresponds  to  cases  of  maximum  energy  coupling.  The 
dashed  lines  in  Fig.  2  shows  the  output  power  Pi(L)  and  the 
FWM-induced  loss  as  functions  of  input  power  Pi(0)  for  a 
fiber  length  of  L  =  50  km.  By  examining  the  figure  we  note 
that  the  signal  depletion  is  more  severe  compared  with  systems 
with  unequal  frequency  spacings  (solid  lines).  In  this  case, 
the  upper  limit  of  input  power  is  about  8  dBm.  Furthermore, 
the  FWM-induced  loss  is  much  larger  due  to  a  stronger 
signal  depletion  in  this  case.  This  verifies  again  the  conclusion 


that  multichannel  lightwave  systems  with  unequal  frequency 
spacings  are  preferable  to  those  with  equal  frequency  spacings 
[7].  To  further  reduce  the  impact  of  FWM  on  the  system 
performance,  all  input  powers  should  be  kept  below  a  certain 
value.  In  addition,  no  channel  frequency  should  be  the  same  as 
(or  even  too  close  to)  the  zero-dispersion  frequency  of  the  fiber 
so  that  phase  matching  condition  for  PDFWM  is  not  satisfied. 
For  a  given  channel  spacing  A  A  and  effective  fiber  length  £. 
the  minimum  separation  /i  -/o  should  be  such  that  A/cL  >  n. 
For  systems  operated  near  the  zero-dispersion  frequency,  Afc 
is  given  by  27r((iDc/dA)(AA)^(/i  —  /o),  where  Dc  is  the  fiber  ' 
chromatic  dispersion  [10].  /i  -  /o  can  be  roughly  estimated  as 
20  GHz  if  the  effective  fiber  length  is  20  km.  In  this  way,  the 
effect  of  FWM  can  be  minimized  due  to  phase  mismatching. 

In  conclusion,  we  have  studied  PDFWM  with  pump  de¬ 
pletion  from  the  physical  point  of  view.  Analytical  solutions 
for  PDFWM  with  all  inputs  present  have  been  derived  for 
multichannel  systems  operated  near  the  zero-dispersion  fre¬ 
quency.  It  has  been  shown  that  energy  can  be  coupled  from 
either  the  center  component  or  the  side  components  as  a  result 
of  interference  of  inputs  and  FWM-generated  components  in 
multichannel  systems  with  equal  frequency  spacing.  It  has  also 
been  shown  that  for  input  power  beyond  a  certain  value,  power 
in  the  pump  channel  will  be  considerably  depleted  ovef'a 
relative  short  distance.  To  ensure  good  system  performance, 
effective  suppression  of  FWM  is  critical. 
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The  performance  of  four  different  mutually  pumped  phase  conjugators  in  barium  titanate  depends  on  the  mu¬ 
tual  coherence  of  the  two  input  beams.  In  three  of  the  conjugators  the  use  of  fully  mutually  coherent  input 
beams  enhances  the  overall  phase-conjugate  reflectivity.  We  find  that  backscattering  and  transmission  grat¬ 
ings  in  the  photorefractive  crystal  both  contribute  to  the  phase-conjugate  signal  and  that  their  relative  strength 
depends  on  the  relative  coherence  of  the  two  input  beams.  We  numerically  solve  coupled- wave  equations  that 
include  all  gratings  and  find  reasonable  agreement  between  our  theory  and  our  experimental  data. 


1.  INTRODUCTION 

Imagine  trying  to  inject  light  into  the  gain  regions  of  a 
present-day  laser-diode  array.  Such  a  laser-diode  array 
might  measure  1  cm  across,  but  its  20  or  so  gain  regions 
are  each  only  200  /xm  across  X  1  ^m  high.  Injecting 
light  precisely  into  these  lasing  regions  is  like  trying  to 
thread  20  needles  in  a  haystack  simultaneously.  How¬ 
ever,  if  the  output  from  the  laser  array  is  sent  into  a 
photorefractive  crystal,  and  if  a  separate  light  beam  is 
also  directed  into  the  same  crystal,  then  the  crystal  could, 
in  principle,  connect  the  two  light  beams,  sending  each 
beam  back  down  the  path  of  the  other.  Acting  as  a  mu¬ 
tually  pumped  phase  conjugator,  the  crystal  would  trans¬ 
form  each  light  beam  into  the  phase-conjugate  replica  of 
the  other  light  beam.  In  the  above  application  the  crys¬ 
tal  would  direct  the  injected  light  beam  precisely  into 
the  lasing  regions  of  the  array.  The  two  incident  light 
beams  need  not  be  coherent  with  each  other,  although 
they  should  have  the  same  nominal  wavelength. 

Using  the  double  phase-conjugate  mirror  (DPCM),^  a 
number  of  authors  have  demonstrated  short-lived  injec¬ 
tion  locking  of  a  laser.^’^  MacCormack  et  al}  used  a  bird¬ 
wing®  conjugator  and  obtained  stable  locking  between  two 
laser  diodes. 

When  two  lasers  are  connected  by  a  phase  conjugator 
and  the  lasers  become  phase  locked,  their  light  beams 
become  mutually  coherent.  Can  the  crystal  continue  to 
act  as  a  mutually  pumped  phase  conjugator  with  mutu¬ 
ally  coherent  light  beams?  We  show  below  that  the  an¬ 
swer  depends  on  the  particular  geometry  of  the  mutually 
pumped  phase  conjugator. 

0740-3224/95/071363-07$06.00 


First  we  give  a  heuristic  explanation  of  how  mutu¬ 
ally  pumped  phase  conjugators  work.®*®  Let  the  two 
input  beams  be  labeled  2  and  4.  (We  use  this  coun¬ 
terintuitive  beam-numbering  scheme  to  stick  to  the 
convention  of  Ref.  7.)  Light  randomly  scattered  from 
beam  2  travels  across  beam  2  and  writes  a  multitude  of 
different  beam-fanning  gratings.  Beam  4  tends  to  erase 
most  of  these  gratings.  However,  beam  4  will  be  incident 
at  the  Bragg  angle  for  some  small  set  of  these  fanning 
gratings,  and  it  will  efficiently  deflect  from  them.  The 
deflected  beam  4  interferes  with  the  incident  beam  4  to 
reinforce  selectively  only  the  deflecting  gratings.  Simul¬ 
taneously  the  same  process  occurs  with  beams  2  and  4 
switched.  The  result  is  that  the  shared  fanning  gratings 
grow  in  strength,  whereas  the  unshared  fanning  gratings 
collapse.  The  shared  fanning  gratings  make  up  a  volume 
hologram  that  converts  beam  2  into  the  phase-conjugate 
replica  of  beam  4,  and  vice  versa.  The  above  description 
applies  to  the  DPCM.  The  bridge,®  the  bird-wing,  and 
the  British  (also  known  as  the  mutual  incoherent  beam 
coupler)®  conjugators,  all  shown  in  Fig.  1,  are  somewhat 
more  complicated;  they  have  more  than  one  interaction 
region,  and  in  the  latter  two  devices  the  light  beams  are 
also  internally  reflected  at  one  or  more  crystal  faces. 

Let  the  two  input  beams  have  wave  vectors  k2  and  ku, 
respectively.  If  these  two  beams  are  at  least  partially 
coherent  with  each  other,  then  they  will  create  gratings 
in  addition  to  the  shared  transmission  grating  described 
above.  These  additional  gratings  affect  the  performance 
of  the  conjugator.  Figure  2  shows  that  input  beam  2  can 
write  a  backscattering  grating,  with  wave  vector  k  =  2k2, 
with  the  counterpropagating  conjugate  beam  that  comes 
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Fig.  1.  Four  different  mutually  pumped  phase  conjugators. 
The  British  and  the  bird-Ns-ing  conjugators  use  internal  reflections 
to  guide  the  beams. 


conjugate  of  the  other  beam,  whereas  the  backscattering  gratings 
deflect  each  input  beam  into  its  own  phase  conjugate.  Backscat¬ 
tering  gratings  form  only  if  the  two  input  beams  are  at  least 
partially  mutually  coherent. 

from  beam  4.  Similarly,  along  the  path  of  beam  4  will 
be  a  backscattering  grating  with  wave  vector  k  =  2k4. 
We  note  that  if  the  input  beams  are  mutually  coherent 
they  will  also  write  a  reflection  grating  of  wave  vector 
k  =  ka  “  k4.  This  reflection  grating  contributes  only 
weakly  and  indirectly  to  the  phase-conjugate  signal,  and 
its  principal  effect  is  to  deflect  a  portion  of  the  input  beams 
out  of  the  interaction  region. 

In  general  there  are  two  primary  contributions  to  the 
phase-conjugate  signal:  (1)  the  beam  deflected  from  the 
shared  transmission  grating  and  (2)  the  beam  reflected 
from  backscattering  gratings.  The  strength  of  only  the 
backscattering  gratings  depends  on  the  degree  of  mu¬ 
tual  coherence  between  the  input  beams.  In  this  paper 
we  measure  how  these  two  contributions  to  the  phase- 
conjugate  signal  vary  as  we  alter  the  mutual  coherence 
of  the  two  incident  light  beams.  Our  measurement  tech¬ 
niques  are  similar  to  those  employed  by  Gnmeisen  et  aL, 
who  measured  the  total  phase-conjugate  signal  in  the 
DPCM.^°  We  also  present  theoretical  fits  to  our  data, 
using  a  recently  developed  theory  of  four-wave  mixing  of 
partially  coherent  beams  in  such  conjugators.^^ 


2.  EXPERIMENT 

We  operated  our  argon-ion  laser  at  A  —  514.5  nm  with 
an  etalon  to  increase  its  coherence  length  to  a  few  me¬ 


ters.  A  beam  splitter  (shown  in  Fig.  3)  generated  the 
two  extraordinary-polarized  input  beams.  We  matched 
the  optical  paths  of  the  two  beams  from  the  laser  to  the 
conjugator  to  within  10  cm.  We  equalized  the  beam  pow¬ 
ers  incident  at  the  entrance  to  the  conjugator  with  a  half¬ 
wave  plate  and  polarizer.  For  all  our  measurements  the 
input  beam  powers  were  typically  --14  mW/beam.  At 
the  conjugator  the  beams'  diameters  {lie  intensity)  were 
0.75  mm.  We  monitored  the  input  beam  power  and  the 
phase-conjugate  signal  for  both  input  beams,  using  reflec¬ 
tions  from  uncoated  optical  wedges.  We  positioned  an 
aperture  in  front  of  each  detector  to  ensure  that  we  mea¬ 
sured  only  the  phase-conjugate  beam.  A  Faraday  isola¬ 
tor  reduced  any  feedback  into  the  argon-ion  laser  by  at 
least  30  dB. 

We  used  nominally  undoped  photorefractive  crystals  of 
BaTiOa  to  form  the  different  conjugators.  Unable  to  find 
one  BaTiOa  crystal  that  could  support  all  four  conjugators, 
we  used  two  crystals  that  exhibited  good  beam  fanning: 
a  crystal  named  Hua  [6.05  mm  x  4.58  mm  x  5.89  mm 
(parallel  to  the  c  axis)],  for  the  bridge  and  the  bird¬ 
wing  conjugators,  and  a  crystal  named  Tvdn  [5.39  mm 
X  5.19  mm  x  5.40  mm  (c  axis)]  for  the  DPCM  and  the 
British  conjugators. 

We  controlled  the  mutual  coherence  of  the  two  input 
beams  as  follows.  We  placed  a  piezoelectric  transducer 
up  against  a  mirror  mount  to  give  beam  4  a  controllable- 
amplitude,  triangle- w'ave  phase  jitter.  We  set  the  jitter 
frequency  at  870  Hz  to  take  advantage  of  a  mechanical 
resonance  in  the  mirror  mount.  The  resulting  phase  jit¬ 
ter  caused  any  optical  interference  pattern  written  be¬ 
tween  the  two  input  beams  to  oscillate  spatially  at  870  Hz. 
This  frequency  was  an  order  of  magnitude  faster  than 
the  measured  response  rate  of  our  crystals  at  the  optical 
intensities  used  in  our  experiment,  so  the  crystal  effec¬ 
tively  time  averaged  the  phase  jitter  and  thereby  reduced 
the  effective  mutual  coherence  of  the  two  input  beams. 
Changing  the  amplitude  of  the  voltage  on  the  piezoelec¬ 
tric  transducer  varied  the  mutual  coherence  of  the  two 
input  beams.  We  observed  that  imparting  a  phase  jitter 
to  an  input  beam  had  no  effect  on  that  beam's  fanning  in 
the  crystal.  This  result  was  expected  because  in  beam 


Detector  for  phase  Detector  for  phase 

conjugate  of  beam  4  conjugate  of  beam  2 

Fig.  3.  Experimental  setup.  The  piezoelectric  pusher  imparts 
a  phase  jitter  to  cause  input  beams  2  and  4  to  be  partially 
mutually  incoherent,  as  measured  by  the  linear  photodiode  array. 
The  half-wave  plate  and  the  polarizer  form  a  variable  attenuator 
to  equalize  the  powers  of  the  two  input  beams  at  the  crystal. 


De  La  Cruz  et  al. 


Vol.  12.  No.  7/July  1995/J.  Opt.  Soc.  Am.  B  1365 


fanning  all  the  rays  originate  from  the  same  beam»  so 
the  same  phase  jitter  was  present  on  all  these  rays  and 
canceled  out. 

We  measured  the  effective  mutual  coherence  of  the 
two  input  beams  by  interfering  them  and  measuring 
their  fringe  visibility  t/  =  (/max  “  /minV(/max  +  /min)*  This 
fringe  visibility  ranges  from  =  0  with  mutually  incoher¬ 
ent  input  beams  to  v  =  I  with  fully  mutually  coherent 
input  beams.  We  monitored  the  interference  fringes  gen¬ 
erated  by  overlapping  equal  intensities  of  both  input 
beams  on  a  512-element  linear  photodiode  array  whose 
20-ms  readout  time  was  comparable  with  the  crystal’s  re¬ 
sponse  time.  We  determined  by  dividing  the  output 
from  the  linear  photodiode  array  by  the  output  produced 
with  mutually  incoherent  beams  =  0)  and  fitting  the 
resulting  curve  to  a  sine  wave  to  determine  its  visibility > 
as  shown  in  Fig.  4.  We  could  vary  the  fringe  visibility 
over  almost  its  full  range,  from  v  =  0.05  to  */  =  0,95. 

We  initially  set  up  four  different  mutually  pumped 
conjugacors,  using  mutually  incoherent  beams  at  angles 
similar  to  those  described  in  Refs.  1,  5,  8,  and  9.  We  then 
adjusted  the  beam  angles  and  positions  to  optimize  the 
stability  and  the  power  reflectivity  of  the  phase-conjugate 
signals  for  both  beams. 

Starting  from  fully  incoherent  beams,  we  decreased  the 
voltage  applied  to  the  piezoelectric  transducer,  measured 
the  corresponding  fringe  visibility,  and  then  allowed  the 
conjugator  to  stabilize  and  to  adapt  to  the  change  in 
mutual  coherence.  At  each  coherence  point  we  measured 
six  quantities: 

(i)  the  total  cw  phase-conjugate  reflectivity  of  beam  2; 

(ii)  the  total  cw  phase-conjugate  reflectivity  of  beam  4; 

(iii)  the  phase-conjugate  transmission  of  beam  4  (we 
blocked  beam  2  and  within  a  few  milliseconds  measured 
any  signal  that  was  phase  conjugate  to  beam  2;  this  was 
a  measure  of  the  fraction  of  beam  4  converted  into  the 
phase  conjugate  of  beam  2  because  the  shutter  momentar¬ 
ily  blocked  any  contribution  to  the  phase-conjugate  signal 
from  the  backscattering  of  beam  2); 

(iv)  the  phase-conjugate  transmission  of  beam  2,  de¬ 
tected  similarly  to  the  measurement  of  (iii)  above; 

(v)  the  purely  backscattered  contribution  to  the  phase 
conjugate  of  beam  2  (we  momentarily  blocked  beam  4  and 
measured  any  signal  that  was  phase  conjugate  to  beam  2); 

(vi)  the  purely  backscattered  contribution  to  the  phase 
conjugate  of  beam  4,  detected  similarly  to  the  measure¬ 
ment  of  (V)  above. 

We  usually  found  no  difference  if  we  approached  a 
particular  value  of  v  from  below  or  above;  however,  for  a 
small  number  of  cases  (not  shown  in  Figs.  5-9  below),  we 
observed  bistable  operation,  with  the  extreme  case  that 
some  conjugators  would  not  self-start  with  fully  coherent 
input  beams. 

When  using  fully  coherent  input  beams,  we  found  that 
tapping  one  of  the  beam-steering  mirrors  caused  the 
phase-conjugate  signal  either  to  drop  immediately  from 
a  maximum  value  or  to  increase  from  a  minimum  value, 
but  never  to  oscillate  about  an  intermediate  value.  This 
effect  implies  that  the  phase  shift  between  the  trans¬ 
mitted  and  the  backscattered  beams  was  always  0*  or 
180*,  as  expected  for  a  backscattering  grating  in  BaTiOa 


with  purely  diffusive  charge  transport  and  no  applied  dc 
electric  field. 

Figures  5-9  show  the  dependence  of  both  the  total 
phase-conjugate  reflectivity  and  the  transmission  on  the 
mutual  coherence  of  the  two  input  beams.  Both  quanti¬ 
ties  are  the  ratio  of  powers,  with  no  corrections  for  Fresnel 
loss.  The  reflectivity  data  include  both  the  transmitted 
and  the  backscattered  contributions  to  the  conjugate 
signals,  whereas  the  transmission  data  reveal  only  the 
transmitted  signal.  The  solid  curves  are  fits  to  the  total 
reflectivity,  as  detailed  in  Section  3.  The  phase- 
conjugate  reflectivities  and  transmissions  of  the  different 
conjugators  depended  on  the  particular  angles  and 
alignment  of  the  two  incident  beams.  Our  data  show  be¬ 
havior  typical  for  each  geometry,  given  that  we  optimized 
the  reflectivity  of  the  conjugators. 

With  mutually  coherent  input  beams  the  stability  of 
the  conjugate  signals  depended  critically  on  the  stability 
of  the  setup,  and  we  eventually  used  smaller  and  more 
stable  38.1-cm-high  pedestal  mounts  to  reduce  fluctua¬ 
tions  greatly.  Such  reduced  fluctuations  can  be  seen 
by  comparison  of  the  British  I  and  the  British  II  conju¬ 
gators  shown  in  Figs.  8  and  9,  respectively.  These  fluc¬ 
tuations  are  likely  caused  by  small  drifts  external  to 
the  crystal.  If  the  reflectivity  fluctuated  (usually  on  a 
time  scale  of  seconds),  we  plotted  its  maximum  value, 
and  we  represent  these  fluctuations  by  vertical  bars  on 
a  few  data  points.  The  horizontal  error  bars  represent 
the  uncertainty  in  the  relative  coherence  values.  This 
uncertainty  was  most  noticeable  at  the  higher  coherence 
values,  where  slight  air  currents  and  vibrations  became 


Fig.  4.  Measurement  of  the  mutual  coherence  v  of  the  two 
input  beams,  (a)  Interference  fringes  produced  by  two  partially 
mutually  coherent  beams  with  a  background  offset  subtracted, 
(b)  No  fringes  produced  by  two  completely  incoherent  beams, 
with  the  same  offset  subtracted.  Division  of  trace  (a)  by  trace 
(b)  produces  the  sinusoidal  fringe  pattern  (c),  which  in  this  case 
has  a  measured  visibility  —  0.54. 
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Mutual  coherence  (v) 


Fig.  5.  Behavior  of  the  DPCM  mutually  pumped  phase  con- 
jugator  as  we  alter  the  mutual  coherence  of  the  input  beams 
from  mutually  incoherent  (i^  =  0)  to  mutually  coherent  =  1). 
The  reflectivity  data  include  contributions  from  both  transmis¬ 
sion  and  backscattering  gratings,  whereas  the  transmission  data 
come  only  from  the  transmission  grating.  The  solid  curves  are 
our  theoretical  fits  to  the  two  total  reflectivities.  The  horizon¬ 
tal  and  the  vertical  bars  on  a  few  representative  data  points 
show  the  uncertainty  in  the  measured  coherence  values  and  the 
fluctuations  in  the  measured  signals.  Note  that  the  DPCM 
conjugator  turns  off  when  the  input  beams  become  too  coherent 
with  each  other. 


Mutual  coherence  (v) 


Fig.  6.  Bird- wing  conjugator.  The  symmetry  of  the  beam 
paths  inside  the  crystal  ensures  similar  data  for  the  two  input 
beams.  The  fits  are  to  the  reflectivity  data. 

important  and  limited  the  fringe  visibility  v  <  0.95.  All 
the  phase-conjugate  signals  were  stable,  with  mutually  in¬ 
coherent  input  beams.  Table  1  lists  our  subjective  com¬ 
parison  of  all  four  conjugators. 

The  British  conjugator  exhibited  different  behaviors  ac¬ 
cording  to  the  position  and  the  angle  of  the  input  beams. 
Figure  8  shows  data  for  British  I,  in  which  we  positioned 
both  input  beams  close  to  the  +c  face  of  the  crystal.  The 
resulting  coupling  channel  was  short,  and  the  curvature 
of  the  two  input  beams  had  a  pronounced  asymmetry  in¬ 
side  the  crystal,  with  beam  2  fanning  more  strongly  than 
beam  4.  In  contrast,  Fig.  9  shows  data  for  British  II,  in 
which  we  positioned  both  input  beams  slighter  further 
from  the  +c  face  and  at  smaller  angles  of  incidence.  This 
conjugator  developed  a  longer  coupling  channel,  with  both 
beams  now  having  similar  curvatxires  inside  the  crystal. 
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Mutual  coherence  (v) 

Fig.  7.  Bridge  conjugator.  The  asymmetr>-  shown  in  the  reflec¬ 
tivities  of  the  bridge  conjugator  is  caused  by  the  difference  in  the 
sign  of  the  backscattering  coupling  strength  for  the  two  input 
beams.  The  fits  are  to  the  reflectivity  data. 


Mutual  Coherence  (v) 

Fig.  8.  British  I  conjugator  has  a  small  loop  in  the  crystal.  The 
fits  are  to  the  reflectirity  data. 
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Mutual  coherence  (v) 

Fig.  9.  British  II  conjugator  has  a  larger  loop  in  the  crystal. 
The  large  signal  fluctuations  shown  in  these  data  were  caused  by 
a  less  stable  optical  setup.  When  we  switched  to  a  more  stable 
setup  the  British  conjugator  was  the  least  noisy  of  the  four  types 
of  conjugator  investigated.  The  fits  are  to  the  reflectivity  data. 
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Table  1.  Ease-of-Use  and  Signal  Fluctuations 
of  the  Mutually  Pumped  Phase  Conjugators 


Conjugator 

Ease  of  Setup 

Signal  Fluctuations 

DPCM 

Easy 

Moderate 

Bird  wing 

Somewhat  easy 

Small 

Bridge 

Somewhat  difficult 

Small 

British 

Difficult 

Very  small 

In  general,  as  we  increased  the  mutual  coherence  of  the 
two  input  beams  at  least  one  conjugate  signal  increased, 
whereas  both  transmissions  decreased.  The  exception 
was  the  DPCM  shown  in  Fig.  5,  which  usually  ceased  to 
operate  when  the  input  beams  became  coherent  with  each 
other  (1/  =  1).  We  could  obtain  phase-conjugate  signals 
from  the  DPCM  with  fully  coherent  input  beams,  but  only 
after  tedious  alignment.  Both  phase-conjugate  outputs 
in  the  bird-wing  conjugator  and  the  British  II  conjuga- 
tor  increased  as  ^  »  1,  whereas  in  the  bridge  and  the 

British  I  conjugators  one  conjugate  signal  grew,  but  the 
other  diminished  in  power.  At  1,  we  checked  that 
the  bridge,  the  bird-wing,  and  the  British  conjugators 
were  not  acting  as  self-pumped  conjugators  by  confirm¬ 
ing  that  both  input  beams  were  required  for  sustained 
conjugation. 

In  the  DPCM  we  occasionally  observed  that  the  phase- 
conjugate  reflectivity  of  beam  2  increased  steadily  as  v  in¬ 
creased  from  0  to  0.25,  where  the  reflectivity  then  dropped 
abruptly,  as  predicted  by  our  theoretical  plots.  However, 
most  of  the  time  the  reflectivity  of  beam  2  in  the  DPCM 
began  to  collapse  at  *=*  0.1,  as  shown  in  Fig.  5.  The 
DPCM  data  presented  here  agree  in  general  with  those 
reported  by  Gruneisen  et  a/.,^°  with  the  exception  that 
their  DPCM  conjugated  weakly  even  with  high  mutual 
coherence.  We  note  that  our  theory  does  not  take  into 
account  the  effects  of  conical  scattering,  which  can  cause 
significant  loss.^"*  Also,  unequal  fanning  losses  experi¬ 
enced  by  beams  2  and  4  can  cause  the  different  reflectivi¬ 
ties  observed  for  the  two  beams  at  i/  =  0. 

As  the  mutual  coherence  of  the  two  input  beams  is  in¬ 
creased,  backscattering  gratings  increase  in  strength  and 
affect  the  performance  of  the  conjugators  in  three  ways. 
First,  each  input  beam  can  now  exchange  energy  with  its 
own  phase-conjugate  beam  through  two-beam  coupling, 
with  the  direction  of  energy  transfer  being  determined 
by  the  orientation  of  the  beams  with  respect  to  the  crys¬ 
tal’s  c  axis.  Second,  this  two-beam  coupling  can  create 
an  imbalance  in  the  two  input  beams’  intensities  inside 
the  crystal;  this  imbalance  reduces  the  cross-erasure  effi¬ 
ciency,  which  weakens  the  transmission  grating.  Third, 
backscattering  gratings  can  decrease  the  overall  trans¬ 
mission  of  a  conjugator  by  simply  reflecting  energy  away 
from  the  transmitted  phase-conjugate  beam. 

Figure  10  shows  how  the  computed  amplitude  cou¬ 
pling  coefficient  of  a  backscattering  grating  varies  as 
a  function  of  the  angle  8  between  the  beams*  propagation 
directions  inside  the  crystal  and  the  c  axis.*^®  For  this 
calculation  we  use  recent  values  for  the  effective  Pockels 
coefficients,  which  include  the  small  influence  of  elastic 
stiffness,  piezoelectric  stress,  and  the  elasto-optic  effect.^® 
We  find  that  input  beams  traveling  at  angles  6“  <  ^  < 
90*  are  depleted  by  backscattering  gratings.  We  esti¬ 
mate  that,  in  the  bird-wing  and  the  British  conjugators, 
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the  angles  of  both  conjugate  beams  with  respect  to  the 
c  axis  lie  in  the  range  95*  <  8  <  150*,  so,  according  to 
Fig.  10,  in  these  conjugators  both  conjugate  beams  gain 
from  backscattering  gratings  at  the  expense  of  the  input 
beams.  To  estimate  the  contribution  of  backscattering 
gratings  to  the  phase-conjugate  signal  we  use  the  relation 

backscatter  reflectivity 

=  (^totai  reflectivity  7:  Vtransmission  , 

where  the  terms  are  added  if  the  components  are  out 
of  phase  or  are  subtracted  if  they  are  in  phase,  and  we 
neglect  the  contribution  of  reflection  grating  to  the  con¬ 
jugate  signal.  For  example,  for  the  bridge  conjugator 
the  computed  backscattering  power  reflectivity  for  beam  2 
with  fully  mutually  coherent  beams  (u  =  1)  is  (>/33^ 
-  V2.4%)^  =  17.6%.  This  result  agrees  with  our  mea¬ 
sured  backscattered  power  reflectivity  of  18%.  We  per¬ 
formed  such  calculations  for  all  the  data  shown  and  found 
good  agreement  between  the  calculated  and  the  mea¬ 
sured  backscattered  contributions.  This  finding  confirms 
that  the  conjugate  signals  are  due  to  transmission  and 
backscattering  gratings  and  that  reflection  gratings  can 
be  neglected  in  the  analysis  of  these  experiments. 

Backscattering  gratings  affect  the  different  conjugators 
differently.  The  effect  of  backscattering  gratings  is  seen 
most  clearly  in  the  data  of  Fig.  7  for  the  bridge  conjugator. 
Here  one  conjugate  beam  travels  somewhat  opposite  the  c 
axis  at  0  --  130*  and  so  gains  from  its  backscattering  grat¬ 
ing,  whereas  the  other  conjugate  beam  is  at  0  -  70*  and 
so  is  depleted  by  its  backscattering  grating.  This  effect 
causes  the  asymmetry  seen  in  the  reflectivity  data.  In 
contrast,  in  the  bird-wing  conjugator  the  path  of  the  two 
beams  are  symmetric  about  the  c  axis,  and  both  conju¬ 
gate  beams  gain  equally  by  backscattering.  Figure  6 
confirms  that  making  the  input  beams  coherent  has 
little  effect  on  the  transmission  grating,  inasmuch  as 
the  strength  of  this  grating  depends  only  on  the  rela- 


Fig.  10.  Amplitude  two-beam  coupling  efficient  yt  for  a  beam 
and  ita  counterpropagating  beam  as  a  function  of  the  internal 
angle  between  the  beam  and  the  BaTiOa  crystal’s  c  axis.  A 
positive  coupling  coefficient  means  that  the  beam  is  gaining 
energy  from  its  counterpropagating  beam.  For  example,  a  beam 
propagating  between  6.7*  and  90*  to  the  +c  axis  is  depleted 
by  its  counterpropagating  beam.  This  calculation  uses  recent 
values  for  the  electro-optic  coefficients^^  and  assumes  an  effective 
charge-carrier  density  of  Af  ■*  5  x  10^®  cm"®  and  e-polarized  rays 
at  A  -  514.5  nm. 
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tive  intensity  of  the  two  input  beams  and  not  on  their 
mutual  coherence.  For  the  DPCM,  as  the  input  beams 
are  made  more  coherent  with  each  other,  the  resulting 
backscattering  gratings  amplify  input  beam  4  and  deplete 
input  beam  2,  and  the  resulting  intensity  imbalance  so 
weakens  the  primary  transmission  grating  that  it  col¬ 
lapses,  as  shown  in  Fig.  5.  We  confirmed  this  effect  by 
increasing  the  power  of  input  beam  2  in  the  DPCM  while 
keeping  the  power  of  beam  4  constant.  We  found  that 
conjugation  would  then  persist  even  with  fully  coherent 
beams  because  now  the  intensities  of  the  input  beams 
were  more  closely  matched  inside  the  crystal.  In  con¬ 
trast,  conjugation  ceased  when  we  decreased  the  power 
of  beam  2.  Evidently  the  bridge  and  the  British  con- 
jugators  fall  somewhere  between  the  bird  wing  and  the 
DPCM.  These  conjugators  experience  an  intensity  im¬ 
balance  with  mutually  coherent  input  beams,  but  not 
enough  to  cause  the  transmission  grating  to  collapse. 


3.  THEORY 

We  write  coupled-wave  equations  that  include  transmis¬ 
sion,  reflection,  and  backscattering  gratings,^  and  we 
also  include  the  mutual  coherence  parameter  between 
beams  i  and  For  a  single  interaction  region  we  obtain 

dAi  AiA,*  +  A2*A3  , 

17-" - h - 
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The  bovindary  conditions  are 


Ai(0)  =  A3(L)  =  0, 

|A4(0)I" 

|A2(L)12  * 

where  A2  and  A*  are  amplitudes  of  the  input  beams  and 


Al  and  A3  are  their  respective  conjugates;  Iq  =  |AiP 
IA2P  lAsl^  +  |A4l^  is  the  total  intensity;  and  7,,  7,., 
yf),  76'  are  the  coupling  coefficients  for  the  transmission, 
the  reflection,  and  the  two  backscattering  gratings,  re¬ 
spectively.  The  backscattering  gratings  for  the  two  in¬ 
put  beams  can  have  different  strengths  in  general.  L  is 
the  interaction  length,  and  z  is  the  direction  of  propaga¬ 
tion  for  beams  1  and  4.  Also,  a  is  the  bulk  absorption 
coefficient,  and  q  is  the  power  ratio  of  the  two  incident  in¬ 
put  pump  beams  just  inside  their  respective  crystal  faces. 
When  the  two  input  beams  are  mutually  incoherent,  i.e., 
when  1^24!  =  I/Sisl  =  1^12!  =  I/334I  =  0,  the  coupled-wave 
equations  reduce  to  the  conventional  equations,  where 
only  the  transmission  grating  is  present.®  Because  the 
photorefractive  effect  in  BaTiOa  crystals  is  dominated  by 
diffusive  charge  transport,  here  the  coupling  coefficients 
are  all  real  numbers.  We  simplify  the  coupled  equations 
by  approximating  1^24!  =  i^ial  =  lySnl  =  t^34l  =  \P\  be¬ 
cause  we  calculate  that  under  experimental  conditions  the 
differences  between  the  |^,^|  values  are  less  than  —5%. 
Although  our  theory  is  based  on  a  one-interaction-region 
model,  the  equations  can  be  extended  to  multiple  regions 
with  different  coupling  strengths.^® 

We  solve  the  coupled- wave  equations  with  standard 
fourth-order  Runge-Kutta  and  shooting  algorithms. 
The  mutual  coherence  is  equal  to  the  measured  fringe 
visibility,  |/3l  =  Py  because  the  two  beams  that  we  inter¬ 
fere  in  our  visibility  measurements  had  equal  intensity.^" 
We  use  our  measured  values  for  the  bulk  absorption  of  the 
BaTiOa  crystal  (a  -  2,6  cm"^  in  Hua  and  a  =  1.5  cm”^ 
in  Twin).  For  each  conjugator  we  calculate  the  ratio  q 
of  the  input  beam  powers  just  inside  the  entrance  face  of 
the  crystal  from  our  measured  incident  beam  power  and 
the  measured  Fresnel  loss  at  the  crystal  surface. 

We  fit  our  experimental  data  with  four  coupling- 
strength  parameters  {yt  •  L,  7r  •  L,  75  •  L,  76'  *  L).  We 
first  choose  7(  •  L  to  make  the  computed  reflectivities 
match  the  experimental  values  measured  with  mutually 
incoherent  beams  (v  =  0).  We  then  adjust  the  values 
of  the  remaining  parameters  to  best  fit  the  measured 
phase-conjugate  reflectivities  with  ^  0,  as  shown  by 
the  solid  curves  in  Figs.  5~9.  For  DPCM,  where  the 
input  beams  are  nearly  counterpropagating,  we  impose 
the  additional  restriction  75  =  -76'-  Table  2  shows  that 
the  resulting  gain -length  products  are  well  within  the 
range  expected  for  BaTiOa  and  have  the  correct  sign  for 
the  given  beam  directions  in  the  crystal.  For  example, 
from  the  measured  angles  of  incidence  of  the  input  beams 
and  the  calculated  backscattering  coefficient  of  Fig.  10 


Table  2.  Values  of  the  Fitting  Parameters® 


Conjugator 

yt-L 

yr-L 

yb-L 

76'  *  L 

a  ’  L 

L  (mm) 

DPCM 

2.8 

2.2 

1.1 

-1.1 

0.5 

6.0 

Bridge 

2.9 

1.0 

4.1 

-0.4 

0.65 

4.5 

Bird  wing 

3.3 

0.5 

2.3 

2.4 

0.9 

8.0 

British  I 

3.9 

0.7 

2.4 

1.2 

1.1 

11.0 

British  II 

3.3 

1.1 

2.1 

0.6 

0.76 

11.5 

“Amplitude  gain  times  length  for  the  transmission,  the  reflection,  the 
two  backscattering  gratings,  and  the  intensity  absorption  used  in  our 
theoretical  fits.  A  negative  value  for  the  backscattering  coefficient  corre¬ 
sponds  to  a  depletion  of  the  conjugate  signal.  Although  not  required  for 
the  fits,  our  estimate  of  the  interaction  length  (2:1  mm)  inside  each  crystal, 
obtained  by  observation  of  the  actual  beam  paths,  is  also  provided. 
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we  expect  that,  in  the  DPCM  and  in  the  bridge  conjuga- 
tor,  backscattering  gratings  will  decrease  the  conjugate 
of  beam  4,  and  not  of  beam  2.  In  these  two  conjugators 
the  phase-conjugate  to  beam  4  travels  at  angles  6"  <  ^  < 
90’  with  respect  to  the  positive  c  axis,  which  makes 
7ft/  <  0.  For  all  other  cases  7i  >  0  smd  7ft/  >  0,  and  the 
backscattering  gratings  couple  energy  from  each  input 
beam  into  its  respective  phase-conjugate  beam.  We  no¬ 
ticed  that  increasing.the  value  of  yr  always  decreases  the 
strength  of  the  conjugate  signal  for  all  four  conjugators. 

4.  CONCLUSION 

With  two  mutually  incoherent  light  beams  incident  upon  a 
mutually  pumped  phase  conjugator,  each  phase-conjugate 
signal  comes  exclusively  from  the  other  light  source. 
However,  when  the  two  incident  light  beams  are  mutu¬ 
ally  coherent,  each  phase-conjugate  signal  is  made  up 
of  a  mixture  of  the  light  from  both  light  sources.  From 
the  viewpoint  of  locking  lasers  with  such  conjugators, 
both  the  transmitted  and  the  backscattered  components 
are  useful  for  injection  locking  because  both  are  at  the 
same  frequency  and  are  phase  conjugate  to  the  output 
of  the  slave  laser.  However,  because  the  backscattering 
grating  cannot  respond  to  any  light  modulations  faster 
than  the  grating’s  response  time,  only  the  transmitted 
component  can  phase  lock  or  phase  modulate  a  laser. 

Of  the  four  mutually  pumped  phase  conjugators  investi¬ 
gated  here,  only  the  bird-wing  and  the  British  conjugators 
provide  two  conjugate  output  beams  when  the  two  input 
beams  are  fully  mutually  coherent.  The  bridge  conjuga¬ 
tor  produces  only  one  output  beam  of  significant  inten¬ 
sity  when  the  beams  are  coherent.  This  result  might  be 
advantageous  in  laser-locking  schemes  in  which  feedback 
into  the  injecting  laser  is  undesirable.  The  DPCM  ceases 
operation  with  fully  mutually  coherent  beams.  A  one- 
region,  one-dimensional  numerical  model  was  able  to  fit 
the  main  trends  of  our  data  with  reasonable  values  for  all 
the  fitting  parameters. 
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We  report  the  results  of  experimental  investigations  of  the  mechanism  transformation  from  four- wave  mixing 
and  stimulated  photorefractive  backscattering  to  four-wave  mixing  and  total  internal  reflection  during  self- 
pumped  phase  conjugation  in  a  photorefractive  crystal  by  varying  the  plane  of  polarization  of  the  input  beam. 

The  use  of  a  mixed-polarization  input  beam  can  change  the  fanning  pattern,  which  leads  to  the  mechanism 
transformation  and  an  enhancement  of  the  phase-conjugation  reflectivity.  The  temporal  dynamics  of  the 
buildup  process  indicate  that  such  an  enhancement  is  due  to  the  presence  of  both  types  of  mechanism  during 
phase  conjugation.  In  addition,  we  compare  the  contributions  of  2k  gratings  with  the  phase  conjugation 
during  both  mechanisms.  ©  1995  Optical  Society  of  America 


The  phenomenon  of  self-pumped  phase  conjugation 
(SPPC)  has  received  a  great  deal  of  attention  be¬ 
cause  of  its  important  potential  applications.  In  1982, 
Feinberg^  demonstrated  a  SPPC  mirror  by  using  a 
single  crystal  without  any  external  mirrors.  Since 
then,  there  has  been  much  interest  in  trying  to  under¬ 
stand  the  physical  mechanisms.  In  1983,  MacDonald 
and  Feinberg^  proposed  a  model  in  which  the  phase- 
conjugate  wave  was  generated  by  means  of  four-wave 
mixing  (FWM)  in  two  separate  interaction  regions. 
These  two  interaction  regions  were  coupled  by  total  in¬ 
ternal  reflections  within  the  crystal.  This  is  known 
as  the  four-wave  mixing  and  total  internal  reflection 
(FWM-TIR),  or  Cat,  model.  Recently  we^"^  proposed 
an  alternative  model  [four-wave  mixing  and  stimulated 
photorefractive  backscattering  (FWM-SPB),  in  which 
the  SPB  is  produced  by  the  formation  of  2k  grat¬ 
ings]  to  explain  the  origin  of  the  SPPC  observed  in 
KTal^Os  and  BaTi03:Ce  crystals.  In  both  FWM-TIR 
and  FWM-SPB  conjugators  the  SPPC  is  established 
through  FWM  interactions  when  the  fanning  beam  is 
retroref lected  by  TIR  or  SPB.  In  the  FWM-TIR  phase 
conjugator,  the  retroreflection  of  the  fanned  beam  is 
due  to  two  successive  total  internal  reflections  from 
one  of  the  crystal's  comers.  Thus  a  feedback  opti¬ 
cal  loop  with  two  FWM  interaction  regions  inside  the 
crystal  can  usually  be  seen.  In  the  FWM-SPB  phase 
conjugator  the  retroreflection  of  the  fanned  beam  is 
generated  by  SPB,  which  may  be  seeded  by  random 
scattering  centers  in  the  crystal  or  on  the  crystal's 
s\irfaces  or  even  by  an  external  retroreflector.  As  a 
result,  individual  filaments  that  contain  coimterpropa- 
gating  beams  can  often  be  seen. 

In  a  typical  SPPC  experimental  configuration  an  ex¬ 
traordinarily  polarized  input  beam  is  incident  upon  the 
a  face  of  a  regular-cut  BaTiOa  photorefractive  crystal. 
The  SPPC  formation  mechanism  is  usually  deter¬ 
mined  by  the  fanning  pattern  and  the  boundary  condi¬ 
tions.  Generally,  SPPC  is  generated  by  the  F\^-SPB 
mechanism  in  some  doped  crystals  but  by  the  FWM- 
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TIR  mechanism  in  undoped  crystals.  This  is  due  to 
stronger  fanning  and  larger  coupling  coefficients  for  2k 
gratings  in  doped  crystals.  Recently  we  observed  that 
the  formation  mechanism  of  SPPC  was  transformed 
from  FWM-SPB  to  FWM-TIR  when  the  wavelength  of 
the  input  beam  was  increased.^  In  this  Letter  we  re¬ 
port,  for  the  first  time  to  our  knowledge,  the  mecha¬ 
nism  transformation  during  SPPC  in  a  given  crystal 
and  at  a  fixed  wavelength  that  is  due  to  the  variation 
of  the  polarization  state  of  the  input  beam.  We  also 
show  that  the  phase-conjugation  (PC)  reflectivity  can 
be  increased  by  changing  the  plane  of  polarization  of 
the  input  beam.  In  this  regime  of  high  reflectivity, 
both  mechanisms  are  responsible  for  the  PC.  In  addi¬ 
tion,  we  compare  the  contributions  of  2k  gratings  to  PC 
during  both  mechanisms  by  using  an  ordinarily  polar¬ 
ized  erasure  beam. 

In  our  experimental  demonstration  we  used  a  Ce- 
doped  BaTiOs  crystal  that  had  a  red-orange  color  and 
a  dopant  concentration  of  approximately  30  parts  in 
10®.  The  crystal  was  poled  and  0*  cut  with  dimen¬ 
sions  6.15  mm  X  5.20  mm  x  8.2  mm  (a  x  6  x  c).  In 
all  our  SPPC  experiments,  we  utilized  the  output  from 
an  Ar"*"  laser  at  514.5  nm,  extraordinarily  polarized, 
directly  incident  upon  the  a  face  of  our  crystal  at  an 
external  angle  of  0  =  60**.  The  experimental  setup  is 
shown  in  Fig.  1.  With  an  input  intensity  of  approxi¬ 
mately  1  W/cm^,  the  plane  of  polarization  of  the  input 
beam  could  be  rotated  by  a  half-wave  plate. 

In  BaTi03:Ce  cryst^s,  extraordinarily  polarized 
beams  produce  the  strongest  fanning  effects.  In 
our  experiment  we  were  able  to  change  the  fanning 
pattern  by  varying  the  plane  of  polarization  of  the 
input  beam.  Figure  2  shows  the  images  of  these 
different  fanning  patterns  at  various  polarization 
states  before  the  phase  conjugation  reaches  its  steady 
state.  It  is  known  that  the  coupling  coefficient  in 
BaTi03:Ce  crystals  for  an  ordinarily  polarized  beam 
is  much  less  than  that  of  an  extraordinarily  polarized 
beam,  and  the  cross  coupling  is  negligible.®  Thus  the 
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polarization-dependent  mechanism  transformation  during 
SPPC.  The  external  incident  angle  was  ^  =  60*,  and  the 
incident  position  d  was  4  mm.  As  shown  in  the  inset, 
the  polarization  state  ^  of  the  input  beam  is  defined  as  the 
angle  of  the  polarization  p  measured  from  yz  plane.  PD, 
photodetector;  BS,  beam  splitter. 


Fig.  2.  Images  of  fanning  patterns  in  o\ir  BaTiOa'Ce 
crystal  at  different  polarization  angles  of  the  incident  beam 
before  SPPC  reaches  its  steady  state:  (a)  <I>  =  0%  (b) 
4)  =  20%  (c)  =  60%  The  c  axis  of  the  crystal  is  directed 

from  the  top  to  the  bottom  of  the  image,  and  the  input  beam 
is  propogating  from  left  to  right. 


extraordinarily  polarized  component  will  contribute 
significantly  to  the  fanning,  whereas  the  ordinarily 
polarized  component,  not  being  phase  matched  to  the 
gratings  written  by  the  extraordinary  component,  will 
only  serve  to  erase  some  of  the  index  gratings  and  to 
reduce  the  effective  coupling  constant.  From  Fig.  2, 
we  also  note  that  the  fanning  patterns  move  to  the 
comer  of  the  crystal  as  <I>  increases.  This  means  that 
the  fanning  (or  the  fanning  coupling  constant)  can  be 
effectively  controlled  by  use  of  a  mixed-polarization 
beam.  We  believe  that  these  changes  of  the  fanning 
pattern  and  the  effective  coupling  constant  are  respon¬ 
sible  for  the  mechanism  transformation  during  SPPC 
in  BaTiOa^Ce  crystals. 

Figure  3(a)  shows  the  beam  path  in  the  crystal  at 
d>  =  0®  when  the  PC  signal  reached  steady  state, 
demonstrating  a  typical  FWM-SPB  mechanism.  The 
input  beam  was  depleted  completely  after  propagating 
through  a  short  distance  and  bent  toward  the  bottom 
surface  of  crystal  because  of  strong  fanning.  The 
beam  pattern  at  steady  state  is  characterized  by  many 
filaments  that  are  produced  through  the  SPB  process 
by  2k  gratings.  This  SPPC  formation  is  significantly 
different  from  the  traditional  FWM-TIR  mechanism. 
There  is  no  closed  optical  path  involving  the  comer 
of  the  crystal.  The  mechanism  of  this  kind  of  phase 
conjugation  was  previously  modeled  and  explained 
in  Ref.  3. 


Figure  3(b)  shows  a  typical  beam  path  in  the  crystal 
at  4>  =  20%  when  two  types  of  mechanism  exist  si¬ 
multaneously.  The  time  response  of  the  self-pumped 
phase-conjugate  reflectivity  was  found  to  consist  of 
a  two-plateau  structure  (Fig.  4).  When  the  input 
beam  was  turned  on,  an  initial  quasi-steady-state 
SPPC  reflectivity  was  reached,  and  two  filaments 
appeared.  This  looked  like  a  FWM-SPB  phase  con- 
jugator.  After  this  stage  a  weak  loop  started  to 
appear  and  the  SPPC  reflectivity  gradually  reached 


Fig.  3.  Images  of  steady-state  beam  paths  during  SPPC 
at  three  different .  polarization  angles  of  the  incident 
beam: .  (a)  <I>  =  0%  (b)  <D  =  20%  (c)  <I>  -  60% 


Fig.  4.  Time  evolution  of  the  phase-conjugate  signal. 
The  pump  intensity  was  approximately  1  W/cm^. 


1.4 
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Polarization  angle  of  input  beam  (<I>) 

Fig.  5.  Normalized  PC  reflectivity  as  a  fimction  of  the 
polarization  angle  <I>  of  the  input  beam.  These  data  are 
not  corrected  for  Fresnel  reflections. 


August  15,  1995  /  VoL  20,  No.  16  /  OPTICS  LETTERS  1685 


Intensity  of  erasing  beam  (W/cm^) 

Fig.  6.  PC  reflectivity  as  a  function  the  intensity  of  the  or¬ 
dinarily  polarized  erase  beam.  Note  the  phase-conjugate 
threshold  at  an  erase  intensity  of  =4  W/cm^.  The  solid 
curve  is  from  a  theoretical  calciilation  for  the  FWM-SPB 
case. 


a  second  steady  state.  This  means  that  the  FWM- 
SPB  mechanism  was  mainly  responsible  for  the 
SPPC  formation  at  first,  and  then  the  FWM-TIR 
mechanism  gradually  settled  in.  When  the  ordi¬ 
narily  polarized  component  of  the  input  beam  was 
increased  by  rotation  of  the  input  plane  of  polarization, 
the  filaments  gradually  moved  to  the  comer  of  the  crys¬ 
tal.  The  optical  beam  paths  supporting  SPPC  were 
gradually  transformed  into  a  FWM-TIR  arrangement. 
Figure  3(c)  shows  the  case  with  —  60%  in  which  the 
optical  path  loop  at  the  bottom-right-hand  comer  is 
clearly  shown. 

The  PC  reflectivity  as  a  function  of  the  plane 
of  polarization  O  of  the  input  beam  is  shown  in 
Fig.  5.  Here  the  reflectivity  is  defined  as 
where  is  the  PC  output  and  h  is  the  intensity 
of  the  extraordinary  component  of  the  input  beam. 
We  note  that  the  phase-conjugate  output  is  always 
extraordinarily  polarized  regardless  of  the  input 
polarization  state.  The  reflectivity  is  normalized 
to  that  of  a  totally  extraordinarily  polarized  state  of 
input  beam  (O  =  0%.  We  observed  that  there  was  a 
reflectivity  enhancement  of  ^20%  for  <I>  =  20"  when 
the  two  SPPC  mechanisms  existed  simultaneously. 
This  leads  to  the  conclusion  that  optimal  SPPC  output 
is  produced  by  utilization  of  both  types  of  mecha¬ 
nism.  The  phenomena  of  mechanism  transformation 
during  SPPC  and  the  enhancement  of  phase-conjugate 
reflectivity  are  due  to  the  presence  of  an  ordinarily 
polarized  component  in  the  input  beam  that  reduces 
the  strength  of  the  fanning  pattern.  Thus  the  energy 
loss  that  is  due  to  fanning  can  be  limited. 

To  aiustrate  further  the  difference  between  the 
FWM-SPB  and  FWM-TIR  mechanisms,  we  uniformly 
illuminated  the  filaments  by  using  a  second  erasure 
beam  incident  from  the  right-hand  side  of  the  crysl^ 
(see  Fig.  1).  The  erasure  beam  was  incoherent  with 
the  input  beam  and  was  ordinarily  polarized.  The 
erasure  beam  reduced  the  modulation  depth  of  the 
index  gyrating  and  the  coupling  coefficients  throughout 


the  region  of  the  filaments  (or  the  loop).  We  measured 
the  PC  reflectivity  as  a  function  of  the  intensity  of  the 
erasure  beam,  /o,  to  compare  the  effect  of  2k  grat¬ 
ings  in  both  FWM-SPB  and  FWM-TIR.  The  back- 
scattering  coupling  coefficient  y2k  ^lay  be  rewritten  as 
y2k{Io)  =  r2*/(l  ^  h/J)*  where  Iq  is  the  intensity  of  the 
erasure  beam,  I  is  the  total  intensity  of  the  forward- 
propagating  beam  (fanning  beam)  and  the  backward- 
propagating  beam  in  the  filament  (or  the  loop), 
and  y2k  is  the  coupling  coefficient  at  /q  =  0.  Figure  6 
shows  the  phase-conjugate  reflectivity  as  a  function 
of  /o  with  an  interaction  length  of  L  =  3  mm  for 
the  2k  gratings.  We  note  that  the  PC  reflectivity 
decreased  as  the  intensity  of  erasure  beam  increased 
and  exhibited  a  sharp  threshold  at  /q  ^  4  W/cm^  in 
the  FWM-SPB  model  [shown  in  Fig.  3(a)].  In  our 
experiment,  I  was  approximately  2  W/cm^,  and  the 
average  coupling  coefficient  ^he  case  shown  in 

Fig.  3(a)  was  approximately  10  cm"%  In  the  FWM- 
TIR  model  of  SPPC,  the  decrease  in  reflectivity  (which 
is  the  ratio  of  the  PC  signal  with  and  without  the 
erase  beam)  is  approximately  10%  when  the  intensity 
of  the  erasure  beam  is  7o  4  W/cm^.  We  conclude 
therefore  that  the  2k  gratings  in  the  loop  make  very 
little  contribution  to  the  phase  conjugation.  This 
result  is  in  agreement  with  the  research  of  Ref.  6.  In 
FWM-SPB,  the  initial  retroreflections  from  crystal 
surface  imperfections  are  very  weak.  Thus  the  pres¬ 
ence  of  2k  gratings  is  essential  for  enhancement  of  the 
retroreflection.  In  the  case  of  FWM-TIR,  the  comer 
cube  is  already  providing  nearly  a  100%  retroreflection 
regardless  of  the  presence  of  gratings. 

In  summary,  we  have  demonstrated  a  polarization- 
dependent  mechanism  transformation  during  SPPC 
in  a  photorefractive  BaTiOa^Ce  crystal.  We  have  also 
shown  that  the  introduction  of  an  ordinarily  polarized 
component  in  the  input  beam  can  modify  the  fanning 
pattern  and  enhance  the  PC  reflectivity.  In  addition, 
we  have  compared  the  contribution  of  the  2k  gratings 
with  the  phase  conjugation  for  both  mechanisms. 
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Abstract.  We  investigate  the  mutually  pumped  phase- 
conjugation  process  in  photorefractive  media  with  par¬ 
tially  coherent  pump  beams.  We  consider  the  effects  of 
transmission,  reflection,  and  2A:-reflection  gratings  as 
well  as  the  depletion  of  both  pump  beams.  The  theoreti¬ 
cal  results  can  be  employed  to  qualitatively  explain  the 
experimental  observations. 

PACS:  42.65.Hw;  42.80;  42.70 


The  Mutually  Pumped  Phase-conjugation  (MPPC) 
process  is  one  of  the  most  important  discoveries  in  the 
field  of  photorefractive  nonlinear  optics  over  the  past 
fifteen  years  [1-3].  Since  it  does  not  require  mutual 
coherence  of  the  two  pump  beams,  the  mutually  pumped 
phase  conjugator  has  great  potentials  in  many  applica¬ 
tions  including  optical  communication,  laser  phase  lock¬ 
ing,  optical  interconnect,  optical  computing,  etc.  [3-6]. 
Experimentally,  it  has  been  observed  that  the  mutual 
coherence  of  the  two  pump  beams  changes  the  perfor¬ 
mance  of  the  mutually  pumped  phase  conjugator  con¬ 
siderably.  This  is  critical  in  applications  such  as  laser 
phase  locking,  and  coherent  optical  communication.  For 
example,  in  the  phase-locking  applications  the  two  free- 
runing  lasers  are  mutually  incoherent  before  they  are 
phase  locked.  These  two  laser  beams  are  incident  into  a 
photorefractive  mutually  pumped  phase  conjugator  as 
the  pump  beams.  Once  the  phase  conjugation  process 
occurs  and  if  the  frequency  of  the  lasers  are  colse  to  each 
other,  these  two  lasers  tend  to  be  phase  locked,  i.e.,  they 
become  partially  or  fully  coherent.  It  would  be  desirable 
to  maintain  an  efficient  coupling  through  the  mutually 
pumped  phase-conjugation  process  with  partially  or  fully 
coherent  laser  beams  so  that  lasers  can  be  stably  phase 
locked.  However,  it  has  been  observed  experimentally 


*  P.  Yeh  is  also  a  Principal  Technical  Advisor  at  Rockwell  Interna¬ 
tional  Science  Center 


that  the  performance  of  some  configurations  of  the 
photorefractive  MPPCs  is  quite  sensitive  to  the  degree  of 
mutual  coherence  of  the  pump  beams  [7].  In  previous 
theoretical  models,  the  pump  beams  are  assumed  mutu¬ 
ally  incoherent  and  only  the  transmission  gratings  are 
considered.  In  fact,  when  the  two  pump  beams  are  par¬ 
tially  coherent,  reflection  and  2^-reflection  gratings  must 
be  taken  into  account.  Although  multiple  grating  effects 
in  the  conventional  four-wave  mixing  geometry  have 
been  investigated  theoretically,  the  results  are  unable  to 
explain  the  experimental  observations  due  to  the  geo¬ 
metry  differences  and  various  limitations  [8, 9].  In  this 
paper,  we  present  theoretical  studies  of  the  mutuallv 
pumped  phase-conjugation  process  with  partial!) 


Fig.  1.  a  Schematic  of  a  photorefractive  mutually  pumped  ph^NC 
conjugator  with  the  presence  of  transmission  (T),  reflection  <  ft  i 
and  2k  reflection  (2A:,  2k')  gratings,  b  Schematic  drawing  of  a  bird 
wing  MPPC 
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coherent  bearns.  The  model  takes  into  account  trans¬ 
mission,  reflection,  and  2/r-reflection  gratings  as  well  as 
the  depletion  of  the  pump  beams.  The  results  are  in 
quanlitative  agreement  with  the  experimental  results 
[7,10]. 

It  is  well  known  that  the  transmission  gratings  are 
responsible  for  the  mutually  pumped  phase-conjugation 
process  in  most  of  the  configurations  [  1 1].  For  simplicity, 
we  consider  the  case  of  one  interactive  region  in  the 
mutually  pumped  phase-conjugation  process,  also 
known  as  Double-Phase  Conjugate  Mirror  (DPCM).  We 
assume  the  diffusion-dominated  photorefractive  effect, 
i.e.,  with  real  coupling  constants.  Referring  to  Fig.  la,  we 
consider  two  plane  waves  of  amplitude  of  A  2  and  A  4.  with 
a  degree  of  mutual  coherence  of  v(0<  v<  1)  incident  in 
a  photorefractive  crystal.  When  the  pump  beams  are 
partially  coherent,  besides  the  transmission  grating 
/li/lj  + /I J/4 3,  reflection  grating  AiA^  +  AfA^,  2k  grat¬ 
ings  and  A3  At  are  all  present  in  the  photorefrac¬ 

tive  medium.  The  nonlinear  coupled  wave  equations  are 
written  as: 
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and  the  boundary  conditions  are 
,  1^4(0)!' 

.4,(0)  =  .43(1)  =  0,  and 


(1) 


(2) 


where  lo  =  \Ai\^  +  \A2\^  +  \A3\^  +  \A^\^  is  the  total  in¬ 
tensity  of  all  the  beams,  y,,  y,,  72k ,  and  y2*'  are  the 


photorefractive  coupling  coefficients  for  transmission, 
reflection,  and  2/c-reflection  gratings,  respectively,  ^  irthe 
intensity  ratio  of  the  pump  beams,  L  is  the  interaction 
length.  It  is  clear  that  when  two  pump  beams  are  mutu¬ 
ally  incoherent,  i.e.,  v  =  0,  the  coupled  wave  equations 
reduce  to  the  conventional  ones,  where  only  the  trans-'v 
mission  grating  is  present  [11]. 

Since  strong  beam  couplings  occur  in  photorefractive 
materials,  pump  depletion  is  common  in  most  of  the 
photorefractive  phase-conjugation  processes,  especially 
in  MPPC.  Therefore,  we  have  to  solve  the  nonlinear 
coupled  wave  equation  taking  into  account  the  depletion 
of  the  pump  beams.  The  coupled  wave  equations  are 
solved  with  numerical  techniques.  We  use  the  conven¬ 
tional  4-th-order  Runge-Kutta  algorithm  to  solve  the 
coupled  equation  in  (1)  first.  We  then  fit  the  correct 
boundary  conditions  (2)  and  allow  the  program  to  iterate 
until  the  boundary  conditions  are  satisfied.  It  is  found 
that  the  iteration  usually  converges,  although  it  may  t.ike 
much  more  iterations  to  converge  for  some  parameter 
regions  than  the  other.  We  investigate  the  transmission 
properties  of  the  mutually  pumped  phase  conjugaior  ..is 


Fig.  2.  Transmissions  of  a  photorefractive  mutually  pumpcO  r*'  • 
conjugator  as  a  function  of  mutual  coherence  of  the  isso 
beams,  where  y,Z,=2.2,  y,L=»0.7,  y2kL=  1.0,  y2*'f-=  -  1  0  < 


Fig.  3.  Transmissions  of  a  photorefractive  mutually  pumpeO  ; ' .  < 
conjugator  as  a  function  of  mutual  coherence  of  the  two  p  •  * . 
beams,  where  y,L»«2.2,  y,L“0.3,  yjtL-O.S,  y2»/L»0.5.  v  '  ' 
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a  function  of  the  degree  of  mutual  coherence  of  the  pump 
beams.  The  transmissions  are  defined  as  the  ratio  of  the 
phase-conjugate  intensity  and  the  pump  intensity,  i.e., 
Tf  =  \Ai{L)\^l\A^{QI)\^  for  the  forward  transmission, 
and  =  MafOjiVi for  the  backward  trans¬ 
mission.  It  has  been  shown  theoretically  that  for  the 
mutually  pumped  phase  conjugators  that  have  only 
transmission  gratings,  the  forward  and  the  backward 
transmissions  are  always  identical  [11].  This 
phenomenon  can  be  interpreted  by  using  the  principle  of 
reciprocity  [12].  However,  in  the  presence  of  reflection 
gratings,  one  of  the  distinct  characteristics  is  that  these 
two  quantities,  Tf  and  r,,,  are  no  longer  the  same.  Figure 
2  shows  the  variation  of  the  transmissions  as  functions 
of  the  degree  of  mutual  coherence  with  ■/,£,  =  2.2, 
=  yikL=\.0,  y2k'L=  1.0,  and  q=  1.  We  note  that 

these  two  transmissions  have  a  quite  different  depen¬ 
dence  on  the  degree  of  mutual  coherence.  Depending  on 
the  sign  and  the  magnitude  of  the  reflection-grating 
coupling  strength  {y,L,  yi^'L),  the  performance  of 

the  mutually  pumped  phase  conjugators  can  be  signifi- 


Fig.  4.  Transmissions  of  a  photorefractive  mutually  pumped  phase 
conjugator  as  a  function  of  mutual  coherence  of  the  two  pump 
beams,  where  y^L-1,2^  72*^“  ”0.5,  72*'^“  ”0.5, 

^=1,0 


Fig.  5.  Coupling  constant  via  Ik  grating  in  the  BaTi03  crystal  as 
a  function  of  the  angle  between  the  pump  beam  and  the  c-axis 


cantly  enhanced  or  reduced  due  to  the  mutual  coherence 
of  the  pump  beams.  Figure  3  shows  an  example  where 
the  phase  conjugator  has  a  weak  reflection  coupling 
(}',L  =  0.3),  and  strong  favored  2k  couplings  (y2t^  =  0.5, 
yi^’L-0.5).  All  other  parameters  remain  the  same  as 
those  in  Fig.  2.  Due  to  the  contribution  of  the  2k  gratings 
where  the  energy  is  coupled  towards  the  directions  of 
phase  conjugation,  the  phase  conjugator  becomes  more 
efficient  when  both  beams  are  coherent.  On  the  other 
hand,  if  we  change  the  sign  of  the  2^  coupling  coef¬ 
ficients,  the  phase  conjugator  will  cease  oscillating  when 
the  pump  beams  become  coherent.  This  example  is  illus¬ 
trated  in  Fig.  4.  We  find  that  the  degree  of  mutual 
coherence  between  two  pump  beams  could  significantly 
affect  the  performance  of  the  mutually  pump  phase  con¬ 
jugators  especially  when  the  crystal  has  strong  reflection¬ 
coupling  coefficients. 

In  a  practical  system,  the  coupling  directions  and  the 
coupling  strength  of  the  reflection  and  2k  gratings  are 
mainly  dependent  on  the  specific  configuration  used.  We 
now  discuss  the  coherence  effect  of  some  of  the  most 
common  mutually  pumped  phase-conjugator  geome¬ 
tries.  For  the  case  of  DPCM  geometry  and  the  bridge 
phase-conjugator  geometry  in  BaTi03  crystals,  two 
pump  beams  are  usually  incident  from  two  sides  of  the 
crystal  with  a  small  intersection  angle.  The  c-axis  of  the 
crystal  has  an  acute  angle  with  one  phase-conjugate 
beam  and  an  obtuse  angle  with  the  other  phase-conjugate 
beam  [1, 13].  As  a  result,  the  2/c-coupling  coefficients  in 
the  directions  of  the  two  pump  beams  always  have  op¬ 
posite  signs  in  these  configurations.  This  fact  becomes 
more  clear  if  we  plot  the  coupling  constant  of  a  2A 
grating  in  BaTiOs  crystals,  as  shown  in  Fig.  5  [3].  From 
Fig.  5  we  notice  that  in  BaTiOj  crystals  the  energy  is 
coupled  toward  the  positive  c-axis  via  the  2k  grating 
occurs  only  in  a  relatively  narrow  range  of  angles,  (i.e  . 
a  <  6”).  In  ploting  Fig.  5  we  assumed  the  trap  density  of 
the  BaTiOs  crystal  .A^a  =  8  x  10‘*/cni’.  Furthermore, 
since  the  two  pump  beams  in  the  DCPM  geometry  are 
directly  overlapped  with  each  other,  the  reflection  grating 
formed  by  the  interference  of  the  pumps  is  also  very 
strong.  Therefore,  according  to  our  analysis,  it  is  expect¬ 
ed  that  the  DPCM  and  the  bridge  phase-conjugator  (in 
BaTiOs)  geometries  arc  very  susceptible  to  the  mutual 
coherence  of  the  pump  beams.  Their  performance  would 
generally  decrease  as  the  pump  beams  become  coherent 
This  is  similar  to  the  case  shown  in  Fig.  2.  For  the  Bird- 
Wing  Phase  Conjugator  (BWPC)  [14]  and  the  Mutually 
Incoherent-Beam  Coupler  (MIBC)  [15]  configurations, 
the  performance  could  be  much  less  susceptible  to  the 
coherence  than  in  DPCM  configuration  and  should  even 
be  enhanced  with  the  coherent  pump  be.ms.  The  reason 
is  that  in  these  two  geometries,  pump  •'earns  are  not 
directly  overlapped  in  space,  instead  they  are  intercon¬ 
nected  by  the  fanning  via  total  internal  reflections. 
Therefore,  they  usually  have  weak  reflection  gratings 
written  by  the  pump  beams.  Most  importantly,  they  hav  e 
the  same  coupling  directions  (same  sign)  for  the  2k  grat¬ 
ings  and  the  energy  coupling  is  also  towards  the  phase- 
conjugate  directions.  Taking  the  BWPC  in  the  BaTiOj 
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crystal  for  example,  as  shown  in  Fig.  lb,  two  pump 
beams  are  incident  in  the  two  sides  of  the  a-face  of  the 
crystal,  where  the  positive  c-axis  is  codirectional  with 
both  pump  beams.  From  the  characteristics  of  the  Ik- 
coupling  coefficient  in  the  BaTi03  crystals  as  shown  in 
Fig.  5,  we  know  that  once  the  angle  between  the  2A:-grat- 
ing  wavevector  and  the  c-axis  direction  is  larger  than 
approximately  6°.  the  energy-coupling  direction  is  always 
opposite  to  the  c-axis.  The  two  pump  beams  have  angles 
larger  than  6“  with  respect  to  the  c-axis  in  both  the 
bird-wing  and  the  MIBC  configurations.  In  these  two 
geometries,  the  phase  conjugate  transmissions  or  reflec¬ 
tivities  should  be  increased  when  pump  beams  become 
mutually  coherent.  This  is  similar  to  the  case  as  shown 
in  Fig.  3.  It  is  also  reasonable  to  predict  that  the  MIBC 
geometry  will  gain  more  than  the  bird-wing  because  the 
MIBC  geometry  tends  to  have  longer  beam  path  inside 
the  crystal  thus  a  higher  gain-length  product  for  the  Ik 
gratings.  Although  our  theory  is  based  on  the  one- 
interactive-region  model,  two-interactive-region  geome¬ 
tries  will  have  similar  behaviors  since  they  may  be  treated 
as  a  generalized  version  of  the  one-region  process.  These 
results  are  also  consistent  with  the  experimental  measure¬ 
ment  [7,  10].  For  the  case  of  complex  coupling  constants, 
the  dynamic  behavior  of  the  phase  conjugator  could  be 
much  more  complicated  especially  when  the  conjugator 
is  seeded.  It  may  be  very  unstable,  even  chaotic  in  these 
situations.  More  detailed  studies  and  comparison  be¬ 
tween  the  theory  and  the  experiment  are  under  current 
investigation. 

In  conclusion,  we  have  developed  a  theoretical  model 
for  mutually  pumped  phase  conjugation  with  partially 
coherent  pump  beams.  The  analysis  shows  that  the  per¬ 
formance  of  the  phase  conjugator  can  be  decreased  or 


enhanced  depending  on  the  contribution  of  the  reflection 
gratings.  This  also  provides  us  a  guide  line  in  takin_ 
advantage  of  the  reflection-grating  couplings  for  en 
hanced  phase  conjugators. 
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We  investigated  2^-grating-assisted  self-pumped  phase  conjugation  in  photorefractive  CeiBaTiOs  crystals. 
The  phase-conjugation  process  involves  a  combination  of  four-\vave  mixing  and  stimulated  photorefractive 
backscattering.  An  approximation  involving  separate  interaction  regions  is- used  tt)  theoretically  calculate 
the  reflectivity  of  phase  conjugation  as  a  function  of  the  coupling  strength  of  four-wave  mixing  and  stimulated 
photorefractive  backscattering.  In  our  experiments,  grating-erasure  techniques  are  employed  at  the  inter¬ 
action  regions  to  investigate  the  dependence  of  phase-conjugate  reflectivity  on  the  coupling  strength  of  four- 
wave  mixing  and  stimulated  photorefractive  backscattering.  The  experimental  results  are  in  good  agreement 
with  the  theoretical  prediction. 

Key  words:  Self-pumped  phase  conjugate,  stimulated  photorefractive  backscattering,  four  wave  mixing. 
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1.  INTRODUCTION 

In  self-pumped  phase  conjugators  the  counterpropagating 
beams  needed  in  a  four-vt^ave-mixing  process  are  gener¬ 
ated  by  the  incident  beam  itself.  By  a  process  that  is 
not  entirely  understood,  the  self-generated  pump  beams 
are  often  formed  in  a  configuration  tending  to  optimize 
the  phase-conjugated  beam.  Self-pumped  phase  conju¬ 
gators  have  been  demonstrated  in  a  variety  of  configu¬ 
rations.  Two  classes  of  self-pumped  phase  conjugator 
have  been  proposed  and  demonstrated:  the  passive 
self-pumped  phase  conjugator,  which  needs  externally 
aligned  mirrors,^  and  the  true  self-pumped  phase  con¬ 
jugator,  which  uses  a  single  cr>^stal  without  external 
elements.^  The  true  self-pumped  phase  conjugator  is  of 
potential  interest  because  of  its  convenience  for  practical 
applications.^  The  first  true  self-pumped  configuration 
was  observed  in  1982  by  Feinberg  in  a  single  BaTiOa 
photorefractive  crystal.^  There  are  three  different  mod¬ 
els  to  explain  the  origin  of  self-pumped  phase  conjugators 
in  a  signal  photorefractive  crystal  (as  shown  in  Fig.  1). 
In  1983,  based  on  the  experimental  observation  of  the 
beam  path  in  BaTiOa  crystals,  MacDonald  and  Feinberg® 
proposed  a  model  that  involved  both  four-wave  mixing 
and  total  internal  reflection  (FWM-TIR)  [see  Fig.  1(a)]. 
In  the  model  the  phase  conjugation  is  a  result  of  FWM 
in  two  regions  that  are  connected  by  TIR’s  from  a  cor¬ 
ner  of  the  crystal.  Hence,  a  closed  optical  path  loop  is 
often  observed  inside  the  crystal.  In  1985,  Chang  and 
Hellwarth®  were  the  first  to  point  out  the  close  analogy 
between  backward-stimulated  Brillouin  scattering  and 
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some  configurations  of  self-pumped  phase  conjugation 
(SPPC)  in  BaTiOa  aind  suggested  that  stimulated  photo- 
refractive  backscattering  (SPB)  or  stimulated  two-wave 
mixing  may  be  responsible  for  the  generation  of  the  SPPC 
[as  shown  in  Fig.  Kb)].  It  has  been  shown  theoretically 
that  the  backward-scattered  wave  is  dominated  by  the 
phase  conjugate  of  the  incident  beam.'  ®  Recently  we®*^® 
proposed  an  alternative  model  of  FWM-SPB  (four-wave 
mixing  and  SPB,  in  'which  the  SPB  is  produced  by  the 
formation  of  2k  gratings)  to  explain  the  SPPC  observed  in 
KTN:Fe  and  BaTiOaiCe  crystals  [as  shown  in  Fig.  1(c)]. 
In  this  case  the  SPPC  formation  relies  on  both  four  FWM 
and  SPB  interactions.  The  phase-conjugate  beam  is  gen¬ 
erated  by  a  four-wave-mixing  process  invoMug  the  inci¬ 
dent  beam,  the  forward-propagating  beam  (the  fanning 
beam),  and  its  backward-stimulated  scattering  beam. 
The  backward-stimulated  scattering  beam  is  generated 
by  the  SPB  process  invol'ving  2k  gratings.  Hence  instead 
of  a  closed  loop  inside  the  crystal,  filaments  representing 
counter  propagating  beams  are  often  observed. 

In  a  t5q)ical  SPPC  experimental  configuration,  an  ex¬ 
traordinarily  polarized  input  beam  is  incident  on  the  a 
face  of  a  regular-cut  BaTiOs  photorefractive  crystal.  The 
SPPC  formation  mechEinism  is  usually  determined  by  the 
fanning  pattern  and  the  boimdary  conditions.  Gener¬ 
ally  speaking,  SPPC  is  often  generated  by  the  SPB  or 
the  FWM-SPB  mechanisms  in  some  doped  crystals  but 
by  the  FWM-TIR  mechanism  in  undoped  crystals.  This 
is  mainly  due  to  a  significantly  stronger  fanning  an 
larger  coupling  coefficients  for  2k  gratings  in  doped  crys¬ 
tals.  Changing  the  dopant  concentration  or  the  operat- 
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Fig.  1.  Different  modes  of  SPPC;  (a)  FWM-TIR,  (b)  SPB,  and 
(c)  FWM-SPB. 

ing  wavelength  can  cause  a  reconfiguration  of  the  initial 
fanning  pattern  that  leads  to  a  change  of  the  physical 
mechanism  and  the  reflectivity  of  SPPC.  Recently,  we 
observed  that  the  formation  mechanism  of  SPPC  was 
transformed  firom  FWM-SPB  to  FWM-TIR  when  the 
wavelength  of  the  input  beam  increased  or  the  dopant 
concentration  decreased.  We  also  observed  that  the 
mechanism  transformation  during  SPPC  and  an  enhance¬ 
ment  of  phase-conjugate  reflectivity  could  be  obtained  by 
the  variation  of  the  polarization  state  of  the  input  beam.^^ 
In  a  typical  experiment  an  extraordinarily  polarized  in¬ 
put  beam  produces  SPPC  by  FWM-SPB  in  a  CeiBaTiOa 
crystal.  As  the  input  polarization  varies,  the  ordinary 
component,  which  serves  as  an  erasing  beam,  decreases 
the  effective  coupling  constant.  This  leads  to  a  weaker 
fanning  pattern.  Depending  on  the  boundary  conditions, 
the  fanned  beam  may  reach  the  comer  of  the  crystal  cube, 
leading  to  SPPC  by  FWM-TIR.  In  this  paper  we  theo¬ 
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retically  and  experimentally  study  the  role  of  the  2k 
grating  of  SPB  and  the  coupling  strength  of  in  the 
FWM-SPB  mode  of  self-pumped  phase  conjugators.  A 
two-interaction-region  approximation  is  used  to  theoreti¬ 
cally  calculate  the  reflectivity  of  phase  conjugation  and 
the  threshold.  We  then  present  the  results  of  our  experi¬ 
mental  studies,  which  validate  the  theoretical  predictions. 
To  vary  the  coupling  constant,  we  also  use  an  ordinarily 
incoherent  erasing  beam  to  reduce  the  index  grating  in 
the  FWM  and  the  SPB  regions.  The  theoretical  and  the 
experimental  results  show  that  the  reflectivity  of  phase 
conjugation  depends  strongly  on  the  coupling  strength  of 
both  the  FW^I  and  the  SPB  regions.  In  addition,  the 
temporal  dynamics  of  the  phase-conjugate  signal  buildup 
process  indicates  that  the  2k  grating  plays  an  important 
role  in  the  FWM-SPB  model.  The  paper  can  be  orga¬ 
nized  as  follows.  In  Section  2  we  discuss  a  theoretical 
model  of  this  new  mode  of  SPPC.  A  closed-form  solution 
of  coupled-wave  equations  is  obtained  by  use  of  an  ap¬ 
proximation.  In  Section  3  we  discuss  the  experimental 
investigation  using  grating-erasure  techniques  at  the  two 
interaction  regions  to  validate  the  theoretical  predictions 
regarding  the  role  of  the  2k  grating  and  the  coupling 
strength  of  FWM  during  SPPC. 

2.  THEORETICAL  ANALYSIS 

Figure  2(a)  shows  a  photograph  of  a  typical  optical  beam 
path  in  the  FW^I-SPB  mode  of  SPPC  in  CeiBaTiOa  crys¬ 
tals.  Based  on  our  experimental  observation,  the  model 
of  the  FWM-SPB  phase  conjugator  can  be  divided  into  two 
interactions  regions,  i.e.,  SPB  and  FWM  interaction  re¬ 
gions  [as  shown  in  Fig.  2(b)].  In  the  SPB  interaction  re¬ 
gion  there  are  two  counterpropagating  beams  coupled  by 
self-induced  2k  gratings.  The  stimulated-backscattering 
beam  A2  is  initially  seeded  by  backscattering  of  fanned 
light  from  imperfection  (pits)  at  the  surface  of  the  crys¬ 
tal.  The  backscattered  light  is  amplified  by  SPB  interac¬ 
tion  with  fanning  beam  Ai .  These  two  beams  provide  the 
counterpropagating  pump  beams  needed  in  the  FWM  for 
the  generation  of  a  phase-conjugate  wave.  Thus  the  SPB 
interaction  region  mainly  acts  as  a  feedback  mirror  (the 
SPB  mirror,  defined  as  the  combination  of  the  2k  grating 
gain  region  plus  the  scattering  center)  that  provides  one  of 
the  pump  beams  for  the  FWM  interaction  region.  In  the 
FWM  interaction  region,  incident  beam  A3,  fanning  beam 
Al,  and  backscattering  beam  A2  are  coupled  by  a  FWM 
interaction  to  generate  a  phase-conjugated  beam  A4.  For 
the  phase-conjugation  process  to  be  sustained,  the  reflec¬ 
tivity  of  SPB  region  must  be  large  enough  to  generate  a 
self-oscillation  between  the  SPB  mirror  and  the  FWM  in¬ 
teraction  region.  Thus  this  phase  conjugator  is  similar 
to  a  semilinear  phase  conjugator,^®  in  which  the  mirror  is 
replaced  by  a  SPB  mirror.  The  threshold  and  the  reflec¬ 
tivity  of  a  FWM-SPB  phase  conjugator  are  related  to  the 
net  reflectivity  of  the  SPB  mirror. 

We  now  treat  this  phase  conjugator  by  using  both  op¬ 
tical  FWM  auid  contradirectional  two-wave  mixing  in  two 
separate  regions.  Referring  to  Fig.  2(b),  we  consider  the 
interactions  of  all  beams  by  using  the  following  two  sets 
of  coupled  equations.  We  consider  only  the  transmission 
grating  in  the  FWM  region  and  the  2k  grating  in  the 
SPB  region.  For  the  FWM  interaction  region  the  coupled 
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Fig.  2.  FWM‘SPB  phase  conjugator  with  a  two-interaction  re¬ 
gion:  (a)  photograph  of  the  optical  beam  path  in  a  CerBaTiOa 
FWM-SPB  phase  conjugator;  (b)  theoretical  model  in  which  the 
forward-going  fanning  beam  Ai  overlaps  well  with  its  backward¬ 
going  scattering  beaun  A2  to  form  a  SPB  region.  The  amplified 
backscattering  beam  A2  enters  the  FWM  regions,  mixes  with  Ai 
and  A4,  and  generates  the  phase  conjugate  A3. 


equations  are  written  as 

^  =  ^{A*A,  +  A2A*)A3  -  ^  AtAzAi, 
az  ZIq  ^io 

^  =  ^(A,At  +  A2A3)A2, 
dz  2/0 
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and  for  the  SPB  interaction  region,  the  coupled  equations 
are  written  as 


t 

s 


dAi 

dz' 


-gA.A5A2. 


dA2 

"dF 


At  A2A,, 


(2) 


where  yi  and  72  are  the  coupling  constants  of  the  trans¬ 
mission  grating  in  the  FWM  region  and  the  2k  grating 
in  the  SPB  region,  respectively..  They  are  real  and  posi¬ 
tive  without  the  external  electrical  field  in  diffusion- 
dominant  crystals,  such  as  BaTiOa,  SBN,  and  KTN 
crystals.  Iq  =  is  intensity  in  the 

FWM  region,  and  1q{z')  =  =  Z;-i  is 

the  total  intensity  in  the  SPB  region.  The  coordinate 
axis  of  the  FWM  region  is  defined  as  2,  and  that  of 
the  SPB  region  is  defined  as  2'.  These  two  regions  are 
connected  at  2'  =  0  or  2  =  Zi,  where  h  is  the  interac¬ 
tion  length  of  the  FWM  region  and  I2  is  the  interaction 
length  of  the  SPB  region.  In  a  typical  FWM-SPB  mode 
of  self-pumped  phase  conjugation  [as  shown  in  Fig.  2(a)] 
the  interaction  length  of  the  FWM  region  (Zi)  is  much 
smaller  than  that  of  the  SPB  region  (Z2),  but  the  yi  is 
much  larger  than  72  (in  our  case,  I2  ^  6Z1  and  71  672). 

Hence  these  two  sets  of  equations  can  be  analytically 
solved  by  neglecting  the  2k  grating  in  the  FWM  region. 
In  what  follows,  we  first  obtain  the  reflection  gener¬ 
ated  by  the  SPB  interaction  by  solving  Eqs.  (2).  We 
then  solve  the  coupled-mode  equations  in  the  FWM  re¬ 
gion  by  neglecting  the  2Zj -grating  terms  of  Eq.  (1).  The 
phase-conjugate  threshold  and  reflectivity  of  the  FWM- 
SPB  phase  conjugator  are  obtained  in  terms  of  the  cou¬ 
pling  strength  (72Z2)  of  SPB  and  that  (71Z1)  of  FWM 
interaction. 


A.  Stimulated  Photoreffactive  Backscattering 
Interaction  Region 

To  examine  the  reflectivity  of  the  SPB  mirror,  we 
first  solve  the  coupled  equations  [Eqs.  (2))  by  follow¬ 
ing  the  procedure  of  contrapropagating  two-wave  mix¬ 
ing  in  Ref.  16  and  by  using  the  boimdary  condition 
h{h)/h{l2)  =  ^0  at  the  surface  of  crystal  (at  2'  =  Z2), 
where  mo  is  the  backscattering  coefficient  of  the  imperfec¬ 
tion  on  the  surface  of  crystal.  The  solutions  of  Eqs.  (2) 
are 


7i(20  =  1Ai(2')I^  =  -C  +  +  B  exp(-7220, 

hiz’)  =  |A2(20I^  =  C  +  VC2  +S  exp(-722'),  (3a) 
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where  B  and  C  are  integration  constants  given  by 


whe 

net 


B  =  /i(Z2)/2(Z2)exp(-72Z2),  C  =  ^^(^2)  ^  •  (3b)  ^ 

i 

The  effective  intensity  reflectivity  of  the  SPB  mirror,  R,  | 
at  2'  =  0  (or  2  =  Zi)  can  therefore  be  written  w 


and 

equi 

the 


(1) 
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Fig.  3.  Reflectivity  of  the  SPB  mirror  as  a  function  of  yzh  with 
various  rriQ. 


tanh  s  =  s .  (S) 

4 

Equations  (7)  and  (8)  provide  some  constraints  on  yiZi  and 
s.  We  note  that  nontrivial  solutions  of  the  transcenden¬ 
tal  equation  exist  only  when  yili  >4.  Once  s  is  solved,  a 
is  determined  by  Eq.  (7).  Being  a  new  power  flow,  cr  has 
to  be  real.  This  requires  >  1/(1  +  H).  Thus  if  s  or 
yili  are  too  small,  there  would  be  no  phase  conjugation. 
These  constraints  lead  to  a  threshold  coupling  strength 
(yi^i)th  and  a  threshold  SPPC  reflectivity  77th  of  the  FWM- 
SPB  phase  conjugator,  according  to  Eqs.  (6)  and  (8),  of 


(9) 

(10) 


(1  -  mp)^  +  2mo  exp(y2^2)  -  (1  ~  rno)y/(l  -  mp)^  +  4mo  exp(y2^2) 
2mo  exp(y2^2) 


When  mo  «  1,  the  reflectivity  R  can  be  written  approxi¬ 
mately  as 


1  +  2mo  exp(y2/2)  ~  Vl  +  imp  exp(y2^ 
2mo  exp(y2^2) 


We  note  that  R  depends  strongly  on  the  product  of  mo 
and  exponential  factor  exp(72^2)-  According  to  Eq.  (5), 
the  presence  of  an  efficient  2k  grating  (72/2  »  1)  is  suf¬ 
ficient  to  achieve  a  high  reflectivity  even  through  the  ini¬ 
tial  seeding  mo  from  crystal-surface  imperfection  is  very 
weak.  The  reflectivity  of  the  SPB  mirror  reaches  100% 
when  72^2  approaches  infinity.  By  use  of  Eq.  (4),  Fig.  3 
plots  the  reflectivity  of  the  SPB  mirror  as  a  function  of 
72/2  for  various  mo-  It  can  be  seen  that  for  a  given  small 
72/2,  the  reflectivity  R  of  the  SPB  mirror  is  very  different 
for  different  mo.  When  72/2  increases,  the  difference  is 
diminished  and  R  reaches  almost  100%  when  72/2  >  5, 
while  mo  varies  from  10'®  to  10”  ^  This  means  that  the 
amplification  of  a  2k  grating  provides  an  efficient  mecha¬ 
nism  to  retroreflect  the  fanning  beam  into  the  FWM  re¬ 
gion  to  facilitate  the  phase-conjugation  process. 


B.  Four-Wave-Mixing  Interaction  Region 
With  the  reflectivity  of  the  SPB  mirror  available  we  can 
now  solve  the  coupled  equations  (1)  in  the  FWM  region  by 
neglecting  the  2k  grating.  Following  a  procedxire  similar 
to  the  one  used  in  Ref.  1,  we  obtain  the  following  expres¬ 
sion  for  the  phase-conjugate  reflectivity: 


|A3(2  =  0)P  _  (  ± 

|A4(z  =  0)P  r’ 


(6) 


where  R  again  is  reflectivity  of  the  SPB  mirror,  cr  is  the 
net  power  flow  in  this  region,  given  by 


and  s  can  be  obtained  firom  the  following  transcendental 
equation,  which  is  related  to  the  coupling  strength  71  Zi  of 
the  FWM  region: 


We  note  that  the  phase-conjugate  reflectivity  is  finite  at 
threshold.  Equation  (9)  shows  that  (7i/i)th  is  a  decreas¬ 
ing  function  of  R,  reaching  a  minimum  of  4.98  at  i?  =  1. 
In  order  to  further  understand  the  concept,  we  plot  (7i/i)th 
as  a  function  of  72/2  for  various  scattering-center  reflec¬ 
tivities  mo  (shown  in  Fig.  4).  The  SPPC  reflectivity  at 
threshold  77th  is  also  plotted  as  a  function  of  72/2  for  vari¬ 
ous  mo  (shown  in  Fig.  5).  It  can  be  seen  that  both  mo 
and  72/2  influence  greatly  the  threshold  of  the  FWM- 
SPB  phase  conjugator.  The  threshold  value  of  (yi/Oth  de¬ 
creases  with  mo  and  72/2-  The  minimum  value  of  4.98  is 
almost  the  same  as  the  threshold  of  the  FWM-TIR  phase 
conjugator  (4.68  without  loss,  in  which  we  should  note 
that  there  is  a  factor  of  2,  according  to  our  definition  of 
the  coupling  constant).  It  is  also  interesting  to  note 
that  the  reflectivity  of  SPPC  at  the  threshold  77th  is  a  de¬ 
creasing  function  of  (7ili)th-  In  general,  the  reflectivity 
of  scattering  center  mo  in  the  crystals  is  —10“^.  An  ef¬ 
fective  method  to  obtain  a  large  SPB-mirror  reflectivity  R 
is  to  increase  the  coupling  strength  of  the  2k  grating  72/2 


Fig.  4,  Threshold  coupling  stren^h  (7i/i)th  of  the  FWM-SPB 
phase  conjugator  versus  the  coupling  strength  of  the  2k  grating 
72/2  for  various  mo. 
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Fig.  5.  SPPC  reflectivity  at  threshold  versus  72^2  for  vari¬ 
ous  mo. 
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Fig.  6.  Phase-conjugate  reflectivity  versus  coupling  strength  of 
the  2k  grating  72^2'  (a)  at  different  coupling  strengths  of  the 
transmission  grating  71/1  and  (b)  for  various  mo. 


(as  shown  in  Fig.  3),  This  can  be  achieved  by  properly 
doping  the  crystals  (e.g.,  CerBaTiOo). 

For  7i/i  >  (7i/i)th  the  reflectivity  of  the  FWM-SPB 
phase  conjugator  77  can  be  plotted  as  a  function  of  72/2  for 
the  various  yih  and  mo  [as  shown  in  Figs.  6(a)  and  6(b)]. 
We  note  that  the  reflectivity  77  obviously  increases  with 
72/2  and  behaves  ver>^  differently  for  different  71Z1.  The 


reflectivity  approaches  a  saturation  value  for  a  given  71/1, 
corresponding  to  a  SPB-mirror  reflectivity  of  R  =  100% 
(at  72Z2  =  Referring  to  Fig.  3,  we  recall  that  SPB- 
mirror  reflectivity  R  is  an  increasing  function  of  72/2, 
reaching  R  =  100%  when  72/2  =  **  We  also  note  that  all 
curves  of  different  mo  will  approach  the  same  saturation 
value  when  72/2  ~  These  figures  confirm  that  the 
presence  of  a  SPB  2k  grating  is  essential  to  enhance 
the  phase-conjugate  reflectivity  during  SPPG.  In  order 
to  increase  the  SPB-mirror  reflectivity  R,  we  can  insert 
a  retroreflecting  screen  into  the  path  of  the  beam  after 
the  fanning  beam  leaves  the  crystal  or  use  an  artificial 
scattering  center  on  the  crystal  face  to  increase  mo.  We 
can  also  increase  the  coupling  constant  of  2k  gratings  by 
choosing  an  appropriate  dopant  and  a  proper  wavelength. 
However,  when  the  reflectivity  of  backscattering  center 
mo  is  large  enough,  the  phase-conjugate  reflectivity  is 
insensitive  to  the  coupling  strength  of  a  2k  grating.  For 
example,  in  the  case  of  FWM-TIR,  the  comer  cube  is 
already  providing  a  near  100%  retroreflection  regardless 
of  the  presence  of  2k  gratings.  A  small  decrease  of  phase- 
conjugate  reflectivity  of  —20%  was  reported  when  the  2k 
grating  was  erased  by  an  erasure  beam  in  Ref.  18. 

4.  EXPERIMENT 

In  this  section  we  report  the  results  of  our  experiments 
to  verify  theoretical  predictions  based  on  the  model  of  the 
FWM-SPB  phase  conjugator.  The  schematic  diagram 
of  the  experimental  arrangement  is  shown  in  Fig.  7. 
Our  Ce-doped  BaTiOa  crystal  has  a  red-orange  color 
and  a  dopant  concentration  of  approximately  30  parts 
in  10®.  The  crystal  is  poled  and  0“  cut  with  dimensions 
6.15  mm  x  5.20  mm  X  8.2  mm  (a  x  6  x  c).  An  extraor¬ 
dinarily  polarized  beam  from  an  argon  laser  (wavelength 
of  514.5  nm)  is  incident  upon  the  crystal  on  the  a  face 
at  an  external  angle  of  ^  =  60*.  The  phase  conjugation 
was  monitored  by  a  calibrated  photodiode  (PDl).  In  or¬ 
der  to  investigate  the  buildup  process  of  2k  gratings,  the 
fanning  beam  was  also  monitored  by  another  photodiode 
(PD2),  and  the  optical  path  in  the  crystal  was  observed  by 
imaging  through  the  top  of  the  crystal  into  a  CCD  camera. 
Figure  2(a)  shows  a  photograph  of  the  optical  beam  path 
inside  the  crystal  during  a  typical  mode  of  FWM-SPB 
SPPC.  The  input  beam  is  depleted  almost  completely 
after  propagating  through  a  short  distance  and  bends  to 
the  bottom  surface  of  the  crystal  owing  to  the  fanning 
effect.  At  steady  state  the  end  of  the  beam  is  charac¬ 
terized  by  a  filament  that  is  sustained  through  the  SPB 


Fig.  7.  Experimental  setup  used  to  demonstrate  SPPC  in  the 
FWM-SPB  mode:  PBS’s,  polarizing  beam  splitters;  BS’s,  beam 


splitters. 
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Fig.  8.  Time  evolution  of  the  reflectivity  of  SPPC  (solid  curve) 
and  intensity  of  the  fanning  beam  (dashed  curve). 


process  by  2k  gratings.  This  SPPC  formation  is  signi¬ 
ficantly  different  from  the  traditional  FWM-TIR  mecha¬ 
nism.  There  is  no  closed  optical  loop  involving  a  comer 
inside  the  crystal. 

Figure  8  shows  the  time  evolution  of  the  phase- 
conjugate  signal  (the  solid  curve)  and  the  fanning  sig¬ 
nal  (the  dotted  curve)  during  the  buildup  of  the  phase 
conjugation.  At  the  first  stage  (t  of  0-110  s,  where 
^  =  0  corresponds  to  the  time  when  the  incident  beam  is 
turned  on),  the  fanning  signal  grows  to  a  steady  state, 
but  the  phase-conjugation  signal  only  begins  to  start, 
as  indicated  by  the  recording  at  PDl.  It  indicates  that 
at  this  moment,  only  the  fanning  beam  is  produced  by 
the  fanning  grating  in  the  FWM  region,  but  the  2k  grat¬ 
ing  has  not  been  established  in  the  SPB  region.  The 
forward  fanning  beam  cannot  be  retroreflected  into  the 
FWM  region  to  generate  the  phase  conjugation.  When 
^  >  110  s,  the  phase-conjugation  signal  appears  suddenly, 
while  the  fanning  signal  decreases  significantly.  This  in¬ 
dicates  the  buildup  of  the  2k  grating.  After  the  fanning 
beam  is  produced  by  the  incident  beam,  a  backscatter- 
ing  can  be  generated  by  some  defects  on  the  surface  of 
crystal,  and  the  2k  grating  can  be  written  by  the  fanning 
beam  and  the  retroreflected  scattering  beam.  Once  the 
2k  grating  is  formed,  energy  will  be  transformed  from 
the  fanning  beam  to  the  backscattering  beam.  With  the 
increase  of  the  backscattering,  the  amplitude  of  2k  grat¬ 
ing  also  increases.  This  leads  to  more  energy  transfer 
from  the  fanning  beam  to  the  backscattering  beam.  The 
positive  feedback  provides  the  stimulated  photorefractive 
backscattering  and  leads  to  some  filaments  at  steady 
state.  At  the  same  time  the  backscattering  beam  will 
generate  the  phase  conjugation  by  the  FWM  region  and 
rewrite  a  new  in-phase  version  of  the  fanning  grating 
with  the  phase-conjugate  beam  to  enhance  the  coupUng 
in  the  FWM  region.  This  explains  the  sharp  rise  of  the 
phase  conjugation  as  indicated  in  Fig.  8. 

To  further  study  the  characteristics  of  the  2k  grating 
in  our  model,  we  use  an  ordinarily  polarized  beam  (see 
Fig.  7)  that  is  incoherent  with  the  input  beam  to  illumi¬ 
nate  the  filaments  uniformly  in  the  SPB  region.  Since 
the  refractive-index  change  of  a  photorefractive  crystal  is 
proportional  to  the  modulation  depth  of  the  input  light, 
the  coupling  coefficient  of  the  2k  grating,  yzih),  can  be 
written  to  incorporate  the  reduction  in  the  modulation 


that  results  from  the  erasing  beam, 


J'zih)  = 


yzo 

1  +  Ir/I  ' 


(11) 


where  I  is  the  total  intensity  in  the  filament,  and  y-zo 
is  the  coupling-coefficient  constant  in  the  absence  of  the 
erasure  beam  7^.  According  to  the  above  equation,  the 
coupling  strength  of  the  photorefractive  grating  can  be 
easily  controlled.  Figure  9  shows  the  time  evolution  and 
the  reflectivity  of  the  phase-conjugate  signal  as  the  era¬ 
sure  beam  is  turned  on.  The  signal  decreases  to  different 
steady-state  values  depending  on  the  intensity  of  erasure 
beams  /r.  Figure  9  also  shows  that  the  reflectivity  of 
the  phase-conjugation  signal  decreases  with  the  intensity 
of  the  erasing  beam  and  exhibits  a  sharp  threshold  at 
7^  =  4  W/cm^.  To  further  illustrate  the  concept,,  we  re¬ 
plot  the  phase-conjugate  reflectivity  as  a  function  of  yzhy 
in  which  the  intensity  of  the  erasure  beam  can  be  calcu¬ 
lated  in  terms  of  the  coupling  strength  by  use  of  Eq.  (10). 
We  also  present  the  theoretical  results  by  the  dashed  line 
in  the  figure  for  the  purpose  of  comparison.  The  results 
are  shown  in  Fig.  10,  in  which  we  estimate  an  average 
of  7  as  4  W/cm^,  a  value  of  yzoh  as  11.5,  a  value  of  yili 
as  7.0,  and  a  reflectivity  of  scattering  center  mo  as  0.001. 
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Fig.  9.  Time-evolution  of  the  phase-conjugate  signal  when  the 
different  erasing  beams  are  applied. 


Fig.  10.  Experimental  phase-conjugate  reflectivity  versus  yzh- 
Note  the  phase-conjugate  threshold  at  72/2  7.5  (correspond¬ 

ing  to  7r  =  4  W/cm^).  The  solid  curve  is  from  a  theoretical 
calculation. 
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Polarization  Angle  of  Input  beam  (<l>) 

Fig.  11.  Normalized  phase-conjugation  reflectivity  as  a  function 
of  the  polarization  angle,  <t>,  of  the  input  beam.  These 
experimental  data  are  not  corrected  for  Fresnel  reflections. 
The  dashed  curve  is  from  a  theoretical  calculation  for  a  pure 
FWM-SPB  mode  of  SPPC. 


This  curve  is  in  agreement  with  the  experimental  results. 
Thus  we  note  again  that  the  presence  of  a  2A-grating 
significantly  affects  the  reflectivity  of  the  backscattering 
beam  and  the  SPPC  in  this  model. 

In  a  different  experiment,  to  study  the  characteristic  of 
the  FWM  region  in  this  model,  we  block  the  erasing  beam 
and  introduce  a  general  linear  polarized  input  beam  by 
inserting  a  half-wave  plate  in  the  path  of  the  input  beam 
(as  shown  in  Fig.  7).  The  plane  of  polarization  of  the  in¬ 
put  beam  could  now  be  changed  by  rotating  the  half-wave 
plate.  This  leads  to  the  presence  of  an  ordinarily  polar¬ 
ized  component  of  the  input  beam.  It  is  known  that  the 
coupling  coefficient  in  CerBaTiOa  crystal  for  an  ordinarily 
polarized  beam  is  much  less  than  that  of  an  extraordinar¬ 
ily  polarized  beam,  and  the  cross  coupling  is  negligible. 
Thus  the  extraordinarily  polarized  beam  will  contribute 
significantly  to  the  fanning  and  the  phase  conjugation, 
whereas  the  ordinarily  polarized  beam  will  serve  only  to 
erase  the  gratings  in  the  interaction  regions.  In  this  case 
the  effective  coupling  coefficient  of  the  grating,  r  :(</>),  can 
be  written  as 


yio 

1  +  tan^  <l> 


(12) 


where  yio  is  the  coupling  coefficient  at  =  0®,  where  <i> 
is  the  state  of  polarization  of  the  input  beam.  It  is  de¬ 
fined  (as  shown  in  Fig.  7)  as  =  0®  for  the  extraordi¬ 
narily  polarized  direction  and  <^  =  90®  for  the  ordinarily 
polarized  direction.  The  experimental  (solid  curve)  and 
the  theoretical  (dashed  curve)  results  of  the  reflectivity  of 
phase  conjugation  are  shown  in  Fig.  11.  Here,  we  note 
that  the  phase-conjugate  output  is  always  extraordinar¬ 
ily  polarized  regardless  of  the  input  polarization  state. 
The  reflectivity  of  phase  conjugation  is  therefore  defined 
as  the  ratio  of  the  conjugated  signal  and  the  extraordi¬ 
nary  component  of  the  incident  beam.  In  other  words, 
the  reflectivity  is  normalized  to  that  of  a  totally  extraor¬ 
dinarily  polarized  state  of  the  input  beam  (</>  =  0®),  -Ac¬ 
cording  to  the  theoretical  prediction,  the  reflectivity  will 
decrease  with  the  ordinarily  component  of  input  beam. 
On  the  contrary,  the  experimental  results  show  an  im¬ 
provement  over  the  theoretical  prediction.  The  reflectiv¬ 
ity  of  phase  conjugation  exceeds  the  value  at  —  0®  in 


the  range  4>  <  35®.  The  increase  in  phase-conjugate  re¬ 
flectivity  can  be  caused  by  a  mechanism  transformation 
during  SPPC.^'*  In  addition,  the  presence  of  an  ordinar¬ 
ily  polarized  component  in  the  input  beam  can  restrain 
the  self-generated  fanning  pattern.  Thus  the  energy  loss 
caused  by  fanning  can  be  limited. 

5.  CONCLUSION 

We  have  carried  out  a  theoretical  analysis  of  a  FWM-SPB 
self-pumped  phase  conjugator  by  using  two-interaction- 
region  approximation.  The  phase-conjugate  wave  is 
generated  by  four-wave  mixing  with  the  assistance  of 
2k  gratings  by  SPB.  From  our  results,  we  find  that 
the  coupling  strength  72^2  of  the  2k  grating  in  the  SPB 
region  and  the  feedback  reflectivity  mo  of  scattering  cen¬ 
ters  play  very  important  roles  in  the  generation  of  the 
phase-conjugate  wave.  We  have  also  presented  the  re¬ 
sults  of  a  comprehensive  experimental  investigation  of  a 
FWM-SPB  phase  conjugator  by  using  two  different  types 
of  erasing  beam.  We  find  that  the  phase  conjugation  of 
the  FWM-SPB  model  relies  on  both  the  FWM  and  the 
SPB  interactions. 
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We  investigate  contradirectiona!  two-wave  mixing  with  partially  coherent  waves  in  photorefractive  crystals  in 
the  nondepleted  pump  regime.  Equations  governing  the  propagation  of  the  self-coherence  function  and  the 
mutual-coherence  function  of  the  signal  wave  and  the  pump  wave  are  derived  and  simulated  numerically. 
Numerical  solutions  of  these  equations  are  in  excellent  agreement  with  the  experimental  measurements. 
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Two-wave  mixing  in  photorefractive  crystals  is  an  in¬ 
teresting  and  useful  nonlinear-optical  phenomenon  for 
many  applications  such  as  image  amplification,  laser 
wave  cleanup,  spatial  light  modulators,  thresholding, 
and  power-limiting  devices.  Most  of  the  theoretical 
study  in  this  area  has  addressed  wave  mixing  with 
monochromatic  waves.  Two-wave  mixing  with  par¬ 
tially  coherent  waves  has  been  studied  for  the  case  of 
transmission  grating  interaction,^*^  in  which  the  opti¬ 
cal  path  difference  between  the  two  interacting  waves 
remains  approximately  the  same  as  the  two  waves 
propagate  codirectionally  through  the  photorefractive 
medium,  especially  when  the  incident  angles  of  the  two 
waves  are  close.  In  the  case  of  reflection  grating  in¬ 
teraction  the  optical  path  difference  between  the  two 
interacting  waves  varies  significantly  as  the  two  waves 
propagate  contradirectionally  through  the  photorefrac- 
tive  medium.  As  a  result,  the  latter  case  is  quite  dif¬ 
ferent  from  the  former  and  is  much  more  complicated  to 
analyze.  In  this  Letter  we  present  a  theoretical  analy¬ 
sis  and  experimental  investigation  for  the  case  of  re¬ 
flection  grating  interaction  in  the  nondepleted  pump 
regime.  The  absorption  effect  is  also  included. 

In  the  case  of  contradirectional  two-wave  mixing  in  a 
purely  diffusive  photorefractive  medium  the  coupled- 
wave  equations  for  the  slowly  var3dng  amplitudes 
El (2,  t)  and  E2(z,  t)  can  be  written  as^ 


aEi  ^ 

r  Q  •  E2  a  „ 

—  «  ...  ■“1  , 

(1) 

dz  V  dt 

2  Ii  +  I2  2 

aE2  1  ^£2 

y  Q*  ■  El  a 

(2) 

dz  V  dt 

2  Ii  +  I2  2 

where  y  is  the  intensity  coupling  constEuit,  a  is  the  in¬ 
tensity  absorption  coefficient,  i;  is  the  group  velocity, 
and  El  exp(— +  ikz)  and  E2  expi-iiot  —  ikz)  are 
the  coupled  quasi-monochromatic  waves  that  interact 
through  a  dynamic  photorefractive  grating  8n(Zyt)  * 
Q  exp(2iA2)  +  c.c.  The  dynamics  of  the  photorefrac- 
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tive  grating  can  be  written  as 

t^  +  Q  =  EiE2*,  (3) 

at 

where  r  is  the  relaxation  time  constant. 

The  temporal  behavior  of  each  wave’s  complex  ampli¬ 
tude  can  be  modeled  with  a  stationary  random  process, 
with  coherence  time  8(0^^  being  substantially  less 
than  the  relaxation  time  of  the  material,  i.e.,  8(jjt  » 
1.^  When  the  optical  path  difference  of  the  two  waves 
is  smaller  than  the  coherence  length  of  the  source  laser 
wave,  a  dynamic  photorefractive  grating  is  recorded  in 
the  mediiun.  Its  position  and  profile  are  neau"ly  tem¬ 
porally  constant.  Then,  as  an  approximation,  we  can 
replace  the  dynamic  grating  amplitude  Q  in  Eqs.  (1) 
and  (2)  with  its  ensemble  average  Q  (Q)  =  (E1E2  ). 

Since  the  complex  amplitudes  'E\(z^t)  and  E2(2, 
are  stationary  random  processes,  we  can  define 
some  of  their  ensemble  averages  as  ri2U,A^)  = 
<Ei(2,^i)E2*(2,^2)>,  rii(2,A0  =  <Ei(2,^i)Ei*(2,^2)>, 
and  r22(2,A0  =  <E2(2,^i)E2*(2,^2)>,  where  M  = 
^1  “  ^2-  We  refer  to  Fn  and  r22  as  the  self-coherence 
of  the  signal  wave  Ei(z,  t)  and  the  pump  wave  £2(2,  t)y 
respectively,  and  ri2  as  the  mutual  coherence  between 
the  signal  wave  and  the  pump  wave.  With  this 
notation  we  can  immediately  write  Q  =  ri2(2,0), 
h  =  rii(2,0),  and  h  =  r22(-2,0).  Equations  (l)-(3) 
can  therefore  be  reduced  to  a  system  for  these  average 
values®: 

ari2(2,A^)  ^  _  2  OTniZyM)  ^  X 

dz  V  dAt  2  h  +  I2 

x[ruU,A0  +  r22(2,A0],  (4) 

dru(zyAt)  y  ri2(2,o)  ♦  .  V  ri2*u,o) 

X  ri2(z,A<)  ~  arii(z,A<),  (5) 
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<^^22(2,  AQ 
dz 


y_ 

2 


ri2(^.0) 

/l  +/2 


ri2*u.-Af)  +  ^ 


ri2*(z.o) 

/l  +  /2  . 


X  ri2(^»Ao  +  ar22izt^t) ,  (6) 

Note  that  Eqs.  (4)-(6)  are  a  set  of  self-consistent 
partial  differential  equations  governing  the  propaga¬ 
tion  of  the  mutual-coherence  and  self-coherence  func¬ 
tions  of  the  two  waves.  We  can  solve  them  numerically 
as  an  initial-value  problem,  if  we  have  the  complete 
boundary  conditions  at  either  one  of  the  two  bound¬ 
aries  that  the  two  waves  are  incident  upon.  In  gen¬ 
eral,  we  know  only  the  self-coherence  functions  and  the 
mutual-coherence  function  of  the  two  waves  before  they 
enter  the  medium.  In  the  case  of  pump  depletion,  ei¬ 
ther  wave  has  changed  significantly  as  it  reaches  the 
second  boundary.  Therefore  the  complete  boundary 
conditions  are  unavailable  at  either  boundary.  In  the 
case  of  the  nondepleted  pump,  we  can  assume  that  the 
pump  wave  passes  through  the  photorefractive  medium 
unaffected  by  the  weak  signal  wave.  In  this  case  we 
obtain  the  complete  boundary  conditions  at  the  bound¬ 
ary  where  the  signal  wave  enters  the  medium.  We 
then  can. solve  Eqs.  (4)-(6)  for  the  self-coherence  func¬ 
tion  of  the  signal  wave  and  the  mutual-coherence  func¬ 
tion  of  the  two  waves. 

To  determine  the  boundary  conditions  in  the  nonde¬ 
pleted  pump  regime,  we  assume  that  both  the  signal 
wave  and  the  pump  wave  are  derived  from  the  same 
source  wave,  as  is  shown  in  Fig.  1.  If  the  source  wave 
has  a  Gaussian  line  shape  with  a  linewidth  of  Au,  then 
the  normalized  self-coherence  function  of  the  source 
wave  can  be  written  as 


rs(5^)  =  exp 


At  the  signal-wave  incident  plane  z  —  0  the  intensity 
ratio  between  the  signal  wave  and  the  pump  wave  is 
assumed  to  be  In  our  simulation  the  intensity  of 
the  pump  wave  is  taken  to  be  1.  Then  the  boundary 
conditions  at  z  =  0  can  be  written  as 


ri2{z  =  +  St),  (8) 

ru(2  =  0,Af)  =  /3r,(Af).  -  (9) 


coupling  between  the  two  waves  decreases,  and  the  sig¬ 
nal  intensity  approaches  a  certain  limit  that  is  much 
less  than  the  pump  intensity.  For  a  large  coupling 
constant  the  mutual-coherence  function  increases  as  z 
increases  because  of  relatively  strong  coupling.  Thus, 
'as  z  increases,  the  coupling  between  the  two  waves 
increases,  and  the  signal  intensity  increases  exponen¬ 
tially.  In  other  words,  for  a  large  coupling  constant 
the  interaction  length  of  the  two  waves  can  be  much 
longer  than  the  coherent  length  of  the  source  wave. 

The  above  theory  is  validated  experimentally.  Re¬ 
ferring  to  Fig.  1,  we  consider  two  partially  coherent 
waves,  obtained  by  splitting  an  argon  laser  wave  with 
a  linewidth  of  1.83  GHz.  The  signal  and  pump  waves 
are  contradirectionally  incident  upon  a  KNbOa^Co 
crystal  (y  =  3.3  cm''^  a  =  0.5  cm"\  and  thickness  d  = 
0.72  cm).  The  optical  path  difference  of  the  two  waves 
at  the  signal-wave  incident  plane  z  =  0  was  set  to 
be  AL  =  L2  —  Li,  To  monitor  the  mutual  coherence 
between  the  signal  wave  and  the  pump  wave  at  the 
output  plane  z  =  d,  we  employ  another  reference 


Fig.  1.  Schematic  of  the  two-wave-mbdng  configuration 
used  in  our  calculations  and  experiments.  The  distances 
L\  and  L2  are  the  optical  path  lengths  of  the  signal  wave 
and  the  pump  wave  from  the  laser  source  to  the  signal-wave 
incident  plane  z  =  0,  respectively.  L2ref  is  the  optical  path 
length  of  reference  wave  from  the  laser  source  to  the 
signal  output  plane  z  —  d. 


^22(^-o,A^)^  r,(Ao,  (10) 

where  the  time  delay  between  the  two  waves. 

In  Fig.  2  we  show  the  mutual-coherence  function 
ri2(z,0)  of  the  two  waves  and  the  self-coherence  func¬ 
tion  rii(z,0)  of  the  signal  wave  as  a  function  of  the 
position  z  inside  the  photorefractive  medium  for  two 
coupling  constants  y  of  3  and  7.  At  the  signal  entrance 
plane  z  =  0  the  time  delay  5^  is  0  s,  and  the  input  in¬ 
tensity  ratio  ^  is  10“"^.  The  linewidth  Ao  of  the  source 
is  10  GHz,  and  the  index  refi*action  of  the  photorefrac¬ 
tive  medium  is  2.3.  As  a  result  of  the  coupling,  part  of 
the  pump  wave  branches  off  in  the  direction  of  the  sig¬ 
nal  wave,  retaining  its  temporal  profile.  Therefore  the 
mutual-coherence  function  will  increase  with  coupling 
constant  y .  For  a  small  coupling  constant  the  mutual- 
coherence  function  eventually  decreases  to  zero  as  z  in¬ 
creases  because  of  the  relatively  rapid  increase  of  the 
time  delay  of  the  two  waves.  Thus,  as  z  increases,  the 


Fig.  2.  Mutual-coherence  function  ri2(z,0)  and  the  self- 
coherence  function  of  the  signal  wave  rii(z,0)  as  a  function 
of  z  for  coupling  constants  y  =*  3  (solid  curves)  and  y  —  1 
(dashed  curves). 
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(a) 

Fig.  3.  Interference  patterns  of  the  signal  wave  and  the 
reference  wave  at  the  output  plane  PI  (a)  without  and 
(b)  with  coupling.  Note  the  increase  of  fringe  visibility 
owing  to  the  coupling. 


1.2 
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Fig.  4.  Mutual  coherence  ri2(ci,0)  as  a  function  of  the 
optical  path  difference  AL. 

wave  £2ref  that  was  split  from  the  pump  wave  £2- 
The  optical  path  difference  of  waves  .Bi  and  B2ref 
was  adjusted  to  be  the  same  as  that  of  waves  Ex 
and  E2  at  the  output  plane  z  —  d.  Using  a  simple 
homodyne  technique,  we  observed  the  interference 
fringes  generated  by  waves  Ex  and  B2ref  with  a  CCD 
camera  at  the  output  plane  PI.  The  normalized 
mutual  coherence  ri2(d,  0)/[ Fii (d ,  0) F 22 {d ,  0)]^  can  be 
estimated  as  (/max  “  /min)/(4>/7i^),  where  (/max  “  /min) 
is  the  amplitude  of  the  fringes.  In  our  experiment  we 
monitored  the  interference  pattern  with  and  without 
pump  beam  B2.  Figure  3  shows  photographs  taken 
with  a  normalized  mutual  coherence  ri2(0,0)  »  0.43  at 
z  -  0  (AL  —  4  cm)  and  an  intensity  ratio  P  =  0.00151. 
The  measured  normalized  mutual  coherence  increases 
from  0.19  to  0.7  at  z  =  d.  We  then  measured  the 
normalized  mutual  coherence  as  a  function  of  the 
optical  path  difference.  Figure  4  shows  the  data 
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Fig.  5.  Intensity  gain  of  the  signal  wave  at  2  =  d  plane  as 
a  function  of  the  path  difference  AL. 

obtaii.ed  from  this  measurement  (filled  circles)  and 
the  theoretical  curve  for  the  same  parameters.  Note 
that  the  jiormalized  mutual  coherence  is  close  to 
1  for  sm^l  optical  path  differences  and  decreases 
quickly  to  zero  as  the  optical  path  difference .  in¬ 
creases  beyond  certain  values.  The  intensity  gain 
of  the  signal  wave  was  also  measured.  Figure  5 
shows  the  measurement  of  the  intensity  gain  (filled 
circles)  of  the  signal  wave  at  the  z  —  d  plane  as 
a  function  of  the  optical  path  difference  AL,  along 
with  the  theoretical  curve  for  the  same  parameters. 
Excellent  agreement  between  theory  and  experi¬ 
ment  has  been  achieved  in  the  case  of  the  signal-wave 
intensity  gain. 

In  conclusion,  we  have  studied  contradirectional  two- 
wave  mixing  with  partially  coherent  waves  in  photore- 
fractive  crystals  both  theoretically  and  experimentally. 
A  set  of  partial  differential  equations  has  been  derived 
to  describe  the  propagation  of  the  mutual-coherence 
and  self-coherence  functions  of  the  two  waves  in  the 
nondepleted  pump  regime.  Excellent  agreement  has 
been  achieved  between  theory  and  experiment. 

This  research  was  supported  by  the  U.S.  Office 
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We  investigate  contradirectional  two-wave  mixing  vidth  partially  coherent  waves  in  photorefractive  crystals. 
By  use  of  a  statistical  theory  on  linear  systems,  a  pneral  formulation  of  the  problem  in  the  space  and  fre¬ 
quency  domain  is  derived  and  implemented  numerically.  We  obtain  results  on  beam  intensity  and  mutual 
coherence.  The  results  on  the  enhancement  of  mutual  coherence  are  compared  with  previous  theoretical  re¬ 
sults  on  simpler  cases  and  with  experimental  measurements.  Excellent  agreements  are  achieved.  The  re¬ 
sults  also  indicate  that  the  effective  interaction  length  can  be  significantly  longer  than  the  coherence  length  of 
the  waves.  ©  1997  Optical  Society  of  America  [30740-3224(97)01906-1] 


1.  INTRODUCTION 

Two-wave  mixing  in  photorefractive  crystals  has  been  in¬ 
vestigated  extensively  for  many  applications  including 
image  amplification,  laser-beam  cleanup,  spatial  light 
modulators,  thresholding,  and  power-limiting  devices.^’^ 
Most  of  the  theoretical  works  in  this  area  are  based  on 
wave-mixing  with  mutually  coherent  waves.  However, 
in  some  applications,  such  as  self-pumped  and  mutually 
pumped  phase-conjugate  mirrors^"®  and  photorefractive 
filters,^  the  effect  of  partial  temporal  coherence  in  a  two- 
wave-mixing  process  cannot  be  ignored.  Two-wave  mix¬ 
ing  with  partially  coherent  weaves  has  been  studied  previ¬ 
ously  for  the  case  of  transmission-grating  interaction.® 
In  the  case  of  transmission-grating  interaction  the  optical 
path  difference  between  the  two  interacting  waves  re¬ 
mains  approximately  the  same  as  the  two  waves  propa¬ 
gating  codirectionally  through  the  photorefractive  me¬ 
dium,  especially  when  the  incident  angles  of  the  two 
waves  are  close  to  each  other.®  In  the  case  of  reflection¬ 
grating  interaction  the  optical  path  difference  between 
the  two  interacting  waves  varies  significantly  as  the  two 
waves  propagate  contradirectionally  through  the  photore¬ 
fractive  medium.  Thus  for  the  case  of  transmission¬ 
grating  interaction,  only  one  free  variable  for  the  position 
is  needed  to  describe  in  a  self-consistent  way  the  second- 
order  statistical  properties  of  the  two  optical  waves,  i.e., 
their  intensities  and  mutual  coherence,  while  at  least  two 
free  variables,  one  for  the  position  and  one  for  the  optical 
path  difference,  will  be  needed  for  the  case  of  reflection- 
grating  interaction.  Another  difficulty  in  studying  the 
reflection-grating  interaction  of  partially  coherent  waves 
is  to  find  a  way  to  incorporate  the  complete  boundary  con¬ 
ditions  into  the  theoretical  formulation  as  a  result  of  the 
two-point  bound  ary- value  problem.  In  such  a  problem  a 
complete  set  of  boundary  conditions  includes  both  the 
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second-order  self-statistical  properties  (e.g.,  self¬ 
coherence)  of  each  wave  at  its  entrance  boundary  and  the 
second-order  mutual  statistical  properties  (e.g.,  mutual 
coherence)  of  the  two  waves  at  their  respective  entrance 
boundaries.  In  a  recent  work  we  provided  a  theoretical 
formulation  of  the  problem  in  the  space  and  time  domain 
for  the  reflection-grating  interaction  in  the  nondepleted- 
pump  regime.®  By  using  the  nondepleted-pump  approxi¬ 
mation,  we  reduced  the  two-point  boundary-value  prob¬ 
lem  to  an  initial  value  problem.  In  this  paper  we  present 
a  general  formulation  of  the  problem  in  the  space  and  fre¬ 
quency  domain  based  on  the  standard  statistical  theory 
on  linear  systems.  The  general  formulation  is  also 
implemented  numerically.  Specifically,  we  investigate 
the  signal-intensity  gain  and  the  mutual  coherence  in  the 
contradirectional  wave  mixing  of  two  partially  coherent 
waves.  Contrary  to  conventional  belief,  we  discover  that 
the  effective  interaction  length  (or  grating  length)  can  be 
significantly  longer  than  the  coherence  length  of  the  inci¬ 
dent  waves.  The  results  are  also  compared  ^^ath  previous 
theoretical  results  on  simpler  cases  and  with  experimen¬ 
tal  measurements. 


2.  THEORETICAL  MODEL 

Photorefractive  two-wave  mixing  is  a  nonlinear  optical 
process.  Because  of  the  mutual  coherence  of  the  two 
waves,  a  dynamic  holographic  grating  is  formed  in  the 
medium.  Its  position  and  index  profile  are  nearly  sta¬ 
tionary  under  the  condition  of  a  cw  illumination.  Both 
waves  are  scattered  into  each  other  by  the  presence  of 
this  index  grating.  Scattering  of  partially  coherent 
waves  by  a  stationary  grating  can  be  modeled  with  a  sta¬ 
tistical  theory  on  linear  systems.^®  An  iterative  proce¬ 
dure  can  subsequently  be  devised  to  obtain  the  final  pho- 
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torefractive  grating  profile  from  an  initially  arbitrary 
grating  profile. 

As  is  shown  in  Fig.  1,  two  counter  propagating  waves 
with  partial  coherence  enter  a  photorefractive  medium  at 
r  =  0  and  z  -  Ly  respectively.  The  electric  field  in  the 
photorefractive  medium  can  be  written  as 

E(2,  t)  =  Ei(2,  ^)exp(“icoo^  +  I^qz) 

+  E2(2,  t)exj>i-i(oot  -  iko2),  (1) 

where  cjq  is  the  center  frequency  of  the  two  partially  co¬ 
herent  waves,  ko  =  is  the  corresponding  wave  vec¬ 

tor,  and  n  is  the  refractive  index  of  the  photorefractive 
medium.  Reflecting  the  partial  coherence,  EiU,  t)  and 
£2! -2^,  t)  are  stationary  random  variables.  They  repre¬ 
sent  the  random  fluctuation  of  the  amplitudes  of  the  two 
waves.  For  the  convenience  of  our  later  discussion  we 
will  now  briefly  describe  some  notations  and  definitions 
for  the  second-order  statistical  properties  of  the  two  opti¬ 
cal  waves.  Let  r)  =  (EjCz,  ^i)Ei*(2,  ^2^)  and 

r22(-^»  =  (^2(2,  ^i)E2*(2,  ^2)}  denote  the  self¬ 

coherence  functions  of  Ei(z,  and  E^fe,  ti)y  respec¬ 
tively,  and  Ti2(Zy  r)  =  (Ejfz,  ^i)E2*(2,  ^2))  the 

mutual-coherence  function  between  Eilz,  and 
£2(2,  t2)y  where  r  =  -  ^2  is  the  time  delay  and  () 

means  ensemble  average.  Let  £11(2,  Aw)  and 
£22(2,  Aw)  denote  the  self-spectral-density  functions  of 
Ei(2,  t)  and  £2(2,  ^),  respectively,  and  £12(21,  22,  Aw) 
be  the  cross-spectral-density  function  between  £1(2 1,  t) 
and  £2(22,  t).  The  spectral-density  functions  and  the 
corresponding  coherence  functions  are  Fourier-transform 
pairs,  i.e., 

Fii(2,  r)  =  f  £11(2,  Aw)exp(“iAwr)dAw,  (2) 


^22(^1  r)  =  J  £22(2,  Aw)exp(-iAwr)dAw, 

(3) 

^  Aw)exp(-fAwT)dAw, 

(4) 


Fig.  1.  Two-wave  mixing  in  photorefractive  crystals  modeled  as 
a  linear  system  with  the  signal-wave  entrance  plane  and  the 
pump-wave  entrance  plane  as  two  input  planes  and  any  arbi- 
frary  plane  inbetween  as  the  output  plane. 


where  Aw  =  w  -  wo»  with  w  being  a  general  frequency 
component  of  the  waves.  With  the  above  notations  and 
relations  the  intensity  of  the  two  waves  can  be  expressed 
as 


Il(z) 

- 

o,=  J 

1*  Eii(z. 

Aa;)dAa>, 

(5) 

12(2 ) 

=  r22(2, 

o'-J 

f  £22(2, 

Aa>)dAaj, 

(6» 

and  the  mutual  coherence  of  the  two  waves  can  be  ex¬ 
pressed  as 

ri2(2,  0)  =  J  Ei2(2,  Z,  Saj)dAw.  (7) 

With  these  equations  the  intensity  of  the  waves  as  well  as 
the  mutual  coherence  can  be  obtained  as  soon  as  the 
spectral-density  functions  are  obtained.  In  what  follows 
we  will  derive  the  spectral-density  functions  by  using  the 
statistical  approach. 

Using  the  above  notations  and  definitions,  we  now  be¬ 
gin  our  discussion  on  the  photorefractive  interaction.  In¬ 
side  the  photorefractive  medium,  a  dynamic  index  grating 
is  generated.  It  can  be  written  as 

y  c  [Q(2,  i) 

Sn  =  “f  jr —  .  exp(2tko2)  +  c.c.  ,  (8) 

where  Q(2,  t)  is  a  measure  of  the  index  grating,  y  is  the 
intensity  coupling  coefficient,  and/o(2)  =  /i(2)  +  72(2)  is 
the  total  intensity  at  position  2.  For  the  purpose  of  our 
discussion  we  will  call  Ej  the  signal  wave  and  £2  the 
pump  wave.  Thus  for  a  photorefractive  grating  with  a 
positive  7,  the  signal  wave  E^  can  be  amplified.  If  the 
photorefractive  effect  is  based  purely  on  carrier  diffusion 
(e.g.,  BaTiOs,)  the  dynamics  of  the  index  grating  is  de¬ 
scribed  by  the  following  equation: 

dQ(2,  t) 

^ph  - -  +  Q(2,  0  =  £1(2:,  t)E2*(Zy  t)y  (9) 

where  is  the  relaxation-time  constant.  By  virtue  of 
photoexcitations,  photorefractive  processes  are  usually 
slow  at  low  intensities.  It  is  reasonable  to  assume  that 
the  coherence  time  of  the  two  partially  coherent 

waves  is  much  smaller  than  the  relaxation  time  Xp^  of  the 
photorefractive  medium,  i.e.,  Scot^^  >  1.  Since  £1(2,  t) 
and  £2(2,  t)  are  stationary  random  variables,  we  can 
make  the  following  approximation  (see  Appendix  A)®: 

Q(2,  t)  =  (Q(2,  t))  =  ri2(2,  0).  (10) 

In  other  words,  two  partially  coherent  waves  with  their 
complex  amplitudes  fluctuating  randomly  with  time  can 
actually  write  a  stationary  grating  in  a  photorefractive 
medium  under  the  appropriate  conditions.  For  simplic¬ 
ity  we  will  denote  Q(2,  t)  with  Q(2)  from  now  on.  Note 
that  Q(2)  is  also  the  mutual  coherence  of  the  two  waves 
at  position  2.  We  note  that  the  approximation  is  also 
valid  when  ^wXph  1. 

Given  arbitrary  functions  of  Q(2)  and  7o(2),  Eq.  (8) 
yields  an  index  grating.  The  propagation  of  a  monochro¬ 
matic  wave  through  such  an  index  grating  can  be  de¬ 
scribed  by  the  coupled  wave  equations.  The  coupled 
monochromatic  waves  can  be  written  as 
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E  =  Ei(2,  a>)exp(-’fco^  +  iVz) 

+  E2(2,  a>)exp(-icot  -  /kz),  (11) 

where  Ej  and  E2  are  the  amplitudes  of  the  monochro¬ 

matic  components,  w  is  the  optical  wave  frequency,  and 
k  =  n  w/c  is  the  optical  wave  vector.  The  coupled  wave 
equations  can  be  written  as 

(?Ei(z,  w)  y  - 

a,)exp(-2fAk2) 

-^£1(2,0.),  (12) 

^Eo(z,  cti)  y  — 

^  Q*U)Ei(r.  c.)exp(2iAk2) 

+  Je2(z,  a>),  (13) 

where  Ak  =  k  -  ko  is  the  phase  mismatch  between  the 
optical  waves  and  the  index  grating  and  a  is  the  intensity 
absorption  coefficient.  With  sufficient  boundary  condi¬ 
tions,  Eqs.  (12)  and  (13)  can  be  solved  either  analytically 
in  some  special  cases  or  numerically  in  general.  When 
the  boundaiy  conditions  are  Ei(z  =  0,  a>)  =  1  and  £2(2 
=  L,  w)  =  0,  the  solutions  (output)  are  denoted  as 
Ei(z,  (o)  =  0))  and  £2(2,  (o)  =  Hi2(2,  a>).  When 

the  boundary  conditions  are  £1(2  ==  0,  a;)  =  0  and  £2(2 
=  L,  to)  =  1,  the  solutions  (output)  are  denoted  as 
£1(2,  w)  =  H2i(2,  (jj)  and  £2(2,  a>)  =  H22(2,  oj).  As  a 
linear  system,  the  general  solutions  are  linear  combina¬ 
tions  of  Hu  ,  Hi2  ,  H21 ,  and  H22 . 

In  each  of  the  iterations  an  arbitrary  stationary  index 
grating  can  be  considered  as  a  linear  system.  Referring 
to  Fig.  1,  we  consider  a  given  stationary  index  grating  in 
the  photorefractive  medium  as  a  linear  system  with  the 
optical  waves  at  the  boundary  planes  2  =  0  and  z  =  L  as 
the  input  and  the  optical  waves  at  an  arbitrary  plane  z  as 
the  output.  According  to  the  theory  on  the  statistical 
properties  of  linear  systems,  the  second-order  statistical 
properties  of  the  optical  waves  at  the  output  plane  can  be 
expressed  in  terms  of  the  second-order  statistical  proper¬ 
ties  of  the  optical  waves  at  the  input  planes  and  the  fre¬ 
quency  response  of  the  linear  system.  To  be  specific,  the 
input  at  the  2  =  0  plane  is  £1(2  =  0,  i),  the  input  at 
2  =  L  plane  is  £2(2  =  L,  t),  and  the  outputs  at  the  2 
plane  are  £1(2,  t)  and  £2(2,  t).  The  frequency  response 
of  this  linear  system  can  be  expressed  in  terms  of  the  so¬ 
lutions  of  Eqs.  (12)  and  (13).  With  the  notations  de¬ 
scribed  in  the  previous  paragraph,  the  frequency  re¬ 
sponses  from  the  input  £1(2  =  0,  t)  to  the  outputs 
£1(2,  t)  and  £2(2,  t)  are 

Hu'(2,  a>)  =  Hii(z,  aj)exp(f  Akz),  (14) 

lii2(2,  0))  =  Hi2(e»  a>)exp(~iAk2), 

(15) 

respectively.  Similarly,  the  frequency  responses  from 
the  input  £2(2  =  L,  O  to  the  outputs  £1(2,  t)  and 
£2(2,  t)  are 

H2i'(2,  cj)  =  H2i(2,  o))exp(iAkL)exp(iAk2),  (16) 


H22'U,  0))  =  H22(2,  aj)exp(iAkL)exp(-iAk2),  (17) 

respectively.  The  additional  phase  terms  in  Eqs.  (14)- 
(17)  account  for  the  difference  Ak  =  k  -  kp  owing  to  a  fi. 
nite  Aoj  =  cuq.  With  these  spectral-response  func¬ 
tions  we  can  express  the  spectral-density  functions  of  the 
two  outputs  £1(2,  t)  and  £2(2,  t)  in  terms  of  the  spectral- 
density  functions  of  the  two  inputs  £1(2  =  0,  0  and 
£2(2  =  L,  t).  Thus,  according  to  their  definitions  and 
Eqs.  (14)-(17),  we  obtain 

Eii(2,  Aa>)  =  Hu(2,  (o)Hu*(2,  a>)£u(2  =  0,  Aa>) 

+  H2i(2,  w)H21*(2,  a>)£22{2  =  L,  AiJt>) 

+  Hn(2,  a;)H2i*(2,  o)) 

X  £12(2  =  0,  2  =  L,  Aa>)exp(-“iAkL) 

+  H2i(2,  £j)Hu*(2,  CO) 

X  [£22(2  =  0,  2  =  Lff  Actj) 

X  exp(“iAkL)]*,  (18) 

£22(2,  Acti)  =  Hi2(2,  w)Hi2*{z,  <*j)£ii(2  =  0,  Aa>) 

+  H22(2,  w)H22*(2,  (d)E22(z  =  L,  Aw) 

-1“  Hi2(2,  w)H22*(2,  w) 

X  £12(2  =  0,  2  =  L,  Aw)exp(-iAkL) 

+  H22(2,  w)Hi2*(2,  w) 

^  [®12(^  =  0,  2  =  L,  Aw) 

X  exp(-f AkL)]*,  (19) 

£12(2,  Zy  Aw)  =  exp(f2Ak2){Hii(2,  w)Hi2*(2,  w) 

X  £11(2  =  0,  Aw) 

. H2i(2,  w)H22*(2,  w) 

X  £22(2^  ” 

+  Hu(2,  w)H22*(2^,  <o) 

X  £^2(2  =  0,  2  =  Ly  Aw) 

X  exp(-fAkL)  +  H2i(2,  w)Hi2*(2,  <d) 

X  [£^2(2  =  0,  2  =  L,  Aw) 

X  exp(~fAkL)]*}.  (20) 

Note  that  the  spectral-density  functions  £11(2  =  0,  Aw), 
£22(2  =  Ly  Aw),  and  £12(^1  =  0,  22  =  Ly  Aw)  of  the  two 
inputs  £1(2  =  0,  t)  and  £2(2  =  L,  t)  are  given  as  the 
boundary  conditions. 

Two  physical  processes  happen  simultaneously  during 
two-wave  mixing  in  a  photorefractive  medium.  First,  the 
two  optical  waves  propagate  through  the  photorefractive 
medium  while  being  scattered  into  each  other  by  the  in¬ 
dex  grating.  Second,  the  scattered  waves  modify  the  in¬ 
dex  grating  through  the  photorefractive  effect  until  a 
steady  state  is  reached.  We  have  provided  above  a  math¬ 
ematical  model  that  describes  these  two  physical  pro¬ 
cesses  separately.  A  steady  state  of  the  two-wave  mixing 
in  the  photorefractive  medium  is  reached  when  the  two 
optical  waves  scattered  by  the  photorefractive  grating  can 
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exactly  sustain  the  same  photorefractive  grating.  A 
steady-state  solution  of  the  two-wave  mixing  in  a  photo- 
refractive  medium  can  thus  be  obtained  by  use  of  the 
mathematical  model  described  above  through  an  iterative 
procedure.  The  procedure  is  outlined  as  follows: 

Step  1:  Give  an  initial  guess  on  the  function 
Q(z)/Iq(z), 

Step  2:  Solve  Eqs.  (12)  and  (13)  for  the  functions 
t*^)  (i,y  -  1»  2)  with  the  function  Q(z)/1q(z)  pro¬ 
vided  in  the  last  step. 

Step  3:  Obtain  the  spectral-density  functions 
Eii(Zy  Aa>),  £22(2,  ^w),  and  £12(2,  Zy  Aw)  by  use  of  Eqs. 
(18)-(20). 

Step  4:  Obtain  Ii(z),  72(2),  Q(a:),  and  Q{z)/Iq(z)  by 
use  of  Eqs.  (5),  (6),  (7),  and  (10). 

Step  5:  Compare  the  new  version  of  the  index  grating 
Q(z)/Iq{z)  and  the  previous  version.  If  they  are  close 
within  a  certain  accuracy  requirement,  the  solution  has 
been  obtained.  Otherwise,  the  iteration  continues  by  use 
of  the  new  version  of  the  index  grating. 

To  determine  the  boundaiy  conditions,  we  assume  that 
both  the  input  optical  waves  are  derived  from  the  same 
laser  source.  If  the  source  laser  wave  has  a  Gaussian 
line  shape  with  a  FWHM  linewidth  of  Sw,  then  the  nor¬ 
malized  spectral-density  function  of  the  source  laser  wave 
can  be  written  as 


Taking  P  as  the  incident-intensity  ratio  /^(z  =  0)/72(2 
=  L)  of  the  two  optical  waves  at  their  respective  en¬ 
trance  boundary  planes,  we  obtain  the  boundary  condi¬ 
tions  as 


£11(2  -  0,  Aw)  =  )3£s3(Aw),  (22) 


£22(^  —  0,  Aw)  “  £gg(Aw), 


(23) 


Ei2(2i  =  0,  Z2  =  L,  Aw)  =  v)3£g3(Aw)exp(“ikoL) 


X  exp(-”iw^rf),  (24) 

where  is  the  time  delay  between  the  optical  waves 
when  they  reach  their  respective  entrance  planes.  In  de¬ 
riving  the  above  boundary  conditions,  we  have  assumed 
that  the  laser  source  has  a  Gaussian  line  shape.  There  is 
no  loss  of  generality  in  this  assumption.  Similar  results 
can  be  obtained  with  a  different  line  shape. 

This  completes  the  general  formulation  to  model  con- 
tradirectional  two-wave  mixing  in  photorefractive  crys¬ 
tals. 


3.  DISCUSSIONS  AND  SIMULATIONS  ^ 

To  clarify  the  complicated  formulation  described  above, 
we  consider  some  simple  cases  first  before  presenting  the 
numerical  and  experimental  results. 

In  the  absence  of  coupling  ( y  =  0)  in  a  lossless  medium 
(a  =  0)  with  p  -  1  the  above  formulation  describes  the 
interference  of  two  coimterpropagating  optical  waves  in  a 


dielectric  medium.  In  this  case  it  takes  only  one  itera¬ 
tion  to  obtain  the  solution,  which  is 

£l2(^,  z.  Aw)  =  £g8(A4ti)exp(-ikL) 

X  exp(“tw^rf)exp(2iAkz).  (25) 

Fourier  transforming  the  above  equation  over  w  and  tak¬ 
ing  into  account  the  extra  phase  term  exp(2fko2)  resulted 
from  the  definition  of  £i(z,  t)  and  £2(2,  ^)  in  Eq.  (1),  we 
obtain  the  mutual  coherence  of  the  two  waves  as  a  func¬ 
tion  of  position  z, 

/  nL  -  2nz\ 

==  rj  r  +  +  - - -  exp(-iwo^d)» 

•  (26) 

where  r5(r)  is  the  normalized  self-coherence  function  of 
the  source  laser  wave  and  n  is  the  index  of  refraction  of 
the  medium.  The  result  has  been  well  established  in  the 
literature.  This  example  can  also  help  us  see  more 
clearly  the  subtle  contribution  of  the  phase  terms  in  the 
above  formulation.  In  our  definition,  =  0  ifz  =  L/2  is 
the  plane  of  zero  path  difference. 

When  the  coherence  time  <5w’^  is  much  larger  than  the 
time  delay  involved  in  the  above  formulation,  i.e.,  Swt^ 
<  1  and  S(x)nL/c  <  1,  the  normalized  spectral-density 
function  of  the  source  laser  wave  can  be  written  as 
From  the  above  formulation  we 
can  obtain  the  set  of  equations  governing  the  intensities 
of  the  two  optical  waves, 


d 

IJ2 

(27) 

d 

hh 

dJ^‘ 

(28) 

These  two  equations  are  exactly  the  same  as  those  ob¬ 
tained  for  two-wave  mixing  of  monochromatic  waves. 
We  note  that  So)t^  1  is  a  condition  for  partially  coher¬ 
ent  waves  to  be  treated  as  monochromatic  waves  in  con- 
tradirectional  two-wave  mixing  in  a  photorefractive  me¬ 
dium.  We  also  recall  that  the  upper  limit  on  the 
coherence  time,  i.e.,  ^wTph  >  1,  still  needs  to  be  satisfied 
to  reach  Eqs.  (27)  and  (28)  from  the  above  general  formu¬ 
lation.  This  is  the  main  difference  between  the  results  of 
Eqs.  (27)  and  (28)  in  this  paper  and  those  obtained  di¬ 
rectly  for  monochromatic  waves.  Note  that  we  did  not 
use  the  iterative  procedure  in  obtaining  Eqs.  (27)  and 
(28).  The  complete  boundary  conditions  for  Eqs.  (27)  and 
(28)  are  simply  7i(z  =  0)  and  72(z  =  L). 

We  now  use  the  above  general  formulation  to  obtain 
the  solution  in  the  nondepleted-pump  regime  (for 
y  ^  0)y  which  has  been  obtained  previously  by  a  different 
method.  Using  Eqs.  (12)  and  (13)  and  the  following  three 
definitions, 

Enizy  ^io)  =  (Ei(z,  a>)Ei*(z,  o>)),  (29) 

£22(2,  Acu)  =  (£2(2,  o>)E2*(z,  (o))y  (30) 

£12(2,  Zy  Aw)  ==  <Ei(z,  w-)E2*(z,  w))exp(2iAkz), 

(31) 

we  can  obtain  a  set  of  equations  governing  the  propaga¬ 
tion  of  the  spectral-density  functions  as 
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<9Eii(Zf  Aco) 
Bz 


6^22^^*  Acu) 

7z 


f?Ei2(2,  2,  Aa>) 
Bz 


+  Q*(2)Ei2(z,  Z,  Ao))] 

-  aEii(2,  Ao),  (32) 

^[Q(2)E:2*(2,2,  Ao,) 

+  Q*(2)Ei2(2,  Aco)] 

+  aE22(2,  Aco),  (33) 

2iAkEi2(2,  Aco)  +  2l[(z) 

X  [EiiiZy  Aco)  +  E22(2,  Aco)] 


--  aEi2{z,  Zy  Aco).  (34) 

Ei(z,  co)  and  £2(2,  co)  in  Eqs.  (29)-(31)  are  related  to  the 
Fourier-transform  coefficients  of  the  two  optical  waves  as 
defined  in  Eq.  (11).  Strictly  speaking,  a  stationary  ran¬ 
dom  process  cannot  be  Fourier  transformed  over  the  time 
variable  t.  However,  we  can  truncate  it  into  a  finite  du¬ 
ration  T  and  Fourier  transform  the  truncated  process. 
After  that,  we  can  first  obtain  a  set  of  equations  similar  to 
Eqs.  (32)-(34)  for  the  truncated  process,  then  let  T  go  to 
infinity  to  obtain  Eqs.  (32)-(34).  This  is  a  standard  pro¬ 
cedure  in  statistical  optics.^®  Fourier  transforming  the 
above  relation  over  Acu,  we  obtain  a  set  of  equations  gov¬ 
erning  the  propagation  of  the  seif  and  mutual  coherence 
of  the  two  optical  waves  as 


<^ri2(z,  r) 
Bz 


r) 

Bz 


2n  BTi2^z,  t)  y  ri2(z,  0) 
"c  ~Bt  2  /i  +  /2 

X  [^11(2,  r)  +  r22(2,  r)], 

y  ri2(2,  0) 


2  /i  +  /2 

y  ri2*(2,  0) 


ri2*(2,  -t) 


C?r22(2,  t) 
Bz 


2  h  +  h 

-  cjrrii(z,  r), 
y  ri2(2,  0) 


ri2(2, 


2  /i  +  /2 


ri2*(2,  -r) 


(35) 


(36) 


+ 


r  ri2*(2, 0) 

2  /1  +  /2 


r  12(^1 


r) 


+  ar22(z,  r).  (37) 

When  there  is  no  absorption  in  the  photorefractive  me¬ 
dium  (i.e.,  a  =  0),  we  can  obtain,  from  Eqs.  (36)  and  (37), 


^[rii(z,  r)  -  r22(z,  r)]  =  0.  (38) 

Further,  letting  r  =  0  in  Eq,  (38),  we  obtain 

^  [h{z)  -  h^z)]  =  0.  (39) 

Equation  (39)  shows  that,  when  there  is  no  absorption  in 
the  photorefractive  medium,  the  intensity  difference  be¬ 


tween  the  signal  wave  and  the  pump  wave  is  a  constant  of 
integration  (conservation  of  power  flow). 

Note  that,  although  Eqs.  (35)-(37)  are  general  and  self- 
consistent,  they  contain  only  the  mutual  coherence  of  the 
two  optical  waves  at  the  same  locations.  They  are  not 
compatible  with  the  kind  of  boundary  conditions  as  given 
by  Eq.  (24),  which  gives  the  mutual  statistical  properties 
of  the  two  waves  at  two  different  points  in  space.  There¬ 
fore  they  are  useful  only  in  the  nondepleted-pump  regime, 
where  we  can  assume  that  the  pump  wave  passes  through 
the  photorefractive  medium  without  changing  its  statisti¬ 
cal  properties.  Under  this  assumption  we  can  obtain 
from  Eqs.  (22)-(24)  the  self-  and  mutual-coherence  func¬ 
tions  of  the  two  waves  at  the  signal-wave  entrance  bound¬ 
ary  plane  (z  =  0)  as 


Tuiz  =  0,  r)  s  ^/^^s8(T  +  St)exp{-imoSt),  (40) 
r„(z  =  0,  r)  *  (41) 


=  0,  r)  ^  r„(T),  (42) 

where  St  -  nL/c  is  the  time  delay  between  the  two 
optical  waves  at  the  signal-wave  entrance  plane  and 


rs8(^)  =  exp^ 


SOJT 

4(ln  2)H  I 


(43) 


is  the  coherence  function  of  the  source  laser  wave.  We 
can  use  Eqs.  (40)‘-(42)  as  the  boundary  conditions  and  in¬ 
tegrate  Eqs,  (35)“(37). 

Equations  (35)-(37)  are  simple  not  only  in  the  sense 
that  they  can  be  implemented  easily  numerically,  but  also 
that  they  can  be  used  to  obtain  approximate  analytical  so¬ 
lutions  to  some  special  cases  within  the  nondepleted  re¬ 
gime.  One  such  case  is  that  the  coherence  length  of  the 
source  laser  wave  is  much  longer  than  the  two-wave¬ 
mixing  interaction  length  and  the  coupling  constant  is 
large.  This  occurs  when  we  use  a  multimode  argon  laser 
and  a  KNbOaiCo  crystal  for  the  two-wave-mixing  experi¬ 
ment.  In  this  case  we  can  neglect  the  term  that  contains 
the  partial  derivative  on  rin  Eq.  (35)  and  reduce  the  set  of 
partial  differential  equations,  i.e.,  Eqs.  (35)-(37),  to  a  set 
of  ordinary  differential  equations  that  can  be  solved  ana¬ 
lytically  under  the  nondepleted-pump  approximation. 
Remember  that  Eq.  (35)  is  derived  directly  from  Eq.  (34). 
We  can  see  from  Eq.  (34)  that  the  approximation  of  ne¬ 
glecting  the  term  that  contains  the  partial  derivative  on  r 
in  Eq.  (35)  implies  that  the  wave-vector  difference  Ak  of 
the  different  frequency  components  of  the  two  partially 
coherent  waves  are  negligible  with  respect  to  the  thick¬ 
ness  of  the  photorefractive  medium  (i.e.,  AkL  1). 
Similar  approximation  has  been  made  previously  by  Sax- 
ena  et  al.  in  the  study  of  multiple-beam  interaction  by 
transmission  gratings  in  the  photorefractive  media. 
Therefore  under  the  nondepleted-pump  approximation 
and  the  approximation  of  AkL  <  1  the  signal-wave  inten¬ 
sity  gain  and  the  normalized  mutual  coherence  of  the  two 
waves  can  be  obtained  as 


7i(z)  _  Tniz,  0)  ri2(0,  0)ri2*(0,  0) 

/i(0)  “  rn(o,o)  ru(0,0)r22(0,0) 

X  [exp(yz)  -  l]exp(*-az)  +  exp(-Qrz),  (44) 
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[r'ii(2, 0)r22('2>  0)] 
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'i2*(o,o)  ru(0,0)r22(0,0) 

[1  -  exp(-y2)]  + 


ri2(0,0) 


ri2^(0, 0) 


exp(-yz) 


1/2* 


(45) 


In  most  photorefractive  crystals  with  exp[(y-  a)z]  >  1 
the  signal-intensity  gain  is  affected  primarily  by 
(y  -  a).  We  also  note  that  the  normalized  mutual  co¬ 
herence  of  the  two  waves  is  affected  only  by  the  photore¬ 
fractive  coupling  constant  y. 

As  a  comparison  of  the  general  formulation,  the  simpli¬ 
fied  formulation  and  the  approximate  anal5d:ical  solu- 
-tions,  we  show  in  Figs.  2  and  3  the  signal-intensity  gain 
and  the  normalized  mutual  coherence  of  the  two  waves  at 
the  signal- wave  exit  plane  {z  =  L)  as  a  function  of  the  op¬ 
tical  path  difference  between  the  two  optical  waves  at  the 
signal-wave  entrance  plane  (2  =  0).  The  parameters 
are  y  =  3.0  cm"\  a  =  0.0  cm~\  n  =  2.3,  L  =  0.72  cm, 
and  =  27r  X  1.8  GHz,  and  P  =  10"'*.  These  are  typi¬ 
cal  parameters  in  our  experiment  when  we  use  a  multi¬ 
mode  argon  laser  and  a  KNbOsiCo  crystal  to  implement 
the  two-wave-mixing  experiment.  In  this  case  the  coher¬ 
ence  length  of  the  source  laser  wave  is  much  longer  than 
the  thickness  of  the  photorefractive  medium.  The  gen¬ 
eral  formulation  and  the  simplified  formulation  produce 
the  same  results  within  the  numerical  accuracy  in  the 
nondepleted  regime.  The  results  of  these  two  formula¬ 
tions  are  represented  by  the  dashed  curves  in  Figs.  2  and 
3.  We  notice  that  the  approximate  analytical  solution  re¬ 
tains  the  major  characteristics  of  the  exact  solution  and 
provides  a  good  understanding  of  the  interaction.  The 
results  of  the  approximate  analytical  solution  are  repre¬ 
sented  by  the  solid  curves  in  the  figures.  Note  that  the 
curves  obtained  from  the  approximate  anai3d:ical  solution 
are  symmetric  about  2©  =  0  owing  to  the  approximation 
of  neglecting  the  partial  derivative  over  r  in  Eq.  (35), 
while  the  curves  obtained  from  the  general  formulation 
are  shifted  to  the  right  side  (see  Figs.  2  and  3).  Here, 
2o  is  defined  as  the  path  difference.  The  difference  be¬ 
tween  the  results  of  the  approximate  analytical  solution 
and  those  of  the  general  formulation  will  decrease  as  the 
ratio  between  the  coherence  length  of  two  waves  and  the 
thickness  of  the  photorefractive  medium  increases. 

Now  let  us  consider  the  case  of  the  depleted  pump. 
Numerical  simulation  with  the  general  formulation  is  the 
only  means  to  analyze  this  case.  We  first  consider  the 
case  in  which  the  coherence  length  of  the  source  laser 
wave  is  finite  but  much  longer  than  the  thickness  of  the 
photorefractive  medium.  Again,  we  use  the  following  set 
of  parameters:  y  =  3.0  cm”^,  a  =  0.0  cm“^,  n  =  2.3, 
L  =  0.72  cm,  and  =  277  x  1,8  GHz.  The  solid 
curves  in  Fig.  4  show  the  signal  intensity  /i(2),  the  pump 
intensity  72(2 ),  and  the  normalized  mutual  coherence 
y{z)  =  Ti2{z,  0)/[rn(2i  0)r22(2:,  0)]^  as  functions  of  the 
position  2  in  the  photorefractive  medium.  The  optical 
path  difference  between  the  signal  wave  and  the  pump 
wave  at  the  signal-wave  entrance  boimdary  (2  =  0)  is 
chosen  to  be  zero.  The  incidence  intensity  ratio  of  the 
two  waves  is  chosen  to  be  one.  We  also  show  in  Fig.  4 
with  dashed  curves  the  signal  intensity  Iiiz)  and  the 


pump  intensity  72(2)  for  two- wave  mixing  with  monochro¬ 
matic  waves  for  the  purpose  of  comparison.  W’e  notice 
that  the  results  of  two-wave  mixing  with  partially  coher¬ 
ent  waves  are  veiy  close  to  those  of  two-wave  mixing  with 


Fig  2.  Signal- wave  intensity  gain  as  a  function  of  the  optical - 
path  difference  at  the  signal-wave  entrance  plane  in  the 
nondepleted-pump  regime.  The  dashed  curve  is  the  numerical 
solution.  The  solid  curve  is  the  approximate  analytical  solution. 


Fig.  3.  Mutual  coherence  of  the  two  waves  at  the  signal- wave 
exit  plane  as  a  function  of  the  optical  path  difference  at  the 
signal- wave  entrance  plane  in  the  nondepleted-pump  regime. 
The  dashed  curve  is  the  numerical  solution.  The  solid  curve  is 
the  approximate  analytical  solution. 


2(cm) 


Fig.  4.  Signal  intensity,  pump  intensity,  and  the  normalized 
mutual  coherence  as  a  fimrtion  of  position  z  in  the  photorefrac¬ 
tive  medium  for  partially  coherent  waves  (solid  curves)  and 
monochromatic  waves  (dashed  curves). 
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monochromatic  waves  when  the  optical  path  difference  of 
the  two  waves  is  small  compared  with  the  coherence 
length  of  the  source  laser  wave.  The  effect  of  partisd  co¬ 
herence  of  the  two  waves  in  this  case  is  significant  only 
for  a  large  optical  path  difference  of  the  two  waves.  We 
show  the  signal-intensity  gain  in  Fig.  5(a)  and  the  mutual 
coherence  between  the  signal  wave  and  the  pump  wave  at 
the  pump-wave  entrance  plane  (z  =  L)  in  Fig.  5(b)  both 
as  functions  of  the  optical  path  difference  between  the  two 
waves  at  the  signal-wave  entrance  plane  (z  =  0)  for  vari¬ 
ous  intensity  ratios  /3.  Figure  5(a)  shows  that  the  signal- 
intensity  gain  decreases  as  the  optical  path  difference  of 
the  two  waves  increases  until  there  is  no  coupling  be¬ 
tween  the  two  waves  when  the  optical  path  difference  of 
the  two  waves  exceeds  the  coherent  length  of  the  source 
laser  wave.  The  same  figure  also  shows  that  the  signal- 
intensity  gain  increases  as  the  intensity  ratio  decreases 
until  the  signal-intensity  gain  saturates  as  the 
nondepleted-pump  regime  is  reached.  Figure  5(b)  shows 
that  the  normalized  mutual  coherence  of  the  two  waves  at 
the  pump-wave  entrance  plane  (z  =  L)  decreases  as  the 
intensity  ratio  P  of  the  two  waves  increases.  When  p  is 
much  larger  than  one,  the  signal  wave  will  pass  through 
the  photorefractive  medium  with  its  statistical  properties 
almost  unchanged  by  coupling.  In  this  limit  the  normal- 


(a) 


(b) 

Fig.  5.  (a)  Signal-intensity  gain  and  (b)  the  normalized  mutual 

coherence  of  the  two  waves  at  the  pump-wave  entrance  plane 
( z  =  L)  as  functions  of  the  optical  path  difference  at  the  signal- 
wave  entrance  plane  (2  =  0)  for  a  coupling  constant  y 
=  3  cm“^  and  various  intensity  ratios  between  the  signal  wave 
and  the  pump  wave.  Note  that  the  curve  for  a  coupling  constant 
y  -  3  cm~^  and  ?>  1  is  the  same  as  that  for  a  coupling  constant 
y  0  cm"^  and  an  arbitrary  p. 


ized  mutual  coherence  of  the  two  waves  at  the  pump-wave 
entrance  plane  (z  -  L)  with  coupling  wiD  be  the  same  as 
that  for  an  arbitrary  p  but  without  coupling.  We  note 
that  the  normalized  mutual  coherence  of  the  two  waves  at 
the  pump-wave  entrance  plane  (z  =  L)  is  enhanced  by 
coupling.  Figure  5(b)  also  shows  that  the  normalized 
mutual  coherence  decreases  quickly  as  the  optical  path 
difference  gets  close  to  and  larger  than  the  coherence 
length  of  the  source  laser  wave. 

Second,  we  consider  the  case  that  the  coherence  length 
of  the  source  laser  wave  is  shorter  than  the  thickness  of 
the  photorefractive  medium.  The  following  parameters 
are  chosen  in  the  simulation:  a  =  0.0  cm"^,  n  =  2.3, 
L  —  2.0  cm,  =  277  X  18  GHz,  and  -  -nLic. 

=  -nL/c  implies  that  the  optical  path  difference  be¬ 
tween  the  signal  wave  and  the  pump  wave  is  zero  at  the 
signal-wave  entrance  plane  (z  =  0).  Since  the  spectral 
line  shape  of  the  source  laser  wave  is  assumed  to  be 
Gaussian,  the  coherence  length  of  the  source  laser  wave 
inside  the  photorefractive  medium  is  =  2  v 
X  0.6QAc/{nS(i))  =  0.48  cm.  In  Fig.  6  we  show  the  sig¬ 
nal  intensity  /i(z)  and  the  grating  profile  Q(z)//o(2)  as 
functions  of  the  position  z  inside  the  photorefractive  me¬ 
dium  for  an  incident  intensity  ratio  p  -  1  and  a  coupling 
constant  y  =  20cm“^.  We  note  that  the  length  of  the 
photorefractive  grating  is  limited  by  the  partial  coherence 
of  the  two  interacting  waves  and  that  the  pump  depletion 
is  moderate  even  for  a  very  large  coupling  constant.  In 
Figs.  7(a)  and  7(b)  we  show  the  signal  intensity  Ii(z)  and 
the  grating  profile  Q(z)//o(2)  as  functions  of  the  position 
z  inside  the  photorefractive  medium  for  an  incident  inten¬ 
sity  ratio  p  =  10““*  and  coupling  constants  y  =  10  cm“^ 
and  y  =  20  cm"\  respectively.  We  note  that  the  length 
of  the  photorefractive  grating  is  increased  but  still  prima¬ 
rily  limited  by  the  partial  coherence  of  the  two  interacting 
waves  for  a  small  incident  intensity  ratio  and  a  small  cou¬ 
pling  constant.  Figure  7(b)  shows  that  the  length  of  the 
photorefractive  grating  is  no  longer  limited  by  the  partial 
coherence  of  the  two  interacting  waves  for  a  small  inci¬ 
dent  intensity  ratio  and  a  large  coupling  constant.  The 
length  of  the  photorefractive  grating  in  this  case  is  lim¬ 
ited  by  the  length  of  the  photorefractive  medium,  which  is 
much  longer  than  the  coherence  length  of  the  interacting 
waves.  The  photorefractive  grating  is  a  temporally  sta¬ 
tionary  index  grating.  When  the  incident  intensity  ratio 
is  small,  the  amplified  signal  wave  at  its  exit  plane  (z 
=  L)  is  primarily  the  incident  pump  wave  reflected  by 
the  photorefractive  grating,  and  the  depleted  pump  wave 
at  its  exit  plsuie  (z  =  0)  is  primarily  the  incident  piimp 
wave  transmitted  through  the  photorefractive  grating. 
When  the  length  of  the  photorefractive  grating  is  compa¬ 
rable  to  or  longer  than  the  coherence  length  of  the  source 
laser  wave  inside  the  photorefractive  medium,  the  output 
waves  Ei(L)  and  £2(0)  may  have  spectra  different  from 
those  of  the  input  waves  Ei(0)  and  E2(L)  owing  to  the 
presence  of  the  photorefractive  grating,  which  acts  as  a 
spectral  filter.  In  Figs.  8(a)  and  8(b)  we  show  the  nor¬ 
malized  spectra  of  the  amplified  signal  wave  and  the  de¬ 
pleted  pump  wave  at  their  respective  exit  planes,  using 
the  same  sets  of  parameters  as  those  used  to  obtain  Figs. 
7(a)  and  7(b),  respectively.  We  note  that  the  amplified 
signal  wave  has  a  bandwidth  narrower  than  that  of  the 
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.z(cm) 

Fig.  6.  Signal  intensity  Iiiz)  (solid  curve)  and  the  grating  pro¬ 
file  Q(2)//o(z)  (dashed  curve)  as  functions  of  the  position  z  inside 
the  photorefractive  medium  for  an  incident  intensity  ratio  p 
=  1  and  a  coupling  constant  y  =  20  cm'^. 


z(cm) 

(a) 


(b) 

Fig.  7.  Signal  intensity  /i(z)  (solid  curve)  and  the  grating  pro¬ 
file  Q(z)//o(z)  (dashed  curve)  as  functions  of  the  position  z  inside 
the  photorefractive  medium  for  an  incident  intensity  ratio  p 
=  10"^  and  coupling  constants  (a)  y  =  10  cm“^  and  (b)  y 
*=  20cm’^ 

incident  pump  wave.  Furthermore,  the  bandwidth  of  the 
amplified  signal  wave  decreases  as  the  length  of  the  pho¬ 
torefractive  grating  increases.  In  Fig.  8(a),  since  the  cou¬ 
pling  constant  is  small,  the  pump-wave  depletion  is  small 
(y  =  10  cm“^)  and  the  spectrum  of  the  transmitted  pump 
wave  is  almost  the  same  as  that  of  the  incident  pump 
wave.  In  Fig.  8(b),  since  the  coupling  constant  is  large 
(y  =  20  cm”^),  the  pump  depletion  is  significant.  The 
central  part  of  the  incident  pump-wave  spectrum  is  de¬ 
pleted  most  significantly,  and  therefore  the  transmitted 


pump  wave  has  a  different  spectrum  from  the  incident 
pump  wave.  According  to  the  spectra  shown  in  Fig.  8, 
the  case  considered  in  Fig.  7(a)  and  Fig.  8(a)  is  in  the 
nondepleted-pump  regime,  and  the  case  considered  in 
Fig.  7(b)  and  Fig.  8(b)  is  in  the  depleted-pump  regime,  al¬ 
though  both  cases  seem  to  be  in  the  nondepleted-pump  re¬ 
gime  when  we  look  only  at  the  signal  intensity.  In  Figs. 
9(a)  and  9(b)  we  show  with  solid  curves  the  normalized 
mutual  coherence  of  the  two  waves  as  a  function  of  posi¬ 
tion  z ,  again  using  the  same  sets  of  parameters  as  those 
used  to  obtain  Figs.  7(a)  and  7(b),  respectively.  For  com¬ 
parison  we  show  in  these  two  figures  with  dashed  curves 
the  normalized  mutual  coherence  of  the  two  waves  as  a 
function  of  position  z  without  coupling.  In  Fig.  9(a)  the 
normalized  mutual  coherence  of  the  two  waves  is  the 
same  with  coupling  as  that  without  coupling  in  the  region 
z  <  0  and  is  increased  in  the  region  z  >  L  owing  to  cou¬ 
pling;  in  Fig.  9(b)  the  normalized  mutual  coherence  of  the 
two  waves  is  decreased  in  the  region  close  to  plane 
z  =  0.  The  increase  of  mutual  coherence  in  both  Figs. 
9(a)  and  9(b)  can  be  attributed  to  the  reflection  of  the 
strong  pump  wave  in  the  direction  of  the  weak  signal 
wave  by  the  stationary  photorefractive  index  grating. 
The  decrease  of  mutual  coherence  in  Fig.  9(b)  can  be  at¬ 
tributed  to  the  spectral-filtering  effect  of  the  photorefrac¬ 
tive  index  grating  on  the  pump  wave.  There  is  no  de¬ 
crease  of  mutual  coherence  in  the  region  close  to  plane  z 
=  0  in  Fig.  9(a)  because  pump  depletion  and  therefore 
the  spectral-filtering  effect  on  the  pump  wave  are  negli¬ 
gible  in  this  case. 

Third,  when  a  laser  beam  enters  a  photorefractive  crys¬ 
tal,  scattering  occurs  because  of  surface  pits,  imperfec- 


Fig.  8.  Normalized  spectra  of  the  amplified  signal  wave  (solid 
curve)  and  the  depleted  pump  wave  (dashed  curve)  at  their  re¬ 
spective  exit  planes  for  an  incident  intensity  ratio  p  =  10““*  and 
coupling  constants  (a)  y  =  10  cm“^  and  (b)  y  =  20  cm“^. 
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the  coherence  length  of  the  incident  beams  for  various 
value  of  the  coupling  constant  y.  The  parameters  in  this 
simulation  are  n  =  2.3,  L  =  0.72  cm,  a  =  0.0  cm’\  and 
=  1.0  X  10“®.  The  time  delay  of  the  two  waves  is  as¬ 
sumed  to  be  zero  at  the  signal-wave  entrance  plane 
(z  =  0).  The  coherence  length  of  the  incident  beams  is 
related  to  the  bandwidth  as  =  277  x  0.664c/ Scj.  We 
note  that  the  phase-conjugation  reflectivity  increases  as 
the  coherence  length  of  the  incident  beam  increases  and 
reaches  to  a  constant  when  the  coherence  length  of  the  in¬ 
cident  beam  is  much  longer  than  the  thickness  of  the  pho- 
torefractive  medium.  We  also  note  that  the  phase- 
conjugation  reflectivity  increases  as  the  coupling  constant 
increases. 


position  z  (cm) 

(b) 

Fig.  9.'  Normalized  mutual  coherence  of  the  two  waves  as  func¬ 
tions  of  position  z  for  an  incident  intensity  ratio  ^  =  10“^  and  for 
coupling  constants  (a)  y  =  10  cm^^  (solid  curve)  and  y 
-  0  cm‘^  (dashed  curv'e)  and  (b)  y  =  20  cm~^  (solid  curve)  and 
y  =  0  cm“^  (dashed  curve). 
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Fig.  10.  Phase-conjugation  reflectivity  as  a  function  of  the  co¬ 
herence  length  of  the  incident  beam  for  various  value  of  the  cou¬ 
pling  constant*  y. 

tion,  and  defects  in  the  crystal.  The  scattered  light  over¬ 
laps  with  the  incident  beam  and  they  undergo  two-wave 
mixing.  Under  the  appropriate  condition  the  scattered 
light  can  be  amplified,  leading  to  phenomena  such  as  fan¬ 
ning  and  stimulated  backscattering.  In  a  manner  very 
similar  to  stimulated  Brillouin  scattering  the  stimulated 
backward  scattering  in  photorefractive  media  is  a  pos¬ 
sible  mechanism  for  self-pumped  phase  conjugation. 

The  general  formulation  developed  in  this  paper  can  be 
employed  to  investigate  the  effect  of  coherence  on  self- 
pumped  phase  conjugation  by  2k  gratings.  We  show  in 
Fig.  10  the  phase-conjugation  reflectivity  as  a  function  of 


4.  EXPERIMENTS 

The  above  theory  is  validated  experimentally.  The  ex¬ 
perimental  setup  is  shown  in  Fig.  11.  We  utilized  a  45®- 
cut  KNbOsiCo  crystal  (the  c  axis  is  in  the  horizontal 
plane  leaning  toward  the  2  =  0  face  of  the  crystal,  and 
the  b  axis  is  in  the  vertical  direction).  The  measured  pa¬ 
rameters  of  the  crystal  are  y  =  3.3  cm‘^,  a  =  0.5  cm"\ 
n  =  2.3,  and  L  -  0.72  cm.  A  multimode  argon  laser  op¬ 
erating  at  514  nm  is  used  as  the  laser  source  with  a  mea¬ 
sured  FWHM  bandwidth  of  1.83  GHz.  The  extraordi¬ 
nary  polarization  of  the  laser  wave  is  used  in  the 
experiment.  As  illustrated  in  Fig.  11,  the  signal  wave 
and  the  pump  wave,  obtained  by  splitting  the  argon  laser 
wave,  propagate  contradirectionally  into  the  KNb03:Co 
crystal.  The  incident  intensity  ratio  of  the  two  waves  is 
p  =  0.00151.  The  power  of  the  pump  wave  is  main¬ 
tained  at  -“50  mW.  The  optical  path  difference  of  the 
two  waves  at  the  signal-wave  incident  plane  2  =  0  is  de¬ 
noted  as  AL  “  L2  “  Lj .  To  monitor  the  mutual  coher¬ 
ence  between  the  signal  wave  and  the  pump  wave  at  the 
output  plane  2  =  L,  we  employed  another  reference  wave 
(E2ref)  l^hat  was  split  from  the  pump  wave  E2 .  The  opti¬ 
cal  path  difference  of  Ej  and  E2ref  waves  was  adjusted  to 
be  the  same  as  that  of  E^  and  E2  waves  at  the  output 
plane  2  =  L.  By  a  simple  homodyne  technique,  the  in¬ 


laser  source  to  the  signal-wave  incident  plane  2  =  0,  respec¬ 
tively.  Loref  is  the  optical  path  length  of  reference  wave  F^ref 
from  the  laser  source  to  the  signal  output  plane  2  =  Z/. 
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(a)  (b) 

Fig.  12.  Interference  pattern  of  the  signal  wave  and  the  refer¬ 
ence  wave  at  the  output  plane  PI  (a)  without  photorefractive  cou¬ 
pling  and  (b)  with  photorefractive  coupling.  Note  the  increase  of 
fringe  visibility  that  is  due  to  the  coupling. 
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optical  path  difference  (cm) 

Fig.  13.  Signal-intensity  gain  as  a  function  of  the  optical  path 
difference  of  the  two  waves  at  the  signal-wave  entrance  plane 
(z  =  0).  The  dots  are  experimental  data,  and  the  solid  curve  is 
the  theoretical  data. 

terference  fringes  generated  by  Ei  and  E2ref  waves  were 
observed  by  a  CCD  camera  at  the  output  plane  PI.  The 
normalized  mutual  coherence  y^iL)  can  be  estimated  as 
(^max  -  where  is  the  ampli¬ 

tude  of  the  fringes.  In  our  experiment  we  monitored  the 
interference  pattern  with  and  without  pump  beam  E2. 
Figure  12  shows  the  interference  patterns  with  a  normal¬ 
ized  mutual  coherence  ri2(0,  0)  0.43  at  2  =  0  (AL 

=  4  cm).  The  measured  normalized  mutual  coherence 
increases  from  0.19  to  0.7  at  2  =  L.  The  intensity  gain 
of  the  signal  wave  was  also  measured.  Figure  13  shows 
the  measurement  of  the  intensity  gain  (dots)  of  the  signal 
wave  at  the  2  =  L  plane  as  a  function  of  the  optical  path 
difference  AL.  Along  with  the  data  is  the  theoretical 
curve  for  the  same  parameters.  An  excellent  agreement 
between  theory  and  experiment  was  achieved. 

5.  CONCLUSIONS 

In  conclusion,  we  have  investigated  theoretically  contra- 
directional  two-wave  mixing  with  partially  coherent 
waves  in  photorefractive  crystals.  A  general  formulation 
based  on  the  theory  of  statistical  properties  of  linear  sys¬ 
tems  is  provided.^®  Previous  results  on  several  simpli¬ 
fied  cases  are  rederived  as  special  cases  of  the  general  for¬ 
mulation  so  that  we  can  get  more  insight  into  the  general 
formulation  as  well  as  into  the  simplified  cases.  Results 
of  numerical  implementation  of  the  general  formulation 


are  also  provided  for  various  coupling  constants  and  vari¬ 
ous  incident  intensity  ratios  between  the  signal  wave  and 
the  pump  wave.  We  found  that  the  mutual  coherence  be¬ 
tween  the  signal  wave  and  the  pump  wave  can  be  both  in¬ 
creased  and  decreased  owing  to  coupling.  We  also  found 
that  both  the  strength  and  the  length  of  the 
photorefractive-index  grating  increases  as  the  coupling 
constant  increases.  A  photorefractive  index  grating 
much  longer  than  the  coherence  length  of  the  incident 
waves  can  be  formed  when  the  coupling  constant  is  large 
and  the  incident  intensity  ratio  is  small.  Owing  to  the 
spectral-filtering  effect  of  the  photorefractive-index  grat¬ 
ing,  the  spectra  of  the  two  interacting  waves  will  be  al¬ 
tered  as  they  pass  through  the  photorefractive  medium. 
The  general  formulation  is  further  used  to  model  the  ef¬ 
fect  of  partial  coherence  on  self-pumped  phase  conjuga¬ 
tion  by  a  2k  grating.  The  theoretical  predictions  are  in 
excellent  agreement  with  experimental  measurements. 

APPENDIX  A:  ERGODICITY  PROPERTY 

In  this  appendix  we  provide  the  derivation  of  Eq.  (10). 
Equation  (9)  gives  the  dynamics  of  the  photorefractive  in¬ 
dex  grating.  It  can  be  rewritten  into  an  integration  form 
as 

Q(2,  0  =  —  f  Ei(2,  r)E2*(2,  r)exp(^ 

^ph  J  \  ^ph 

(Al) 

Note  that  l/rp^Z-x  exp[(i'  -  O/rphJd^'  =  1.  According  to 
Eq.  (Al),  the  grating  amplitude  Q(2,  t)  is  approximately 
the  average  value  of  Ei(2,  ^')E2*(2,  t')  over  a  time  pe¬ 
riod  Tph  . 

Since  the  optical  wave  amplitudes  £^(2,  t)  and 
£2(2,  t)  are  stationary  random  processes,  the  ensemble 
average  {£1(2,  OE2*(2,  t))  is  independent  of  the  time 
variable  t.  T-aking  the  ensemble  average  of  Eq.  (Al),  we 
can  obtain 

(Q(2,  t))  =  (£1(2,  0E2*(2,  t))  ^  ri2(2,  0).  (A2) 

Equation  (A2)  shows  that  the  ensemble  average  of  the 
grating  amplitude  Q(2,  t)  is  equal  to  the  mutual  coher¬ 
ence  of  the  two  waves. 

In  general,  the  grating  amplitude  Q(2,  t)  is  also  a  ran¬ 
dom  variable.  It  fluctuates  around  its  ensemble  average. 
The  mean  square  value  of  this  random  fluctuation,  also 
called  the  variance  of  the  random  vEiriable  Q(2,  0,  can  be 
written  as 

<IQ(e,  t)  -  <Q(2,  tm  =  (1Q(2,  t)\^)  -  KQ(^, 

(A3) 

By  using  Eq.  (Al),  we  can  write  the  first  term  on  the  right 
side  of  Eq.  (A3)  explicitly  as 

(IQ(^,  t)\^)  =  A  f 

’■ph 

It'-t]  lt"-t\ 

X  exp  -  exp  -  df'di", 

\  ^ph  /  \  ^ph  / 


(A4) 
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wlierxj  F(z,  t)  ^  EiCr,  t)F,2*{z,  t)  is  a  shorthand  nota¬ 
tion.  Since  the  optical  wave  amplitudes  Ei(z,  t)  and 
E2(z,  t)  are  stationary  random  processes,  F(2,  t)  is  also  a 
stationary  random  process.  Therefore  the  ensemble  av¬ 
erage  (F(2,  ^')F*(2,  t'"))  is  a  function  of  only  two  vari¬ 
ables  z  and  t*  -  t”.  We  define  R(2,  t’  - 
s  {F(z,  ^')F*(2,  ^"))  as  a  shorthand  notation.  With 
this  relation  we  can  change  the  integration  arguments  in 
Eq.  (A4)  from  and  t"  to  ti  —  t"  -  2t  and  t2  —  t* 
-  t"  and  rewrite  Eq.  (A4)  as 


Ei(2,  0E2*(2,  t).  This  type  of  property  is  referred  to  as 
ergodicity  in  statistics. 
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(A5) 


Note  that  l/(2rph^)/®ocd^i  exp(^i/Tph)/^~^^d^2  =  1*  Sub¬ 
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